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.

Epidemiology of carbapenem resistant Acinetobacter baumannii in New Zealand.

AIM: Carbapenem resistant Acinetobacter baumannii have limited treatment options and a propensity to cause hospital outbreaks. In recent years an increase in their detection has been observed in New Zealand. This study aimed to describe the molecular epidemiology of these isolates. METHOD: This study utilised carbapenem resistant A. baumannii complex isolates identified across New Zealand between January 2010 to April 2018. Whole genome sequence analysis and associated demographic information was used to contextualise local isolates within the global epidemiology and establish the relationship between isolates. RESULTS: Thirty-three carbapenem resistant A. baumannii complex isolates (31 A. baumannii sensu stricto) were identified. Twenty-four (73%) were from January 2015 onwards. Twenty-four (73%) had an identifiable epidemiological link to overseas hospitalisation. Twenty-three (74%) of 31 A. baumannii sensu stricto were sequence type (ST) 2 (Pasteur scheme). Phylogenetic analysis identified three ST2 clusters. The largest cluster, of 12 isolates, was from 2015 onwards; with nine (75%) associated with recent hospitalisation in Fiji or Samoa. CONCLUSION: Increasing numbers of carbapenem resistant A. baumannii are being identified in New Zealand. Our data show that this is in large part associated with transnational spread of a single A. baumannii sensu stricto ST 2 strain between Fiji, Samoa and New Zealand.

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

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.

Reconstruction of the Genomes of Drug-Resistant Pathogens for Outbreak Investigation through Metagenomic Sequencing.

Culture-independent methods that target genome fragments have shown promise in identifying certain pathogens, but the holy grail of comprehensive pathogen genome detection from microbiologically complex samples for subsequent forensic analyses remains a challenge. In the context of an investigation of a nosocomial outbreak, we used shotgun metagenomic sequencing of a human fecal sample and a neural network algorithm based on tetranucleotide frequency profiling to reconstruct microbial genomes and tested the same approach using rectal swabs from a second patient. The approach rapidly and readily detected the genome of Klebsiella pneumoniae carbapenemase (KPC)-producing K. pneumoniae in the patient fecal specimen and in the rectal swab sample, achieving a level of strain resolution that was sufficient for confident transmission inference during a highly clonal outbreak. The analysis also detected previously unrecognized colonization of the patient by vancomycin-resistant Enterococcus faecium, another multidrug-resistant bacterium.IMPORTANCE The study results reported here perfectly demonstrate the power and promise of clinical metagenomics to recover genome sequences of important drug-resistant bacteria and to rapidly provide rich data that inform outbreak investigations and treatment decisions, independently of the need to culture the organisms.

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

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.

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

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

More evolution underground: Accelerated mitochondrial substitution rate in Australian burrowing freshwater crayfishes (Decapoda: Parastacidae).

To further understand the evolutionary history and mitogenomic features of Australia's highly distinctive freshwater crayfish fauna, we utilized a recently described rapid mitogenome sequencing pipeline to generate 24 new crayfish mitogenomes including a diversity of burrowing crayfish species and the first for Astacopsis gouldi, the world's largest freshwater invertebrate. Whole mitogenome-based phylogeny estimates using both Bayesian and Maximum Likelihood methods substantially strengthen existing hypotheses for systematic relationships among Australian freshwater crayfish with evidence of pervasive diversifying selection and accelerated mitochondrial substitution rate among the members of the clade representing strongly burrowing crayfish that may reflect selection pressures for increased energy requirement for adaptation to terrestrial environment and a burrowing lifestyle. Further, gene rearrangements are prevalent in the burrowing crayfish mitogenomes involving both tRNA and protein coding genes. In addition, duplicated control regions were observed in two closely related Engaeus species, together with evidence for concerted evolution. This study significantly adds to the understanding of Australian freshwater crayfish evolutionary relationships and suggests a link between mitogenome evolution and adaptation to terrestrial environments and a burrowing lifestyle in freshwater crayfish.

Topical Antibiotic Use Coselects for the Carriage of Mobile Genetic Elements Conferring Resistance to Unrelated Antimicrobials in Staphylococcus aureus.

Topical antibiotics, such as mupirocin and fusidic acid, are commonly used in the prevention and treatment of skin infections, particularly those caused by staphylococci. However, the widespread use of these agents is associated with increased resistance to these agents, potentially limiting their efficacy. Of particular concern is the observation that resistance to topical antibiotics is often associated with multidrug resistance, suggesting that topical antibiotics may play a role in the emergence of multidrug-resistant (MDR) strains. New Zealand (NZ) has some of the highest globally recorded rates of topical antibiotic usage and resistance. Using a combination of Pacific Biosciences single-molecule real-time (SMRT) whole-genome sequencing, Illumina short-read sequencing, and Bayesian phylogenomic modeling on 118 new multilocus sequence type 1 (ST1) community Staphylococcus aureus isolates from New Zealand and 61 publically available international ST1 genome sequences, we demonstrate a strong correlation between the clinical introduction of topical antibiotics and the emergence of MDR ST1 S. aureus We also provide in vitro experimental evidence showing that exposure to topical antibiotics can lead to the rapid selection of MDR S. aureus isolates carrying plasmids that confer resistance to multiple unrelated antibiotics, from within a mixed population of competitor strains. These findings have important implications regarding the impact of the indiscriminate use of topical antibiotics.

Donor-Derived Mycoplasma hominis and an Apparent Cluster of M. hominis Cases in Solid Organ Transplant Recipients.

BACKGROUND: Invasive and disseminated Mycoplasma hominis infections are well recognized but uncommon complications in solid organ transplant recipients. In a single center, a cluster of M. hominis infections were identified in lung transplant recipients from the same thoracic intensive care unit (ICU). We sought to determine the source(s) of these infections. METHODS: Medical records of the donor and infected transplant recipients were reviewed for clinical characteristics. Clinical specimens underwent routine processing with subculture on Mycoplasma-specific Hayflick agar. Mycoplasma hominis identification was confirmed using sequencing of the 16S ribosomal RNA gene. Mycoplasma hominis isolates were subjected to whole-genome sequencing on the Illumina NextSeq platform. RESULTS: Three lung transplant recipients presented with invasive M. hominis infections at multiple sites characterized by purulent infections without organisms detected by Gram staining. Each patient had a separate donor; however, pretransplant bronchoalveolar lavage fluid was only available from the donor for patient 1, which subsequently grew M. hominis. Phylo- and pangenomic analyses indicated that the isolates from the donor and the corresponding recipient (patient 1) were closely related and formed a distinct single clade. In contrast, isolates from patients 2 and 3 were unrelated and divergent from one another. CONCLUSIONS: Mycoplasma hominis should be considered a cause of donor-derived infection. Genomic data suggest donor-to-recipient transmission of M. hominis. Additional patients co-located in the ICU were found to have genetically unrelated M. hominis isolates, excluding patient-to-patient transmission.

A Supervised Statistical Learning Approach for Accurate Legionella pneumophila Source Attribution during Outbreaks.

Public health agencies are increasingly relying on genomics during Legionnaires' disease investigations. However, the causative bacterium (Legionella pneumophila) has an unusual population structure, with extreme temporal and spatial genome sequence conservation. Furthermore, Legionnaires' disease outbreaks can be caused by multiple L. pneumophila genotypes in a single source. These factors can confound cluster identification using standard phylogenomic methods. Here, we show that a statistical learning approach based on L. pneumophila core genome single nucleotide polymorphism (SNP) comparisons eliminates ambiguity for defining outbreak clusters and accurately predicts exposure sources for clinical cases. We illustrate the performance of our method by genome comparisons of 234 L. pneumophila isolates obtained from patients and cooling towers in Melbourne, Australia, between 1994 and 2014. This collection included one of the largest reported Legionnaires' disease outbreaks, which involved 125 cases at an aquarium. Using only sequence data from L. pneumophila cooling tower isolates and including all core genome variation, we built a multivariate model using discriminant analysis of principal components (DAPC) to find cooling tower-specific genomic signatures and then used it to predict the origin of clinical isolates. Model assignments were 93% congruent with epidemiological data, including the aquarium Legionnaires' disease outbreak and three other unrelated outbreak investigations. We applied the same approach to a recently described investigation of Legionnaires' disease within a UK hospital and observed a model predictive ability of 86%. We have developed a promising means to breach L. pneumophila genetic diversity extremes and provide objective source attribution data for outbreak investigations.IMPORTANCE Microbial outbreak investigations are moving to a paradigm where whole-genome sequencing and phylogenetic trees are used to support epidemiological investigations. It is critical that outbreak source predictions are accurate, particularly for pathogens, like Legionella pneumophila, which can spread widely and rapidly via cooling system aerosols, causing Legionnaires' disease. Here, by studying hundreds of Legionella pneumophila genomes collected over 21 years around a major Australian city, we uncovered limitations with the phylogenetic approach that could lead to a misidentification of outbreak sources. We implement instead a statistical learning technique that eliminates the ambiguity of inferring disease transmission from phylogenies. Our approach takes geolocation information and core genome variation from environmental L. pneumophila isolates to build statistical models that predict with high confidence the environmental source of clinical L. pneumophila during disease outbreaks. We show the versatility of the technique by applying it to unrelated Legionnaires' disease outbreaks in Australia and the UK.

Target-Specific Assay for Rapid and Quantitative Detection of Mycobacterium chimaera DNA.

Mycobacterium chimaera is an opportunistic environmental mycobacterium belonging to the Mycobacterium avium-M. intracellulare complex. Although most commonly associated with pulmonary disease, there has been growing awareness of invasive M. chimaera infections following cardiac surgery. Investigations suggest worldwide spread of a specific M. chimaera clone, associated with contaminated hospital heater-cooler units used during the surgery. Given the global dissemination of this clone, its potential to cause invasive disease, and the laboriousness of current culture-based diagnostic methods, there is a pressing need to develop rapid and accurate diagnostic assays specific for M. chimaera Here, we assessed 354 mycobacterial genome sequences and confirmed that M. chimaera is a phylogenetically coherent group. In silico comparisons indicated six DNA regions present only in M. chimaera We targeted one of these regions and developed a TaqMan quantitative PCR (qPCR) assay for M. chimaera with a detection limit of 100 CFU/ml in whole blood spiked with bacteria. In vitro screening against DNA extracted from 40 other mycobacterial species and 22 bacterial species from 21 diverse genera confirmed the in silico-predicted specificity for M. chimaera Screening 33 water samples from heater-cooler units with this assay highlighted the increased sensitivity of PCR compared to culture, with 15 of 23 culture-negative samples positive by M. chimaera qPCR. We have thus developed a robust molecular assay that can be readily and rapidly deployed to screen clinical and environmental specimens for M. chimaera.

Microbial mercury methylation in Antarctic sea ice.

Atmospheric deposition of mercury onto sea ice and circumpolar sea water provides mercury for microbial methylation, and contributes to the bioaccumulation of the potent neurotoxin methylmercury in the marine food web. Little is known about the abiotic and biotic controls on microbial mercury methylation in polar marine systems. However, mercury methylation is known to occur alongside photochemical and microbial mercury reduction and subsequent volatilization. Here, we combine mercury speciation measurements of total and methylated mercury with metagenomic analysis of whole-community microbial DNA from Antarctic snow, brine, sea ice and sea water to elucidate potential microbially mediated mercury methylation and volatilization pathways in polar marine environments. Our results identify the marine microaerophilic bacterium Nitrospina as a potential mercury methylator within sea ice. Anaerobic bacteria known to methylate mercury were notably absent from sea-ice metagenomes. We propose that Antarctic sea ice can harbour a microbial source of methylmercury in the Southern Ocean.

Genomic Insights into a Sustained National Outbreak of Yersinia pseudotuberculosis.

In 2014, a sustained outbreak of yersiniosis due to Yersinia pseudotuberculosis occurred across all major cities in New Zealand (NZ), with a total of 220 laboratory-confirmed cases, representing one of the largest ever reported outbreaks of Y. pseudotuberculosis. Here, we performed whole genome sequencing of outbreak-associated isolates to produce the largest population analysis to date of Y. pseudotuberculosis, giving us unprecedented capacity to understand the emergence and evolution of the outbreak clone. Multivariate analysis incorporating our genomic and clinical epidemiological data strongly suggested a single point-source contamination of the food chain, with subsequent nationwide distribution of contaminated produce. We additionally uncovered significant diversity in key determinants of virulence, which we speculate may help explain the high morbidity linked to this outbreak.