High performance Legionella pneumophila source attribution using genomics-based machine learning classification

Fundamental to effective Legionnaires' disease outbreak control is the ability to rapidly identify the environmental source(s) of the causative agent, Legionella pneumophila. Genomics has revolutionized pathogen surveillance, but L. pneumophila has a complex ecology and population structure that can limit source inference based on standard core genome phylogenetics. Here, we present a powerful machine learning approach that assigns the geographical source of Legionnaires' disease outbreaks more accurately than current core genome comparisons. Models were developed upon 534 L. pneumophila genome sequences, including 149 genomes linked to 20 previously reported Legionnaires' disease outbreaks through detailed case investigations. Our classification models were developed in a cross-validation framework using only environmental L. pneumophila genomes. Assignments of clinical isolate geographic origins demonstrated high predictive sensitivity and specificity of the models, with no false positives or false negatives for 13 out of 20 outbreak groups, despite the presence of within-outbreak polyclonal population structure. Analysis of the same 534-genome panel with a conventional phylogenomic tree and a core genome multi-locus sequence type allelic distance-based classification approach revealed that our machine learning method had the highest overall classification performance-agreement with epidemiological information. Our multivariate statistical learning approach maximizes the use of genomic variation data and is thus well-suited for supporting Legionnaires' disease outbreak investigations.IMPORTANCEIdentifying the sources of Legionnaires' disease outbreaks is crucial for effective control. Current genomic methods, while useful, often fall short due to the complex ecology and population structure of Legionella pneumophila, the causative agent. Our study introduces a high-performing machine learning approach for more accurate geographical source attribution of Legionnaires' disease outbreaks. Developed using cross-validation on environmental L. pneumophila genomes, our models demonstrate excellent predictive sensitivity and specificity. Importantly, this new approach outperforms traditional methods like phylogenomic trees and core genome multi-locus sequence typing, proving more efficient at leveraging genomic variation data to infer outbreak sources. Our machine learning algorithms, harnessing both core and accessory genomic variation, offer significant promise in public health settings. By enabling rapid and precise source identification in Legionnaires' disease outbreaks, such approaches have the potential to expedite intervention efforts and curtail disease transmission.

Genomic surveillance for antimicrobial resistance - a One Health perspective

Antimicrobial resistance (AMR) - the ability of microorganisms to adapt and survive under diverse chemical selection pressures - is influenced by complex interactions between humans, companion and food-producing animals, wildlife, insects and the environment. To understand and manage the threat posed to health (human, animal, plant and environmental) and security (food and water security and biosecurity), a multifaceted 'One Health' approach to AMR surveillance is required. Genomic technologies have enabled monitoring of the mobilization, persistence and abundance of AMR genes and mutations within and between microbial populations. Their adoption has also allowed source-tracing of AMR pathogens and modelling of AMR evolution and transmission. Here, we highlight recent advances in genomic AMR surveillance and the relative strengths of different technologies for AMR surveillance and research. We showcase recent insights derived from One Health genomic surveillance and consider the challenges to broader adoption both in developed and in lower- and middle-income countries.

High-performance enrichment-based genome sequencing to support the investigation of hepatitis A virus outbreaks

IMPORTANCE

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.

Bioinformatic investigation of discordant sequence data for SARS-CoV-2: insights for robust genomic analysis during pandemic surveillance

The COVID-19 pandemic has necessitated the rapid development and implementation of whole-genome sequencing (WGS) and bioinformatic methods for managing the pandemic. However, variability in methods and capabilities between laboratories has posed challenges in ensuring data accuracy. A national working group comprising 18 laboratory scientists and bioinformaticians from Australia and New Zealand was formed to improve data concordance across public health laboratories (PHLs). One effort, presented in this study, sought to understand the impact of the methodology on consensus genome concordance and interpretation. SARS-CoV-2 WGS proficiency testing programme (PTP) data were retrospectively obtained from the 2021 Royal College of Pathologists of Australasia Quality Assurance Programmes (RCPAQAP), which included 11 participating Australian laboratories. The submitted consensus genomes and reads from eight contrived specimens were investigated, focusing on discordant sequence data and findings were presented to the working group to inform best practices. Despite using a variety of laboratory and bioinformatic methods for SARS-CoV-2 WGS, participants largely produced concordant genomes. Two participants returned five discordant sites in a high-Cτ replicate, which could be resolved with reasonable bioinformatic quality thresholds. We noted ten discrepancies in genome assessment that arose from nucleotide heterogeneity at three different sites in three cell-culture-derived control specimens. While these sites were ultimately accurate after considering the participants' bioinformatic parameters, it presented an interesting challenge for developing standards to account for intrahost single nucleotide variation (iSNV). Observed differences had little to no impact on key surveillance metrics, lineage assignment and phylogenetic clustering, while genome coverage <90 % affected both. We recommend PHLs bioinformatically generate two consensus genomes with and without ambiguity thresholds for quality control and downstream analysis, respectively, and adhere to a minimum 90 % genome coverage threshold for inclusion in surveillance interpretations. We also suggest additional PTP assessment criteria, including primer efficiency, detection of iSNVs and minimum genome coverage of 90 %. This study underscores the importance of multidisciplinary national working groups in informing guidelines in real time for bioinformatic quality acceptance criteria. It demonstrates the potential for enhancing public health responses through improved data concordance and quality control in SARS-CoV-2 genomic analysis during pandemic surveillance.

The two-component system WalKR provides an essential link between cell wall homeostasis and DNA replication in Staphylococcus aureus

Among the 16 two-component systems in the opportunistic human pathogen Staphylococcus aureus, only WalKR is essential. Like the orthologous systems in other Bacillota, S. aureus WalKR controls autolysins involved in peptidoglycan remodeling and is therefore intimately involved in cell division. However, despite the importance of WalKR in S. aureus, the basis for its essentiality is not understood and the regulon is poorly defined. Here, we defined a consensus WalR DNA-binding motif and the direct WalKR regulon by using functional genomics, including chromatin immunoprecipitation sequencing, with a panel of isogenic walKR mutants that had a spectrum of altered activities. Consistent with prior findings, the direct regulon includes multiple autolysin genes. However, this work also revealed that WalR directly regulates at least five essential genes involved in lipoteichoic acid synthesis (ltaS): translation (rplK), DNA compaction (hup), initiation of DNA replication (dnaA, hup) and purine nucleotide metabolism (prs). Thus, WalKR in S. aureus serves as a polyfunctional regulator that contributes to fundamental control over critical cell processes by coordinately linking cell wall homeostasis with purine biosynthesis, protein biosynthesis, and DNA replication. Our findings further address the essentiality of this locus and highlight the importance of WalKR as a bona fide target for novel anti-staphylococcal therapeutics. IMPORTANCE The opportunistic human pathogen Staphylococcus aureus uses an array of protein sensing systems called two-component systems (TCS) to sense environmental signals and adapt its physiology in response by regulating different genes. This sensory network is key to S. aureus versatility and success as a pathogen. Here, we reveal for the first time the full extent of the regulatory network of WalKR, the only staphylococcal TCS that is indispensable for survival under laboratory conditions. We found that WalKR is a master regulator of cell growth, coordinating the expression of genes from multiple, fundamental S. aureus cellular processes, including those involved in maintaining cell wall metabolism, protein biosynthesis, nucleotide metabolism, and the initiation of DNA replication.

A statistical genomics framework to trace bacterial genomic predictors of clinical outcomes in Staphylococcus aureus bacteremia

Outcomes of severe bacterial infections are determined by the interplay between host, pathogen, and treatments. While human genomics has provided insights into host factors impacting Staphylococcus aureus infections, comparatively little is known about S. aureus genotypes and disease severity. Building on the hypothesis that bacterial pathoadaptation is a key outcome driver, we developed a genome-wide association study (GWAS) framework to identify adaptive mutations associated with treatment failure and mortality in S. aureus bacteremia (1,358 episodes). Our research highlights the potential of vancomycin-selected mutations and vancomycin minimum inhibitory concentration (MIC) as key explanatory variables to predict infection severity. The contribution of bacterial variation was much lower for clinical outcomes (heritability <5%); however, GWASs allowed us to identify additional, MIC-independent candidate pathogenesis loci. Using supervised machine learning, we were able to quantify the predictive potential of these adaptive signatures. Our statistical genomics framework provides a powerful means to capture adaptive mutations impacting severe bacterial infections.

The importance of utilizing travel history metadata for informative phylogeographical inferences: a case study of early SARS-CoV-2 introductions into Australia

Inferring the spatiotemporal spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) via Bayesian phylogeography has been complicated by the overwhelming sampling bias present in the global genomic dataset. Previous work has demonstrated the utility of metadata in addressing this bias. Specifically, the inclusion of recent travel history of SARS-CoV-2-positive individuals into extended phylogeographical models has demonstrated increased accuracy of estimates, along with proposing alternative hypotheses that were not apparent using only genomic and geographical data. However, as the availability of comprehensive epidemiological metadata is limited, many of the current estimates rely on sequence data and basic metadata (i.e. sample date and location). As the bias within the SARS-CoV-2 sequence dataset is extensive, the degree to which we can rely on results drawn from standard phylogeographical models (i.e. discrete trait analysis) that lack integrated metadata is of great concern. This is particularly important when estimates influence and inform public health policy. We compared results generated from the same dataset, using two discrete phylogeographical models: one including travel history metadata and one without. We utilized sequences from Victoria, Australia, in this case study for two unique properties. Firstly, the high proportion of cases sequenced throughout 2020 within Victoria and the rest of Australia. Secondly, individual travel history was collected from returning travellers in Victoria during the first wave (January to May) of the coronavirus disease 2019 (COVID-19) pandemic. We found that the implementation of individual travel history was essential for the estimation of SARS-CoV-2 movement via discrete phylogeography models. Without the additional information provided by the travel history metadata, the discrete trait analysis could not be fit to the data due to numerical instability. We also suggest that during the first wave of the COVID-19 pandemic in Australia, the primary driving force behind the spread of SARS-CoV-2 was viral importation from international locations. This case study demonstrates the necessity of robust genomic datasets supplemented with epidemiological metadata for generating accurate estimates from phylogeographical models in datasets that have significant sampling bias. For future work, we recommend the collection of metadata in conjunction with genomic data. Furthermore, we highlight the risk of applying phylogeographical models to biased datasets without incorporating appropriate metadata, especially when estimates influence public health policy decision making.

Integrative omics identifies conserved and pathogen-specific responses of sepsis-causing bacteria

Even in the setting of optimal resuscitation in high-income countries severe sepsis and septic shock have a mortality of 20-40%, with antibiotic resistance dramatically increasing this mortality risk. To develop a reference dataset enabling the identification of common bacterial targets for therapeutic intervention, we applied a standardized genomic, transcriptomic, proteomic and metabolomic technological framework to multiple clinical isolates of four sepsis-causing pathogens: Escherichia coli, Klebsiella pneumoniae species complex, Staphylococcus aureus and Streptococcus pyogenes. Exposure to human serum generated a sepsis molecular signature containing global increases in fatty acid and lipid biosynthesis and metabolism, consistent with cell envelope remodelling and nutrient adaptation for osmoprotection. In addition, acquisition of cholesterol was identified across the bacterial species. This detailed reference dataset has been established as an open resource to support discovery and translational research.

Phylogenomic investigation of an outbreak of fluoroquinolone-resistant Salmonella enterica subsp. enterica serovar Paratyphi A in Phnom Penh, Cambodia

In early 2020, the Medical Biology Laboratory of the Pasteur Institute of Cambodia isolated an unusually high number of fluoroquinolone-resistant Salmonella enterica subspecies enterica serovar Paratyphi A strains during its routine bacteriological surveillance activities in Phnom Penh, Cambodia. A public-health investigation was supported by genome sequencing of these Paratyphi A strains to gain insights into the genetic diversity and population structure of a potential outbreak of fluoroquinolone-resistant paratyphoid fever. Comparative genomic and phylodynamic analyses revealed the 2020 strains were descended from a previously described 2013-2015 outbreak of Paratyphi A infections. Our analysis showed sub-lineage 2.3.1 had remained largely susceptible to fluoroquinolone drugs until 2015, but acquired chromosomal resistance to these drugs during six separate events between late 2012 and 2015. The emergence of fluoroquinolone resistance was rapidly followed by the replacement of the original susceptible Paratyphi A population, which led to a dramatic increase of fluoroquinolone-resistant blood-culture-confirmed cases in subsequent years (2016-2020). The rapid acquisition of resistance-conferring mutations in the Paratyphi A population over a 3 year period is suggestive of a strong selective pressure on that population, likely linked with fluoroquinolone use. In turn, emergence of fluoroquinolone resistance has led to increased use of extended-spectrum cephalosporins like ceftriaxone that are becoming the drug of choice for empirical treatment of paratyphoid fever in Cambodia.

An ISO-certified genomics workflow for identification and surveillance of antimicrobial resistance

Realising the promise of genomics to revolutionise identification and surveillance of antimicrobial resistance (AMR) has been a long-standing challenge in clinical and public health microbiology. Here, we report the creation and validation of abritAMR, an ISO-certified bioinformatics platform for genomics-based bacterial AMR gene detection. The abritAMR platform utilises NCBI's AMRFinderPlus, as well as additional features that classify AMR determinants into antibiotic classes and provide customised reports. We validate abritAMR by comparing with PCR or reference genomes, representing 1500 different bacteria and 415 resistance alleles. In these analyses, abritAMR displays 99.9% accuracy, 97.9% sensitivity and 100% specificity. We also compared genomic predictions of phenotype for 864 Salmonella spp. against agar dilution results, showing 98.9% accuracy. The implementation of abritAMR in our institution has resulted in streamlined bioinformatics and reporting pathways, and has been readily updated and re-verified. The abritAMR tool and validation datasets are publicly available to assist laboratories everywhere harness the power of AMR genomics in professional practice.

State-wide genomic epidemiology investigations of COVID-19 in healthcare workers in 2020 Victoria, Australia: Qualitative thematic analysis to provide insights for future pandemic preparedness

Background

COVID-19 has affected many healthcare workers (HCWs) globally. We performed state-wide SARS-CoV-2 genomic epidemiological investigations to identify HCW transmission dynamics and provide recommendations to optimise healthcare system preparedness for future outbreaks.

Methods

Genome sequencing was attempted on all COVID-19 cases in Victoria, Australia. We combined genomic and epidemiologic data to investigate the source of HCW infections across multiple healthcare facilities (HCFs) in the state. Phylogenetic analysis and fine-scale hierarchical clustering were performed for the entire dataset including community and healthcare cases. Facilities provided standardised epidemiological data and putative transmission links.

Findings

Between March-October 2020, approximately 1,240 HCW COVID-19 infection cases were identified; 765 are included here, requested for hospital investigations. Genomic sequencing was successful for 612 (80%) cases. Thirty-six investigations were undertaken across 12 HCFs. Genomic analysis revealed that multiple introductions of COVID-19 into facilities (31/36) were more common than single introductions (5/36). Major contributors to HCW acquisitions included mobility of staff and patients between wards and facilities, and characteristics and behaviours of patients that generated numerous secondary infections. Key limitations at the HCF level were identified.

Interpretation

Genomic epidemiological analyses enhanced understanding of HCW infections, revealing unsuspected clusters and transmission networks. Combined analysis of all HCWs and patients in a HCF should be conducted, supported by high rates of sequencing coverage for all cases in the population. Established systems for integrated genomic epidemiological investigations in healthcare settings will improve HCW safety in future pandemics.

Funding

The Victorian Government, the National Health and Medical Research Council Australia, and the Medical Research Future Fund.

Genomic Evidence of In-Flight SARS-CoV-2 Transmission, India to Australia, April 2021

Epidemiologic and genomic investigation of SARS-CoV-2 infections associated with 2 repatriation flights from India to Australia in April 2021 indicated that 4 passengers transmitted SARS-CoV-2 to >11 other passengers. Results suggest transmission despite mandatory mask use and predeparture testing. For subsequent flights, predeparture quarantine and expanded predeparture testing were implemented.

Niche-specific genome degradation and convergent evolution shaping Staphylococcus aureus adaptation during severe infections

During severe infections, Staphylococcus aureus moves from its colonising sites to blood and tissues and is exposed to new selective pressures, thus, potentially driving adaptive evolution. Previous studies have shown the key role of the agr locus in S. aureus pathoadaptation; however, a more comprehensive characterisation of genetic signatures of bacterial adaptation may enable prediction of clinical outcomes and reveal new targets for treatment and prevention of these infections. Here, we measured adaptation using within-host evolution analysis of 2590 S. aureus genomes from 396 independent episodes of infection. By capturing a comprehensive repertoire of single nucleotide and structural genome variations, we found evidence of a distinctive evolutionary pattern within the infecting populations compared to colonising bacteria. These invasive strains had up to 20-fold enrichments for genome degradation signatures and displayed significantly convergent mutations in a distinctive set of genes, linked to antibiotic response and pathogenesis. In addition to agr-mediated adaptation, we identified non-canonical, genome-wide significant loci including sucA-sucB and stp1. The prevalence of adaptive changes increased with infection extent, emphasising the clinical significance of these signatures. These findings provide a high-resolution picture of the molecular changes when S. aureus transitions from colonisation to severe infection and may inform correlation of infection outcomes with adaptation signatures.

The bacterium Staphylococcus aureus lives harmlessly on our skin and noses. However, occasionally, it gets into our blood and internal organs, such as our bones and joints, where it causes severe, long-lasting infections that are difficult to treat.

Over time, S. aureus acquire characteristics that help them to adapt to different locations, such as transitioning from the nose to the blood, and avoid being killed by antibiotics. Previous studies have identified changes, or ‘mutations’, in genes that are likely to play an important role in this evolutionary process. One of these genes, called accessory gene regulator (or agr for short), has been shown to control the mechanisms S. aureus use to infect cells and disseminate in the body. However, it is unclear if there are changes in other genes that also help S. aureus adapt to life inside the human body.

To help resolve this mystery, Giulieri et al. collected 2,500 samples of S. aureus from almost 400 people. This included bacteria harmlessly living on the skin or in the nose, as well as strains that caused an infection. Gene sequencing revealed a small number of genes, referred to as ‘adaptive genes’, that often acquire mutations during infection. Of these, agr was the most commonly altered. However, mutations in less well-known genes were also identified: some of these genes are related to resistance to antibiotics, while others are involved in chemical processes that help the bacteria to process nutrients.

Most mutations were caused by random errors being introduced in to the bacteria’s genetic code which stopped genes from working. However, in some cases, genes were turned off by small fragments of DNA moving around and inserting themselves into different parts of the genome.

This study highlights a group of genes that help S. aureus to thrive inside the body and cause severe and prolonged infections. If these results can be confirmed, it may help to guide which antibiotics are used to treat different infections. Furthermore, understanding which genes are important for infection could lead to new strategies for eliminating this dangerous bacterium.

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).

SARS-CoV-2 infection results in immune responses in the respiratory tract and peripheral blood that suggest mechanisms of disease severity

Respiratory tract infection with SARS-CoV-2 results in varying immunopathology underlying COVID-19. We examine cellular, humoral and cytokine responses covering 382 immune components in longitudinal blood and respiratory samples from hospitalized COVID-19 patients. SARS-CoV-2-specific IgM, IgG, IgA are detected in respiratory tract and blood, however, receptor-binding domain (RBD)-specific IgM and IgG seroconversion is enhanced in respiratory specimens. SARS-CoV-2 neutralization activity in respiratory samples correlates with RBD-specific IgM and IgG levels. Cytokines/chemokines vary between respiratory samples and plasma, indicating that inflammation should be assessed in respiratory specimens to understand immunopathology. IFN-α2 and IL-12p70 in endotracheal aspirate and neutralization in sputum negatively correlate with duration of hospital stay. Diverse immune subsets are detected in respiratory samples, dominated by neutrophils. Importantly, dexamethasone treatment does not affect humoral responses in blood of COVID-19 patients. Our study unveils differential immune responses between respiratory samples and blood, and shows how drug therapy affects immune responses during COVID-19.

Genomic Characterisation Reveals a Dominant Lineage of SARS-CoV-2 in Papua New Guinea

The coronavirus disease pandemic has highlighted the utility of pathogen genomics as a key part of comprehensive public health response to emerging infectious diseases threats, however, the ability to generate, analyse, and respond to pathogen genomic data varies around the world. Papua New Guinea (PNG), which has limited in-country capacity for genomics, has experienced significant outbreaks of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with initial genomics data indicating a large proportion of cases were from lineages that are not well defined within the current nomenclature. Through a partnership between in-country public health agencies and academic organisations, industry, and a public health genomics reference laboratory in Australia a system for routine SARS-CoV-2 genomics from PNG was established. Here we aim to characterise and describe the genomics of PNG’s second wave and examine the sudden expansion of a lineage that is not well defined but very prevalent in the Western Pacific region. We generated 1797 sequences from cases in PNG and performed phylogenetic and phylodynamic analyses to examine the outbreak and characterise the circulating lineages and clusters present. Our results reveal the rapid expansion of the B.1.466.2 and related lineages within PNG, from multiple introductions into the country. We also highlight the difficulties that unstable lineage assignment causes when using genomics to assist with rapid cluster definitions.

AusTrakka: Fast-tracking nationalized genomics surveillance in response to the COVID-19 pandemic

The COVID-19 pandemic has driven demand for integrated genomics, resulting in fast-tracked development of AusTrakka, Australia’s pathogen genomics platform. This facilitated rapid data sharing, democratised access to computational and bioinformatic resources and expertise, and achieved national real-time genomic surveillance.

The COVID-19 pandemic has accelerated the demand for near real-time analysis and dissemination of pathogen genomic data. In this comment, the authors describe how Australia has developed and rolled out its SARS-CoV-2 genomics platform, AusTrakka, and used it to support public health action.

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.

Genomic diversity of antimicrobial resistance in non-typhoidal Salmonella in Victoria, Australia

Non-typhoidal Salmonella (NTS) is the second most common cause of foodborne bacterial gastroenteritis in Australia with antimicrobial resistance (AMR) increasing in recent years. Whole-genome sequencing (WGS) provides opportunities for in silico detection of AMR determinants. The objectives of this study were two-fold: (1) establish the utility of WGS analyses for inferring phenotypic resistance in NTS, and (2) explore clinically relevant genotypic AMR profiles to third generation cephalosporins (3GC) in NTS lineages. The concordance of 2490 NTS isolates with matched WGS and phenotypic susceptibility data against 13 clinically relevant antimicrobials was explored. In silico serovar prediction and typing was performed on assembled reads and interrogated for known AMR determinants. The surrounding genomic context, plasmid determinants and co-occurring AMR patterns were further investigated for multidrug resistant serovars harbouring bla_CMY-2, bla_CTX-M-55 or bla_CTX-M-65. Our data demonstrated a high correlation between WGS and phenotypic susceptibility testing. Phenotypic-genotypic concordance was observed between 2440/2490 (98.0 %) isolates, with overall sensitivity and specificity rates >98 % and positive and negative predictive values >97 %. The most common AMR determinants were bla_TEM-1, sul2, tet(A), strA-strB and floR. Phenotypic resistance to cefotaxime and azithromycin was low and observed in 6.2 % (151/2486) and 0.9 % (16/1834) of the isolates, respectively. Several multi-drug resistant NTS lineages were resistant to 3GC due to different genetic mechanisms including bla_CMY-2, bla_CTX-M-55 or bla_CTX-M-65. This study shows WGS can enhance existing AMR surveillance in NTS datasets routinely produced in public health laboratories to identify emerging AMR in NTS. These approaches will be critical for developing capacity to detect emerging public health threats such as resistance to 3GC.

Search and Contain: Impact of an Integrated Genomic and Epidemiological Surveillance and Response Program for Control of Carbapenemase-producing Enterobacterales

BACKGROUND

Multiresistant organisms (MROs) pose a critical threat to public health. Population-based programs for control of MROs such as carbapenemase-producing Enterobacterales (CPE) have emerged and evaluation is needed. We assessed the feasibility and impact of a statewide CPE surveillance and response program deployed across Victoria, Australia (population 6.5 million).

METHODS

A prospective multimodal intervention including active screening, carrier isolation, centralized case investigation, and comparative pathogen genomics was implemented. We analyzed trends in CPE incidence and clinical presentation, risk factors, and local transmission over the program's first 3 years (2016-2018).

RESULTS

CPE case ascertainment increased over the study period to 1.42 cases/100 000 population, linked to increased screening without a concomitant rise in active clinical infections (0.45-0.60 infections/100 000 population, P = .640). KPC-2 infection decreased from 0.29 infections/100 000 population prior to intervention to 0.03 infections/100 000 population in 2018 (P = .003). Comprehensive case investigation identified instances of overseas community acquisition. Median time between isolate referral and genomic and epidemiological assessment for local transmission was 11 days (IQR, 9-14). Prospective surveillance identified numerous small transmission networks (median, 2; range, 1-19 cases), predominantly IMP and KPC, with median pairwise distance of 8 (IQR, 4-13) single nucleotide polymorphisms; low diversity between clusters of the same sequence type suggested genomic cluster definitions alone are insufficient for targeted response.

CONCLUSIONS

We demonstrate the value of centralized CPE control programs to increase case ascertainment, resolve risk factors, and identify local transmission through prospective genomic and epidemiological surveillance; methodologies are transferable to low-prevalence settings and MROs globally.

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.

Tuberculosis in Australia's tropical north: a population-based genomic epidemiological study

BACKGROUND

The Northern Territory (NT) has the highest tuberculosis (TB) rate of all Australian jurisdictions. We combined TB public health surveillance data with genomic sequencing of Mycobacterium tuberculosis isolates in the tropical 'Top End' of the NT to investigate trends in TB incidence and transmission.

METHODS

This retrospective observational study included all 741 culture-confirmed cases of TB in the Top End over three decades from 1989-2020. All 497 available M. tuberculosis isolates were sequenced. We used contact tracing data to define a threshold pairwise SNP distance for hierarchical single linkage clustering, and examined putative transmission clusters in the context of epidemiologic information.

FINDINGS

There were 359 (48%) cases born overseas, 329 (44%) cases among Australian First Nations peoples, and 52 (7%) cases were Australian-born and non-Indigenous. The annual incidence in First Nations peoples from 1989-2019 fell from average 50.4 to 11.0 per 100,000 (P<0·001). First Nations cases were more likely to die from TB (41/329, 12·5%) than overseas-born cases (11/359, 3·1%; P<0·001). Using a threshold of ≤12 SNPs, 28 clusters of between 2-64 individuals were identified, totalling 250 cases; 214 (86%) were First Nations cases and 189 (76%) were from a remote region. The time between cases and past epidemiologically- and genomically-linked contacts ranged from 4·5 months to 24 years.

INTERPRETATION

Our findings support prioritisation of timely case detection, contact tracing augmented by genomic sequencing, and latent TB treatment to break transmission chains in Top End remote hotspot regions.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2 with grinch

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

Accessible Platform for High-Throughput COVID-19 Molecular Diagnostics and Genome Sequencing Using a Repurposed 3D Printer for RNA Extraction

The COVID-19 pandemic has exposed the dependence of diagnostic laboratories on a handful of large corporations with market monopolies on the worldwide supply of reagents, consumables, and hardware for molecular diagnostics. Global shortages of key consumables for RT-qPCR detection of SARS-CoV-2 RNA have impaired the ability to run essential, routine diagnostic services. Here, we describe a workflow for rapid detection of SARS-CoV-2 RNA in upper respiratory samples including nasal swabs and saliva, utilizing low-cost equipment and readily accessible reagents. Using repurposed Creality3D Ender-3 three-dimensional (3D) printers, we built a semiautomated paramagnetic bead RNA extraction platform. The hardware for the system was built for $300 USD, and the material cost per reaction was $1 USD. Named the Ender VX500, instrument performance when paired with RT-qPCR for SARS-CoV-2 detection in nasal and saliva specimens was two virus copies per microliter. There was a high-performance agreement (assessed using 458 COVID-19 nasal swab specimens) with the Aptima SARS-CoV-2 assay run on the Hologic Panther, a commercial automated RNA extraction and detection platform. Inter- and intrainstrument precision was excellent (coefficients of variation (CoV) of 1.10 and 0.66-1.32%, respectively) across four instruments. The platform is scalable with throughput ranging from 23 specimens on a single instrument run by one user in 50 min to 364 specimens on four instruments run by four users in 190 min. Step-by-step instructions and protocols for building and running the Ender VX500 have been made available without restriction.

Defective Severe Acute Respiratory Syndrome Coronavirus 2 Immune Responses in an Immunocompromised Individual With Prolonged Viral Replication

We describe severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-specific immune responses in a patient with lymphoma and recent programmed death 1 (PD-1) inhibitor therapy with late onset of severe coronavirus disease 2019 disease and prolonged SARS-CoV-2 replication, in comparison to age-matched and immunocompromised controls. High levels of HLA-DR+/CD38+ activation, interleukin 6, and interleukin 18 in the absence of B cells and PD-1 expression was observed. SARS-CoV-2-specific antibody responses were absent and SARS-CoV-2-specific T cells were minimally detected. This case highlights challenges in managing immunocompromised hosts who may fail to mount effective virus-specific immune responses.

Immune responses in COVID-19 respiratory tract and blood reveal mechanisms of disease severity

Although the respiratory tract is the primary site of SARS-CoV-2 infection and the ensuing immunopathology, respiratory immune responses are understudied and urgently needed to understand mechanisms underlying COVID-19 disease pathogenesis. We collected paired longitudinal blood and respiratory tract samples (endotracheal aspirate, sputum or pleural fluid) from hospitalized COVID-19 patients and non-COVID-19 controls. Cellular, humoral and cytokine responses were analysed and correlated with clinical data. SARS-CoV-2-specific IgM, IgG and IgA antibodies were detected using ELISA and multiplex assay in both the respiratory tract and blood of COVID-19 patients, although a higher receptor binding domain (RBD)-specific IgM and IgG seroconversion level was found in respiratory specimens. SARS-CoV-2 neutralization activity in respiratory samples was detected only when high levels of RBD-specific antibodies were present. Strikingly, cytokine/chemokine levels and profiles greatly differed between respiratory samples and plasma, indicating that inflammation needs to be assessed in respiratory specimens for the accurate assessment of SARS-CoV-2 immunopathology. Diverse immune cell subsets were detected in respiratory samples, albeit dominated by neutrophils. Importantly, we also showed that dexamethasone and/or remdesivir treatment did not affect humoral responses in blood of COVID-19 patients. Overall, our study unveils stark differences in innate and adaptive immune responses between respiratory samples and blood and provides important insights into effect of drug therapy on immune responses in COVID-19 patients.

Genomics-informed responses in the elimination of COVID-19 in Victoria, Australia: an observational, genomic epidemiological study

BACKGROUND

A cornerstone of Australia's ability to control COVID-19 has been effective border control with an extensive supervised quarantine programme. However, a rapid recrudescence of COVID-19 was observed in the state of Victoria in June, 2020. We aim to describe the genomic findings that located the source of this second wave and show the role of genomic epidemiology in the successful elimination of COVID-19 for a second time in Australia.

METHODS

In this observational, genomic epidemiological study, we did genomic sequencing of all laboratory-confirmed cases of COVID-19 diagnosed in Victoria, Australia between Jan 25, 2020, and Jan 31, 2021. We did phylogenetic analyses, genomic cluster discovery, and integrated results with epidemiological data (detailed information on demographics, risk factors, and exposure) collected via interview by the Victorian Government Department of Health. Genomic transmission networks were used to group multiple genomic clusters when epidemiological and genomic data suggested they arose from a single importation event and diversified within Victoria. To identify transmission of emergent lineages between Victoria and other states or territories in Australia, all publicly available SARS-CoV-2 sequences uploaded before Feb 11, 2021, were obtained from the national sequence sharing programme AusTrakka, and epidemiological data were obtained from the submitting laboratories. We did phylodynamic analyses to estimate the growth rate, doubling time, and number of days from the first local infection to the collection of the first sequenced genome for the dominant local cluster, and compared our growth estimates to previously published estimates from a similar growth phase of lineage B.1.1.7 (also known as the Alpha variant) in the UK.

FINDINGS

Between Jan 25, 2020, and Jan 31, 2021, there were 20 451 laboratory-confirmed cases of COVID-19 in Victoria, Australia, of which 15 431 were submitted for sequencing, and 11 711 met all quality control metrics and were included in our analysis. We identified 595 genomic clusters, with a median of five cases per cluster (IQR 2-11). Overall, samples from 11 503 (98·2%) of 11 711 cases clustered with another sample in Victoria, either within a genomic cluster or transmission network. Genomic analysis revealed that 10 426 cases, including 10 416 (98·4%) of 10 584 locally acquired cases, diagnosed during the second wave (between June and October, 2020) were derived from a single incursion from hotel quarantine, with the outbreak lineage (transmission network G, lineage D.2) rapidly detected in other Australian states and territories. Phylodynamic analyses indicated that the epidemic growth rate of the outbreak lineage in Victoria during the initial growth phase (samples collected between June 4 and July 9, 2020; 47·4 putative transmission events, per branch, per year [1/years; 95% credible interval 26·0-85·0]), was similar to that of other reported variants, such as B.1.1.7 in the UK (mean approximately 71·5 1/years). Strict interventions were implemented, and the outbreak lineage has not been detected in Australia since Oct 29, 2020. Subsequent cases represented independent international or interstate introductions, with limited local spread.

INTERPRETATION

Our study highlights how rapid escalation of clonal outbreaks can occur from a single incursion. However, strict quarantine measures and decisive public health responses to emergent cases are effective, even with high epidemic growth rates. Real-time genomic surveillance can alter the way in which public health agencies view and respond to COVID-19 outbreaks.

FUNDING

The Victorian Government, the National Health and Medical Research Council Australia, and the Medical Research Future Fund.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

Complete microbial genomes for public health in Australia and the Southwest Pacific

Complete genomes of microbial pathogens are essential for the phylogenomic analyses that increasingly underpin core public health laboratory activities. Here, we announce a BioProject (PRJNA556438) dedicated to sharing complete genomes chosen to represent a range of pathogenic bacteria with regional importance to Australia and the Southwest Pacific; enriching the catalogue of globally available complete genomes for public health while providing valuable strains to regional public health microbiology laboratories. In this first step, we present 26 complete high-quality bacterial genomes. Additionally, we describe here a framework for reconstructing complete microbial genomes and highlight some of the challenges and considerations for accurate and reproducible genome reconstruction.

Case Report: Confirmation by Metagenomic Sequencing of Visceral Leishmaniasis in an Immunosuppressed Returned Traveler

There has been increased interest in using metagenomic next-generation sequencing as an unbiased approach for diagnosing infectious diseases. We describe a 61-year-old man on fingolimod therapy for multiple sclerosis with an extensive travel history who presented with 7 months of fevers, night sweats, and weight loss. Peripheral blood tests showed pancytopenia and abnormal acute phase reactants. A bone marrow aspirate showed the presence of numerous intracellular and extracellular amastigotes consistent with visceral leishmaniasis (VL). Metagenomic sequencing of the bone marrow aspirate confirmed Leishmania infantum, a species widely reported in the Mediterranean region. This correlated with acquisition of VL infection during the patient's most recent epidemiological exposure in southern Italy 12 months prior. This case demonstrates the potential application of metagenomic sequencing for identification and speciation of Leishmania in cases of VL; however, further assessment is required using other more readily obtained clinical samples such as blood.

Viral genomics to inform infection control response in occupational COVID-19 transmission

Healthcare workers are at increased risk of occupational transmission of SARS-CoV-2. We report two instances of healthcare workers contracting SARS-CoV-2 despite no known breach of personal protective equipment. Additional specific equipment cleaning was initiated. Viral genomic sequencing supported this transmission hypothesis and our subsequent response.

Co-circulation of Multidrug-resistant Shigella Among Men Who Have Sex With Men in Australia

BACKGROUND

In urban Australia, the burden of shigellosis is either in returning travelers from shigellosis-endemic regions or in men who have sex with men (MSM). Here, we combine genomic data with comprehensive epidemiological data on sexual exposure and travel to describe the spread of multidrug-resistant Shigella lineages.

METHODS

A population-level study of all cultured Shigella isolates in the state of Victoria, Australia, was undertaken from 1 January 2016 through 31 March 2018. Antimicrobial susceptibility testing, whole-genome sequencing, and bioinformatic analyses of 545 Shigella isolates were performed at the Microbiological Diagnostic Unit Public Health Laboratory. Risk factor data on travel and sexual exposure were collected through enhanced surveillance forms or by interviews.

RESULTS

Rates of antimicrobial resistance were high, with 17.6% (95/541) and 50.6% (274/541) resistance to ciprofloxacin and azithromycin, respectively. There were strong associations between antimicrobial resistance, phylogeny, and epidemiology. Specifically, 2 major MSM-associated lineages were identified: a Shigellasonnei lineage (n = 159) and a Shigella flexneri 2a lineage (n = 105). Of concern, 147/159 (92.4%) of isolates within the S. sonnei MSM-associated lineage harbored mutations associated with reduced susceptibility to recommended oral antimicrobials: namely, azithromycin, trimethoprim-sulfamethoxazole, and ciprofloxacin. Long-read sequencing demonstrated global dissemination of multidrug-resistant plasmids across Shigella species and lineages, but predominantly associated with MSM isolates.

CONCLUSIONS

Our contemporary data highlight the ongoing public health threat posed by resistant Shigella, both in Australia and globally. Urgent multidisciplinary public health measures are required to interrupt transmission and prevent infection.

Tracking the COVID-19 pandemic in Australia using genomics

Genomic sequencing has significant potential to inform public health management for SARS-CoV-2. Here we report high-throughput genomics for SARS-CoV-2, sequencing 80% of cases in Victoria, Australia (population 6.24 million) between 6 January and 14 April 2020 (total 1,333 COVID-19 cases). We integrate epidemiological, genomic and phylodynamic data to identify clusters and impact of interventions. The global diversity of SARS-CoV-2 is represented, consistent with multiple importations. Seventy-six distinct genomic clusters were identified, including large clusters associated with social venues, healthcare and cruise ships. Sequencing sequential samples from 98 patients reveals minimal intra-patient SARS-CoV-2 genomic diversity. Phylodynamic modelling indicates a significant reduction in the effective viral reproductive number (R_e) from 1.63 to 0.48 after implementing travel restrictions and physical distancing. Our data provide a concrete framework for the use of SARS-CoV-2 genomics in public health responses, including its use to rapidly identify SARS-CoV-2 transmission chains, increasingly important as social restrictions ease globally.

Validation of a single-step, single-tube reverse transcription loop-mediated isothermal amplification assay for rapid detection of SARS-CoV-2 RNA

Introduction. The SARS-CoV-2 pandemic of 2020 has resulted in unparalleled requirements for RNA extraction kits and enzymes required for virus detection, leading to global shortages. This has necessitated the exploration of alternative diagnostic options to alleviate supply chain issues.Aim. To establish and validate a reverse transcription loop-mediated isothermal amplification (RT- LAMP) assay for the detection of SARS-CoV-2 from nasopharyngeal swabs.Methodology. We used a commercial RT-LAMP mastermix from OptiGene in combination with a primer set designed to detect the CDC N1 region of the SARS-CoV-2 nucleocapsid (N) gene. A single-tube, single-step fluorescence assay was implemented whereby 1 µl of universal transport medium (UTM) directly from a nasopharyngeal swab could be used as template, bypassing the requirement for RNA purification. Amplification and detection could be conducted in any thermocycler capable of holding 65 °C for 30 min and measure fluorescence in the FAM channel at 1 min intervals.Results. Assay evaluation by assessment of 157 clinical specimens previously screened by E-gene RT-qPCR revealed assay sensitivity and specificity of 87 and 100%, respectively. Results were fast, with an average time-to-positive (Tp) for 93 clinical samples of 14 min (sd±7 min). Using dilutions of SARS-CoV-2 virus spiked into UTM, we also evaluated assay performance against FDA guidelines for implementation of emergency-use diagnostics and established a limit-of-detection of 54 Tissue Culture Infectious Dose 50 per ml (TCID50 ml-1), with satisfactory assay sensitivity and specificity. A comparison of 20 clinical specimens between four laboratories showed excellent interlaboratory concordance; performing equally well on three different, commonly used thermocyclers, pointing to the robustness of the assay.Conclusion. With a simplified workflow, The N1 gene Single Tube Optigene LAMP assay (N1-STOP-LAMP) is a powerful, scalable option for specific and rapid detection of SARS-CoV-2 and an additional resource in the diagnostic armamentarium against COVID-19.

Incursions of Candida auris into Australia, 2018

Candida auris is an emerging global healthcare-associated pathogen. During July–December 2018, four patients with C. auris were identified in Victoria, Australia, all with previous overseas hospitalization. Phylogenetic analysis revealed putative transmission between 2 patients and suspected overseas acquisition in the others. Vigilant screening of at-risk patients is required.

Isolation and rapid sharing of the 2019 novel coronavirus (SARS-CoV-2) from the first patient diagnosed with COVID-19 in Australia

Objectives

To describe the first isolation and sequencing of SARS‐CoV‐2 in Australia and rapid sharing of the isolate.

Setting

SARS‐CoV‐2 was isolated from a 58‐year‐old man from Wuhan, China who arrived in Melbourne on 19 January 2020 and was admitted to the Monash Medical Centre, Melbourne from the emergency department on 24 January 2020 with fever, cough, and progressive dyspnoea.

Major outcomes

Clinical course and laboratory features of the first reported case of COVID‐19 (the illness caused by SARS‐CoV‐2) in Australia; isolation, whole genome sequencing, imaging, and rapid sharing of virus from the patient.

Results

A nasopharyngeal swab and sputum collected when the patient presented to hospital were each positive for SARS‐CoV‐2 (reverse transcription polymerase chain reaction). Inoculation of Vero/hSLAM cells with material from the nasopharyngeal swab led to the isolation of SARS‐CoV‐2 virus in culture. Electron microscopy of the supernatant confirmed the presence of virus particles with morphology characteristic of viruses of the family Coronaviridae. Whole genome sequencing of the viral isolate and phylogenetic analysis indicated the isolate exhibited greater than 99.99% sequence identity with other publicly available SARS‐CoV‐2 genomes. Within 24 hours of isolation, the first Australian SARS‐CoV‐2 isolate was shared with local and overseas reference laboratories and major North American and European culture collections.

Conclusions

The ability to rapidly identify, propagate, and internationally share our SARS‐CoV‐2 isolate is an important step in collaborative scientific efforts to deal effectively with this international public health emergency by developing better diagnostic procedures, vaccine candidates, and antiviral agents.

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.

The changing landscape of vancomycin-resistant Enterococcus faecium in Australia: a population-level genomic study

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.

Increasing tolerance of hospital Enterococcus faecium to handwash alcohols

Alcohol-based disinfectants and particularly hand rubs are a key way to control hospital infections worldwide. Such disinfectants restrict transmission of pathogens, such as multidrug-resistant Staphylococcus aureus and Enterococcus faecium Despite this success, health care infections caused by E. faecium are increasing. We tested alcohol tolerance of 139 hospital isolates of E. faecium obtained between 1997 and 2015 and found that E. faecium isolates after 2010 were 10-fold more tolerant to killing by alcohol than were older isolates. Using a mouse gut colonization model of E. faecium transmission, we showed that alcohol-tolerant E. faecium resisted standard 70% isopropanol surface disinfection, resulting in greater mouse gut colonization compared to alcohol-sensitive E. faecium We next looked for bacterial genomic signatures of adaptation. Alcohol-tolerant E. faecium accumulated mutations in genes involved in carbohydrate uptake and metabolism. Mutagenesis confirmed the roles of these genes in the tolerance of E. faecium to isopropanol. These findings suggest that bacterial adaptation is complicating infection control recommendations, necessitating additional procedures to prevent E. faecium from spreading in hospital settings.

Bionitio: demonstrating and facilitating best practices for bioinformatics command-line software

BACKGROUND

Bioinformatics software tools are often created ad hoc, frequently by people without extensive training in software development. In particular, for beginners, the barrier to entry in bioinformatics software development is high, especially if they want to adopt good programming practices. Even experienced developers do not always follow best practices. This results in the proliferation of poorer-quality bioinformatics software, leading to limited scalability and inefficient use of resources; lack of reproducibility, usability, adaptability, and interoperability; and erroneous or inaccurate results.

FINDINGS

We have developed Bionitio, a tool that automates the process of starting new bioinformatics software projects following recommended best practices. With a single command, the user can create a new well-structured project in 1 of 12 programming languages. The resulting software is functional, carrying out a prototypical bioinformatics task, and thus serves as both a working example and a template for building new tools. Key features include command-line argument parsing, error handling, progress logging, defined exit status values, a test suite, a version number, standardized building and packaging, user documentation, code documentation, a standard open source software license, software revision control, and containerization.

CONCLUSIONS

Bionitio serves as a learning aid for beginner-to-intermediate bioinformatics programmers and provides an excellent starting point for new projects. This helps developers adopt good programming practices from the beginning of a project and encourages high-quality tools to be developed more rapidly. This also benefits users because tools are more easily installed and consistent in their usage. Bionitio is released as open source software under the MIT License and is available at https://github.com/bionitio-team/bionitio.

Zinc-binding to the cytoplasmic PAS domain regulates the essential WalK histidine kinase of Staphylococcus aureus

WalKR (YycFG) is the only essential two-component regulator in the human pathogen Staphylococcus aureus. WalKR regulates peptidoglycan synthesis, but this function alone does not explain its essentiality. Here, to further understand WalKR function, we investigate a suppressor mutant that arose when WalKR activity was impaired; a histidine to tyrosine substitution (H271Y) in the cytoplasmic Per-Arnt-Sim (PAS^CYT) domain of the histidine kinase WalK. Introducing the WalK^H271Y mutation into wild-type S. aureus activates the WalKR regulon. Structural analyses of the WalK PAS^CYT domain reveal a metal-binding site, in which a zinc ion (Zn²⁺) is tetrahedrally-coordinated by four amino acids including H271. The WalK^H271Y mutation abrogates metal binding, increasing WalK kinase activity and WalR phosphorylation. Thus, Zn²⁺-binding negatively regulates WalKR. Promoter-reporter experiments using S. aureus confirm Zn²⁺ sensing by this system. Identification of a metal ligand recognized by the WalKR system broadens our understanding of this critical S. aureus regulon.

Emergence and divergence of major lineages of Shiga-toxin-producing Escherichia coli in Australia

Shiga-toxin-producing Escherichia coli (STEC) infection is an important global cause of foodborne disease. To date however, genomics-based studies of STEC have been predominately focused upon STEC collected in the Northern Hemisphere. Here, we demonstrate the population structure of 485 STEC isolates in Australia, and show that several clonal groups (CGs) common to Australia were infrequently detected in a representative selection of contemporary STEC genomes from around the globe. Further, phylogenetic analysis demonstrated that lineage II of the global O157:H7 STEC was most prevalent in Australia, and was characterized by a frameshift mutation in flgF, resulting in the H-non-motile phenotype. Strong concordance between in silico and phenotypic serotyping was observed, along with concordance between in silico and conventional detection of stx genes. These data represent the most comprehensive STEC analysis from the Southern Hemisphere, and provide a framework for future national genomics-based surveillance of STEC in Australia.

Remodeling of pSK1 Family Plasmids and Enhanced Chlorhexidine Tolerance in a Dominant Hospital Lineage of Methicillin-Resistant Staphylococcus aureus

Staphylococcus aureus is a significant human pathogen whose evolution and adaptation have been shaped in part by mobile genetic elements (MGEs), facilitating the global spread of extensive antimicrobial resistance. However, our understanding of the evolutionary dynamics surrounding MGEs, in particular, how changes in the structure of multidrug resistance (MDR) plasmids may influence important staphylococcal phenotypes, is incomplete. Here, we undertook a population and functional genomics study of 212 methicillin-resistant S. aureus (MRSA) sequence type 239 (ST239) isolates collected over 32 years to explore the evolution of the pSK1 family of MDR plasmids, illustrating how these plasmids have coevolved with and contributed to the successful adaptation of this persistent MRSA lineage. Using complete genomes and temporal phylogenomics, we reconstructed the evolution of the pSK1 family lineage from its emergence in the late 1970s and found that multiple structural variants have arisen. Plasmid maintenance and stability were linked to IS256- and IS257-mediated chromosomal integration and disruption of the plasmid replication machinery. Overlaying genomic comparisons with phenotypic susceptibility data for gentamicin, trimethoprim, and chlorhexidine, it appeared that pSK1 has contributed to enhanced resistance in ST239 MRSA isolates through two mechanisms: (i) acquisition of plasmid-borne resistance mechanisms increasing the rates of gentamicin resistance and reduced chlorhexidine susceptibility and (ii) changes in the plasmid configuration linked with further enhancement of chlorhexidine tolerance. While the exact mechanism of enhanced tolerance remains elusive, this research has uncovered a potential evolutionary response of ST239 MRSA to biocides, one of which may contribute to the ongoing persistence and adaptation of this lineage within health care institutions.

Wave 2 strains of atypical Vibrio cholerae El Tor caused the 2009-2011 cholera outbreak in Papua New Guinea

Vibrio cholerae is the causative agent of cholera, a globally important human disease for at least 200 years. In 2009–2011, the first recorded cholera outbreak in Papua New Guinea (PNG) occurred. We conducted genetic and phenotypic characterization of 21 isolates of V. cholerae, with whole-genome sequencing conducted on 2 representative isolates. The PNG outbreak was caused by an atypical El Tor strain harbouring a tandem repeat of the CTX prophage on chromosome II. Whole-genome sequence data, prophage structural analysis and the absence of the SXT integrative conjugative element was indicative that the PNG isolates were most closely related to strains previously isolated in South-East and East Asia with affiliations to global wave 2 strains. This finding suggests that the cholera outbreak in PNG was caused by an exotic (non-endemic) strain of V. cholerae that originated in South-East Asia.

Genomic epidemiology and antimicrobial resistance of Neisseria gonorrhoeae in New Zealand

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.

Global spread of three multidrug-resistant lineages of Staphylococcus epidermidis

Staphylococcus epidermidis is a conspicuous member of the human microbiome, widely present on healthy skin. Here we show that S. epidermidis has also evolved to become a formidable nosocomial pathogen. Using genomics, we reveal that three multidrug-resistant, hospital-adapted lineages of S. epidermidis (two ST2 and one ST23) have emerged in recent decades and spread globally. These lineages are resistant to rifampicin through acquisition of specific rpoB mutations that have become fixed in the populations. Analysis of isolates from 96 institutions in 24 countries identified dual D471E and I527M RpoB substitutions to be the most common cause of rifampicin resistance in S. epidermidis, accounting for 86.6% of mutations. Furthermore, we reveal that the D471E and I527M combination occurs almost exclusively in isolates from the ST2 and ST23 lineages. By breaching lineage-specific DNA methylation restriction modification barriers and then performing site-specific mutagenesis, we show that these rpoB mutations not only confer rifampicin resistance, but also reduce susceptibility to the last-line glycopeptide antibiotics, vancomycin and teicoplanin. Our study has uncovered the previously unrecognized international spread of a near pan-drug-resistant opportunistic pathogen, identifiable by a rifampicin-resistant phenotype. It is possible that hospital practices, such as antibiotic monotherapy utilizing rifampicin-impregnated medical devices, have driven the evolution of this organism, once trivialized as a contaminant, towards potentially incurable infections.

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 ).

Genomic exploration of sequential clinical isolates reveals a distinctive molecular signature of persistent Staphylococcus aureus bacteraemia

BACKGROUND

Large-scale genomic studies of within-host diversity in Staphylococcus aureus bacteraemia (SAB) are needed to understanding bacterial adaptation underlying persistence and thus refining the role of genomics in management of SAB. However, available comparative genomic studies of sequential SAB isolates have tended to focus on selected cases of unusually prolonged bacteraemia, where secondary antimicrobial resistance has developed.

METHODS

To understand bacterial genetic diversity during SAB more broadly, we applied whole genome sequencing to a large collection of sequential isolates obtained from patients with persistent or relapsing bacteraemia. After excluding genetically unrelated isolates, we performed an in-depth genomic analysis of point mutations and chromosome structural variants arising within individual SAB episodes.

RESULTS

We show that, while adaptation pathways are heterogenous and episode-specific, isolates from persistent bacteraemia have a distinctive molecular signature, characterised by a low mutation frequency and high proportion of non-silent mutations. Analysis of structural genomic variants revealed that these often overlooked genetic events are commonly acquired during SAB. We discovered that IS256 insertion may represent the most effective driver of within-host microevolution in selected lineages, with up to three new insertion events per isolate even in the absence of other mutations. Genetic mechanisms resulting in significant phenotypic changes, such as increases in vancomycin resistance, development of small colony phenotypes, and decreases in cytotoxicity, included mutations in key genes (rpoB, stp, agrA) and an IS256 insertion upstream of the walKR operon.

CONCLUSIONS

This study provides for the first time a large-scale analysis of within-host genomic changes during invasive S. aureus infection and describes specific patterns of adaptation that will be informative for both understanding S. aureus pathoadaptation and utilising genomics for management of complicated S. aureus infections.