Genomic epidemiology of community-associated Staphylococcus aureus in northern Australia and Papua New Guinea

Eike Steinig investigated the transmission dynamics of community-associated MRSA in Far North Queensland and Papua New Guinea. He developed bioinformatics methods based on portable nanopore sequencing technology to infer epidemiological outbreak parameters and clinical genotypes of bacterial pathogens, with applications in rural and remote hospitals that lack sequencing infrastructure

Taking the next-gen step: Comprehensive antimicrobial resistance detection from Burkholderia pseudomallei

Background: Antimicrobial resistance (AMR) poses a major threat to human health. Whole-genome sequencing holds great potential for AMR identification; however, there remain major gaps in accurately and comprehensively detecting AMR across the spectrum of AMR-conferring determinants and pathogens. Methods: Using 16 wild-type Burkholderia pseudomallei and 25 with acquired AMR, we first assessed the performance of existing AMR software (ARIBA, CARD, ResFinder, and AMRFinderPlus) for detecting clinically relevant AMR in this pathogen. B. pseudomallei was chosen due to limited treatment options, high fatality rate, and AMR caused exclusively by chromosomal mutation (i.e. single-nucleotide polymorphisms [SNPs], insertions-deletions [indels], copy-number variations [CNVs], inversions, and functional gene loss). Due to poor performance with existing tools, we developed ARDaP (Antimicrobial Resistance Detection and Prediction) to identify the spectrum of AMR-conferring determinants in B. pseudomallei. Findings: CARD, ResFinder, and AMRFinderPlus failed to identify any clinically-relevant AMR in B. pseudomallei; ARIBA identified AMR encoded by SNPs and indels that were manually added to its database. However, none of these tools identified CNV, inversion, or gene loss determinants, and ARIBA could not differentiate AMR determinants from natural genetic variation. In contrast, ARDaP accurately detected all SNP, indel, CNV, inversion, and gene loss AMR determinants described in B. pseudomallei (n≈50). Additionally, ARDaP accurately predicted three previously undescribed determinants. In mixed strain data, ARDaP identified AMR to as low as ~5% allelic frequency. Interpretation: Existing AMR software packages are inadequate for chromosomal AMR detection due to an inability to detect resistance conferred by CNVs, inversions, and functional gene loss. ARDaP overcomes these major shortcomings. Further, ARDaP enables AMR prediction from mixed sequence data down to 5% allelic frequency, and can differentiate natural genetic variation from AMR determinants. ARDaP databases can be constructed for any microbial species of interest for comprehensive AMR detection

Staphylococcus aureus from patients with chronic rhinosinusitis show minimal genetic association between polyp and non-polyp phenotypes

Background: Staphylococcus aureus has a high prevalence in chronic rhinosinusitis (CRS) patients and is suggested to play a more etiopathogenic role in CRS patients with nasal polyps (CRSwNP), a severe form of the CRS spectrum with poorer surgical outcomes. We performed a microbial genome-wide association study (mGWAS) to investigate whether S. aureus isolates from CRS patients have particular genetic markers associated with CRS with nasal polyps (CRSwNP) or CRS without nasal polyps (CRSsNP). Methods: Whole genome sequencing was performed on S. aureus isolates collected from 28 CRSsNP and 30 CRSwNP patients. A mGWAS approach was employed using large-scale comparative genomics to identify genetic variation within our dataset. Results: Considerable genetic variation was observed, with >90,000 single nucleotide polymorphisms (SNPs) sites identified. There was little correlation with CRS subtype based on SNPs and Insertion/Delection (Indels). One indel was found to significantly correlate with CRSwNP and occurred in the promoter region of a bacitracin transport system ATP-binding protein. Additionally, two variants of the highly variable superantigen-like (SSL) proteins were found to significantly correlate with each CRS phenotype. No significant association with other virulence or antibiotic resistance genes were observed, consistent with previous studies. Conclusion: To our knowledge this study is the first to use mGWAS to investigate the contribution of microbial genetic variation to CRS presentations. Utilising the most comprehensive genome-wide analysis methods available, our results suggest that CRS phenotype may be influenced by genetic factors other than specific virulence mechanisms within the S. aureus genome

Global Scale Dissemination of ST93: A Divergent Staphylococcus aureus Epidemic Lineage That Has Recently Emerged From Remote Northern Australia.

Background: In Australia, community-associated methicillin-resistant Staphylococcus aureus (MRSA) lineage sequence type (ST) 93 has rapidly risen to dominance since being described in the early 1990s. We examined 459 ST93 genome sequences from Australia, New Zealand, Samoa, and Europe to investigate the evolutionary history of ST93, its emergence in Australia and subsequent spread overseas. Results: Comparisons with other S. aureus genomes indicate that ST93 is an early diverging and recombinant lineage, comprising of segments from the ST59/ST121 lineage and from a divergent but currently unsampled Staphylococcal population. However, within extant ST93 strains limited genetic diversity was apparent with the most recent common ancestor dated to 1977 (95% highest posterior density 1973-1981). An epidemic ST93 population arose from a methicillin-susceptible progenitor in remote Northern Australia, which has a proportionally large Indigenous population, with documented overcrowded housing and a high burden of skin infection. Methicillin-resistance was acquired three times in these regions, with a clade harboring a staphylococcal cassette chromosome mec (SCCmec) IVa expanding and spreading to Australia's east coast by 2000. We observed sporadic and non-sustained introductions of ST93-MRSA-IVa to the United Kingdom. In contrast, in New Zealand, ST93-MRSA-IVa was sustainably transmitted with clonal expansion within the Pacific Islander population, who experience similar disadvantages as Australian Indigenous populations. Conclusion: ST93 has a highly recombinant genome including portions derived from an early diverging S. aureus population. Our findings highlight the need to understand host population factors in the emergence and spread of antimicrobial resistant community pathogens

Phylodynamic Inference of Bacterial Outbreak Parameters Using Nanopore Sequencing

Nanopore sequencing and phylodynamic modeling have been used to reconstruct the transmission dynamics of viral epidemics, but their application to bacterial pathogens has remained challenging. Cost-effective bacterial genome sequencing and variant calling on nanopore platforms would greatly enhance surveillance and outbreak response in communities without access to sequencing infrastructure. Here, we adapt random forest models for single nucleotide polymorphism (SNP) polishing developed by Sanderson and colleagues (2020. High precision Neisseria gonorrhoeae variant and antimicrobial resistance calling from metagenomic nanopore sequencing. Genome Res. 30(9):1354–1363) to estimate divergence and effective reproduction numbers (Re) of two methicillin-resistant Staphylococcus aureus (MRSA) outbreaks from remote communities in Far North Queensland and Papua New Guinea (PNG; n = 159). Successive barcoded panels of S. aureus isolates (2 × 12 per MinION) sequenced at low coverage (>5× to 10×) provided sufficient data to accurately infer genotypes with high recall when compared with Illumina references. Random forest models achieved high resolution on ST93 outbreak sequence types (>90% accuracy and precision) and enabled phylodynamic inference of epidemiological parameters using birth–death skyline models. Our method reproduced phylogenetic topology, origin of the outbreaks, and indications of epidemic growth (Re > 1). Nextflow pipelines implement SNP polisher training, evaluation, and outbreak alignments, enabling reconstruction of within-lineage transmission dynamics for infection control of bacterial disease outbreaks on portable nanopore platforms. Our study shows that nanopore technology can be used for bacterial outbreak reconstruction at competitive costs, providing opportunities for infection control in hospitals and communities without access to sequencing infrastructure, such as in remote northern Australia and PNG

Exploring the evolution and epidemiology of European CC1-MRSA-IV: tracking a multidrug-resistant community-associated meticillin-resistant Staphylococcus aureus clone

This study investigated the evolution and epidemiology of the community-associated and multidrug-resistant Staphylococcus aureus clone European CC1-MRSA-IV. Whole-genome sequences were obtained for 194 European CC1-MRSA-IV isolates (189 of human and 5 of animal origin) from 12 countries, and 10 meticillin-susceptible precursors (from North-Eastern Romania; all of human origin) of the clone. Phylogenetic analysis was performed using a maximum-likelihood approach, a time-measured phylogeny was reconstructed using Bayesian analysis, and in silico microarray genotyping was performed to identify resistance, virulence-associated and SCCmec (staphylococcal cassette chromosome mec) genes. Isolates were typically sequence type 1 (190/204) and spa type t127 (183/204). Bayesian analysis indicated that European CC1-MRSA-IV emerged in approximately 1995 before undergoing rapid expansion in the late 1990s and 2000s, while spreading throughout Europe and into the Middle East. Phylogenetic analysis revealed an unstructured meticillin-resistant S. aureus (MRSA) population, lacking significant geographical or temporal clusters. The MRSA were genotypically multidrug-resistant, consistently encoded seh, and intermittently (34/194) encoded an undisrupted hlb gene with concomitant absence of the lysogenic phage-encoded genes sak and scn. All MRSA also harboured a characteristic ~5350 nt insertion in SCCmec adjacent to orfX. Detailed demographic data from Denmark showed that there, the clone is typically (25/35) found in the community, and often (10/35) among individuals with links to South-Eastern Europe. This study elucidated the evolution and epidemiology of European CC1-MRSA-IV, which emerged from a meticillin-susceptible lineage prevalent in North-Eastern Romania before disseminating rapidly throughout Europe

Phylodynamic signatures in the emergence of community-associated MRSA

Community-associated, methicillin-resistant Staphylococcus aureus (MRSA) lineages have emerged in many geographically distinct regions around the world during the past 30 y. Here, we apply consistent phylodynamic methods across multiple community-associated MRSA lineages to describe and contrast their patterns of emergence and dissemination. We generated whole-genome sequencing data for the Australian sequence type (ST) ST93-MRSA-IV from remote communities in Far North Queensland and Papua New Guinea, and the Bengal Bay ST772-MRSA-V clone from metropolitan communities in Pakistan. Increases in the effective reproduction number (Re) and sustained transmission (Re > 1) coincided with spread of progenitor methicillin-susceptible S. aureus (MSSA) in remote northern Australian populations, dissemination of the ST93-MRSA-IV genotype into population centers on the Australian East Coast, and subsequent importation into the highlands of Papua New Guinea and Far North Queensland. Applying the same phylodynamic methods to existing lineage datasets, we identified common signatures of epidemic growth in the emergence and epidemiological trajectory of community-associated S. aureus lineages from America, Asia, Australasia, and Europe. Surges in Re were observed at the divergence of antibiotic-resistant strains, coinciding with their establishment in regional population centers. Epidemic growth was also observed among drug-resistant MSSA clades in Africa and northern Australia. Our data suggest that the emergence of community-associated MRSA in the late 20th century was driven by a combination of antibiotic-resistant genotypes and host epidemiology, leading to abrupt changes in lineage-wide transmission dynamics and sustained transmission in regional population centers

Phylodynamic signatures in the emergence of community-associated MRSA

Community-associated, methicillin-resistant Staphylococcus aureus (MRSA) lineages have emerged in many geographically distinct regions around the world during the past 30 y. Here, we apply consistent phylodynamic methods across multiple community-associated MRSA lineages to describe and contrast their patterns of emergence and dissemination. We generated whole-genome sequencing data for the Australian sequence type (ST) ST93-MRSA-IV from remote communities in Far North Queensland and Papua New Guinea, and the Bengal Bay ST772-MRSA-V clone from metropolitan communities in Pakistan. Increases in the effective reproduction number (Re) and sustained transmission (Re > 1) coincided with spread of progenitor methicillin-susceptible S. aureus (MSSA) in remote northern Australian populations, dissemination of the ST93-MRSA-IV genotype into population centers on the Australian East Coast, and subsequent importation into the highlands of Papua New Guinea and Far North Queensland. Applying the same phylodynamic methods to existing lineage datasets, we identified common signatures of epidemic growth in the emergence and epidemiological trajectory of community-associated S. aureus lineages from America, Asia, Australasia, and Europe. Surges in Re were observed at the divergence of antibiotic-resistant strains, coinciding with their establishment in regional population centers. Epidemic growth was also observed among drug-resistant MSSA clades in Africa and northern Australia. Our data suggest that the emergence of community-associated MRSA in the late 20th century was driven by a combination of antibiotic-resistant genotypes and host epidemiology, leading to abrupt changes in lineage-wide transmission dynamics and sustained transmission in regional population centers

Evolution and Global Transmission of a Multidrug-Resistant, Community-Associated Methicillin-Resistant Staphylococcus aureus Lineage from the Indian Subcontinent

The evolution and global transmission of antimicrobial resistance have been well documented for Gram-negative bacteria and health care-associated epidemic pathogens, often emerging from regions with heavy antimicrobial use. However, the degree to which similar processes occur with Gram-positive bacteria in the community setting is less well understood. In this study, we traced the recent origins and global spread of a multidrug-resistant, community-associated Staphylococcus aureus lineage from the Indian subcontinent, the Bengal Bay clone (ST772). We generated whole-genome sequence data of 340 isolates from 14 countries, including the first isolates from Bangladesh and India, to reconstruct the evolutionary history and genomic epidemiology of the lineage. Our data show that the clone emerged on the Indian subcontinent in the early 1960s and disseminated rapidly in the 1990s. Short-term outbreaks in community and health care settings occurred following intercontinental transmission, typically associated with travel and family contacts on the subcontinent, but ongoing endemic transmission was uncommon. Acquisition of a multidrug resistance integrated plasmid was instrumental in the emergence of a single dominant and globally disseminated clade in the early 1990s. Phenotypic data on biofilm, growth, and toxicity point to antimicrobial resistance as the driving force in the evolution of ST772. The Bengal Bay clone therefore combines the multidrug resistance of traditional health care-associated clones with the epidemiological transmission of community-associated methicillin-resistant S. aureus (MRSA). Our study demonstrates the importance of whole-genome sequencing for tracking the evolution of emerging and resistant pathogens. It provides a critical framework for ongoing surveillance of the clone on the Indian subcontinent and elsewhere

Evolution and Global Transmission of a Multidrug-Resistant, Community-Associated Methicillin-Resistant Staphylococcus aureus Lineage from the Indian Subcontinent.

The evolution and global transmission of antimicrobial resistance have been well documented for Gram-negative bacteria and health care-associated epidemic pathogens, often emerging from regions with heavy antimicrobial use. However, the degree to which similar processes occur with Gram-positive bacteria in the community setting is less well understood. In this study, we traced the recent origins and global spread of a multidrug-resistant, community-associated Staphylococcus aureus lineage from the Indian subcontinent, the Bengal Bay clone (ST772). We generated whole-genome sequence data of 340 isolates from 14 countries, including the first isolates from Bangladesh and India, to reconstruct the evolutionary history and genomic epidemiology of the lineage. Our data show that the clone emerged on the Indian subcontinent in the early 1960s and disseminated rapidly in the 1990s. Short-term outbreaks in community and health care settings occurred following intercontinental transmission, typically associated with travel and family contacts on the subcontinent, but ongoing endemic transmission was uncommon. Acquisition of a multidrug resistance integrated plasmid was instrumental in the emergence of a single dominant and globally disseminated clade in the early 1990s. Phenotypic data on biofilm, growth, and toxicity point to antimicrobial resistance as the driving force in the evolution of ST772. The Bengal Bay clone therefore combines the multidrug resistance of traditional health care-associated clones with the epidemiological transmission of community-associated methicillin-resistant S. aureus (MRSA). Our study demonstrates the importance of whole-genome sequencing for tracking the evolution of emerging and resistant pathogens. It provides a critical framework for ongoing surveillance of the clone on the Indian subcontinent and elsewhere.IMPORTANCE The Bengal Bay clone (ST772) is a community-associated and multidrug-resistant Staphylococcus aureus lineage first isolated from Bangladesh and India in 2004. In this study, we showed that the Bengal Bay clone emerged from a virulent progenitor circulating on the Indian subcontinent. Its subsequent global transmission was associated with travel or family contact in the region. ST772 progressively acquired specific resistance elements at limited cost to its fitness and continues to be exported globally, resulting in small-scale community and health care outbreaks. The Bengal Bay clone therefore combines the virulence potential and epidemiology of community-associated clones with the multidrug resistance of health care-associated S. aureus lineages. This study demonstrates the importance of whole-genome sequencing for the surveillance of highly antibiotic-resistant pathogens, which may emerge in the community setting of regions with poor antibiotic stewardship and rapidly spread into hospitals and communities across the world