Search

Plasmid vectors based on bacteriophage integrases are important tools in molecular microbiology for the introduction of foreign DNA, especially into bacterial species where other systems for genetic manipulation are limited. Site specific integrases catalyze recombination between phage and bacterial attachment sites (attP and attB, respectively) and the best studied integrases in the actinomycetes are the serine integrases from the Streptomyces bacteriophages ΦC31 and ΦBT1. As this reaction is unidirectional and highly stable, vectors containing phage integrase systems have been used in a number of genetic engineering applications. Plasmids bearing the ΦBT1 integrase have been used to introduce DNA into Streptomyces and Amycolatopsis strains; however, they have not been widely studied in other actinobacterial genera. Here, we show that vectors based on ΦBT1 integrase can stably integrate into the chromosomes of a range of Nocardia species, and that this integration occurs despite the absence of canonical attB sites in these genomes. Furthermore, we show that a ΦBT1 integrase-based vector can insert at multiple pseudo-attB sites within a single strain and we determine the sequence of a pseudo-attB motif. These data suggest that ΦBT1 integrase-based vectors can be used to readily and semi-randomly introduce foreign DNA into the genomes of a range of Nocardia species. However, the precise site of insertion will likely require empirical determination in each species to avoid unexpected off-target effects.

Abstract Serial intervals – the time between symptom onset in infector and infectee – are a fundamental quantity in infectious disease control. However, their estimation requires knowledge of individuals’ exposures, typically obtained through resource-intensive contact tracing efforts. We introduce an alternate framework using virus sequences to inform who infected whom and thereby estimate serial intervals. We apply our technique to SARS-CoV-2 sequences from case clusters in the first two COVID-19 waves in Victoria, Australia. We find that our approach offers high resolution, cluster-specific serial interval estimates that are comparable with those obtained from contact data, despite requiring no knowledge of who infected whom and relying on incompletely-sampled data. Compared to a published serial interval, cluster-specific serial intervals can vary estimates of the effective reproduction number by a factor of 2–3. We find that serial interval estimates in settings such as schools and meat processing/packing plants are shorter than those in healthcare facilities.

Abstract Background The COVID-19 pandemic remains a global public health concern. Advances in sequencing technologies has allowed for high numbers of SARS-CoV-2 whole genome sequence (WGS) data and rapid sharing of sequences through global repositories to enable almost real-time genomic analysis of the pathogen. WGS data has been used previously to group genetically similar viral pathogens to reveal evidence of transmission, including methods that identify distinct clusters on a phylogenetic tree. Identifying clusters of linked cases can aid in the regional surveillance and management of the disease. In this study, we present a novel method for producing stable genomic clusters of SARS-CoV-2 cases, cov2clusters, and compare the accuracy and stability of our approach to previous methods used for phylogenetic clustering using real-world SARS-CoV-2 sequence data obtained from British Columbia, Canada. Results We found that cov2clusters produced more stable clusters than previously used phylogenetic clustering methods when adding sequence data through time, mimicking an increase in sequence data through the pandemic. Our method also showed high accuracy when predicting epidemiologically informed clusters from sequence data. Conclusions Our new approach allows for the identification of stable clusters of SARS-CoV-2 from WGS data. Producing high-resolution SARS-CoV-2 clusters from sequence data alone can a challenge and, where possible, both genomic and epidemiological data should be used in combination.

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.

Salmonella enterica serovar 4,[5],12:i:- (Salmonella 4,[5],12:i:-) is a major pathogen of humans and animals with a reported incidence in Australia three times higher than the UK and USA. Here, the authors report the circulation, antimicrobial resistance signatures, and effects on host cells, of three Salmonella4,[5],12:i:- lineages within Australia.

Abstract Background Pathogen whole genome sequencing (WGS) is being incorporated into public health surveillance and disease control systems worldwide and has the potential to make significant contributions to infectious disease surveillance, outbreak investigation and infection prevention and control. However, to date, there are limited data regarding (i) the optimal models for integration of genomic data into epidemiological investigations and (ii) how to quantify and evaluate public health impacts resulting from genomic epidemiological investigations. Methods We developed the Pathogen Genomics in Public HeAlth Surveillance Evaluation (PG-PHASE) Framework to guide examination of the use of WGS in public health surveillance and disease control. We illustrate the use of this framework with three pathogens as case studies: Listeria monocytogenes, Mycobacterium tuberculosis and SARS-CoV-2. Results The framework utilises an adaptable whole-of-system approach towards understanding how interconnected elements in the public health application of pathogen genomics contribute to public health processes and outcomes. The three phases of the PG-PHASE Framework are designed to support understanding of WGS laboratory processes, analysis, reporting and data sharing, and how genomic data are utilised in public health practice across all stages, from the decision to send an isolate or sample for sequencing to the use of sequence data in public health surveillance, investigation and decision-making. Importantly, the phases can be used separately or in conjunction, depending on the need of the evaluator. Subsequent to conducting evaluation underpinned by the framework, avenues may be developed for strategic investment or interventions to improve utilisation of whole genome sequencing. Conclusions Comprehensive evaluation is critical to support health departments, public health laboratories and other stakeholders to successfully incorporate microbial genomics into public health practice. The PG-PHASE Framework aims to assist public health laboratories, health departments and authorities who are either considering transitioning to whole genome sequencing or intending to assess the integration of WGS in public health practice, including the capacity to detect and respond to outbreaks and associated costs, challenges and facilitators in the utilisation of microbial genomics and public health impacts.

Genome sequencing can be used to infer pathogen transmission dynamics and inform public health responses. Here, the authors sequence >1,200 SARS-CoV-2 samples from Victoria, Australia and find genomic support for the effectiveness of social restrictions in reducing transmission.

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

Summary: 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.

Malarial and other haemosporidian parasites are widespread; however, their temporal dynamics are ill-understood. Longitudinal sampling of a threatened riparian bird revealed a consistently very low prevalence over 13 years (∼5%) despite infections persisting and prevalence increasing with age. In contrast, three key species within this tropical community were highly infected (∼20–75% prevalence) and these differences were stable. Although we found novel lineages and phylogenetic structure at the local level, there was little geographic structuring within Australasia. This study suggests that malarial parasite susceptibility is determined by host factors and that species can maintain low levels despite high community prevalence. Keywords: Avian malaria, Wildlife diseases, Purple-crowned fairy-wrens, Haemoproteus, Plasmodium

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

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

ABSTRACT Staphylococcus epidermidis is a significant opportunistic pathogen of humans. Molecular studies in this species have been hampered by the presence of restriction-modification (RM) systems that limit introduction of foreign DNA. Here, we establish the complete genomes and methylomes for seven clinically significant, genetically diverse S. epidermidis isolates and perform the first systematic genomic analyses of the type I RM systems within both S. epidermidis and Staphylococcus aureus. Our analyses revealed marked differences in the gene arrangement, chromosomal location, and movement of type I RM systems between the two species. Unlike S. aureus, S. epidermidis type I RM systems demonstrate extensive diversity even within a single genetic lineage. This is contrary to current assumptions and has important implications for approaching the genetic manipulation of S. epidermidis. Using Escherichia coli plasmid artificial modification (PAM) to express S. epidermidis hsdMS, we readily overcame restriction barriers in S. epidermidis and achieved electroporation efficiencies equivalent to those of modification-deficient mutants. With these functional experiments, we demonstrated how genomic data can be used to predict both the functionality of type I RM systems and the potential for a strain to be electroporation proficient. We outline an efficient approach for the genetic manipulation of S. epidermidis strains from diverse genetic backgrounds, including those that have hitherto been intractable. Additionally, we identified S. epidermidis BPH0736, a naturally restriction-defective, clinically significant, multidrug-resistant ST2 isolate, as an ideal candidate for molecular studies. IMPORTANCE Staphylococcus epidermidis is a major cause of hospital-acquired infections, especially those related to implanted medical devices. Understanding how S. epidermidis causes disease and devising ways to combat these infections have been hindered by an inability to genetically manipulate clinically significant hospital-adapted strains. Here, we provide the first comprehensive analyses of the barriers to the uptake of foreign DNA in S. epidermidis and demonstrate that these are distinct from those described for S. aureus. Using these insights, we demonstrate an efficient approach for the genetic manipulation of S. epidermidis to enable the study of clinical isolates for the first time.

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

ABSTRACT Mutations in the beta-subunit of bacterial RNA polymerase (RpoB) cause resistance to rifampin (Rifr), a critical antibiotic for treatment of multidrug-resistant Staphylococcus aureus. In vitro studies have shown that RpoB mutations confer decreased susceptibility to other antibiotics, but the clinical relevance is unknown. Here, by analyzing 7,099 S. aureus genomes, we demonstrate that the most prevalent RpoB mutations promote clinically relevant phenotypic plasticity resulting in the emergence of stable S. aureus lineages, associated with increased risk of therapeutic failure through generation of small-colony variants (SCVs) and coresistance to last-line antimicrobial agents. We found eight RpoB mutations that accounted for 93% (469/505) of the total number of Rifr mutations. The most frequently selected amino acid substitutions affecting residue 481 (H481N/Y) were associated with worldwide expansions of Rifr clones spanning decades. Recreating the H481N/Y mutations confirmed no impact on S. aureus growth, but the H481N mutation promoted the emergence of a subpopulation of stable Rifr SCVs with reduced susceptibility to vancomycin and daptomycin. Recreating the other frequent RpoB mutations showed similar impacts on resistance to these last-line agents. We found that 86% of all Rifr isolates in our global sample carried the mutations promoting cross-resistance to vancomycin and 52% to both vancomycin and daptomycin. As four of the most frequent RpoB mutations confer only low-level Rifr, equal to or below some international breakpoints, we recommend decreasing these breakpoints and reconsidering the appropriate use of rifampin to reduce the fixation and spread of these clinically deleterious mutations. IMPORTANCE Increasing antibiotic resistance in the major human pathogen Staphylococcus aureus is threatening the ability to treat patients with these infections. Recent laboratory studies suggest that mutations in the gene commonly associated with rifampin resistance may also impact susceptibility to other last-line antibiotics in S. aureus; however, the overall frequency and clinical impact of these mutations are unknown. By mining a global collection of clinical S. aureus genomes and by mutagenesis experiments, this work reveals that common rifampin-induced rpoB mutations promote phenotypic plasticity that has led to the global emergence of stable, multidrug-resistant S. aureus lineages that are associated with increased risk of therapeutic failure through coresistance to other last-line antimicrobials. We recommend decreasing susceptibility breakpoints for rifampin to allow phenotypic detection of critical rpoB mutations conferring low resistance to rifampin and reconsidering the appropriate use of rifampin to reduce the fixation and spread of these deleterious mutations globally.

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.

Background The emergence and evolution of community-acquired methicillin resistant Staphylococcus aureus (CA-MRSA) strains in Africa is poorly understood. However, one particular MRSA lineage called ST88, appears to be rapidly establishing itself as an “African” CA-MRSA clone. In this study, we employed whole genome sequencing to provide more information on the genetic background of ST88 CA-MRSA isolates from Ghana and to describe in detail ST88 CA-MRSA isolates in comparison with other MRSA lineages worldwide. Methods We first established a complete ST88 reference genome (AUS0325) using PacBio SMRT sequencing. We then used comparative genomics to assess relatedness among 17 ST88 CA-MRSA isolates recovered from patients attending Buruli ulcer treatment centres in Ghana, three non-African ST88s and 15 other MRSA lineages. Results We show that Ghanaian ST88 forms a discrete MRSA lineage (harbouring SCCmec-IV [2B]). Gene content analysis identified five distinct genomic regions enriched among ST88 isolates compared with the other S. aureus lineages. The Ghanaian ST88 isolates had only 658 core genome SNPs and there was no correlation between phylogeny and geography, suggesting the recent spread of this clone. The lineage was also resistant to multiple classes of antibiotics including β-lactams, tetracycline and chloramphenicol. Discussion This study reveals that S. aureus ST88-IV is a recently emerging and rapidly spreading CA-MRSA clone in Ghana. The study highlights the capacity of small snapshot genomic studies to provide actionable public health information in resource limited settings. To our knowledge this is the first genomic assessment of the ST88 CA-MRSA clone.

We tested the hypothesis that tapirs tolerate individuals from adjacent and overlapping home ranges if they are related. We obtained genetic data from fecal samples collected in the Balbina reservoir landscape, central Amazon. Samples were genotyped at 14 microsatellite loci, of which five produced high quality informative genotypes. Based on an analysis of 32 individuals, we inferred a single panmictic population with high levels of heterozygosity. Kinship analysis identified 10 pairs of full siblings or parent-offspring, 10 pairs of half siblings and 25 unrelated pairs. In 10 cases, the related individuals were situated on opposite margins of the reservoir, suggesting that tapirs are capable of crossing the main river, even after damming. The polygamous model was the most likely mating system for Tapirus terrestris. Moran's I index of allele sharing between pairs of individuals geographically close (3 km). Confirming this result, the related individuals were not geographically closer than unrelated ones (W = 188.5; p = 0.339). Thus, we found no evidence of a preference for being close to relatives and observed a tendency for dispersal. The small importance of relatedness in determining spatial distribution of individuals is unusual in mammals, but not unheard of. Finally, non-invasive sampling allowed efficient access to the genetic data, despite the warm and humid climate of the Amazon, which accelerates DNA degradation.