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.

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 ΦBT1 large serine recombinase catalyzes DNA integration at pseudo-attB sites in the genus Nocardia.

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.

Comparative Genomics Shows That Mycobacterium ulcerans Migration and Expansion Preceded the Rise of Buruli Ulcer in Southeastern Australia.

Since 2000, cases of the neglected tropical disease Buruli ulcer, caused by infection with Mycobacterium ulcerans, have increased 100-fold around Melbourne (population 4.4 million), the capital of Victoria, in temperate southeastern Australia. The reasons for this increase are unclear. Here, we used whole-genome sequence comparisons of 178 M. ulcerans isolates obtained primarily from human clinical specimens, spanning 70 years, to model the population dynamics of this pathogen from this region. Using phylogeographic and advanced Bayesian phylogenetic approaches, we found that there has been a migration of the pathogen from the east end of the state, beginning in the 1980s, 300 km west to the major human population center around Melbourne. This move was then followed by a significant increase in M. ulcerans population size. These analyses inform our thinking around Buruli ulcer transmission and control, indicating that M. ulcerans is introduced to a new environment and then expands, rather than it being from the awakening of a quiescent pathogen reservoir.IMPORTANCE Buruli ulcer is a destructive skin and soft tissue infection caused by Mycobacterium ulcerans and is characterized by progressive skin ulceration, which can lead to permanent disfigurement and long-term disability. Despite the majority of disease burden occurring in regions of West and central Africa, Buruli ulcer is also becoming increasingly common in southeastern Australia. Major impediments to controlling disease spread are incomplete understandings of the environmental reservoirs and modes of transmission of M. ulcerans The significance of our research is that we used genomics to assess the population structure of this pathogen at the Australian continental scale. We have then reconstructed a historical bacterial spread and modeled demographic dynamics to reveal bacterial population expansion across southeastern Australia. These findings provide explanations for the observed epidemiological trends with Buruli ulcer and suggest possible management to control disease spread.

Convergent Evolution Driven by Rifampin Exacerbates the Global Burden of Drug-Resistant Staphylococcus aureus.

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.

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.

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.

A nanomechanical study of the effects of colistin on the Klebsiella pneumoniae AJ218 capsule.

Atomic force microscopy measurements of capsule thickness revealed that that the wild-type Klebsiella pneumoniae AJ218 capsular polysaccharides were rearranged by exposure to colistin. The increase in capsule thickness measured near minimum inhibitory/bactericidal concentration (MIC/MBC) is consistent with the idea that colistin displaces the divalent cations that cross-bridge adjacent lipopolysaccharide (LPS) molecules through the capsule network. Cryo-electron microscopy demonstrated that the measured capsule thickness at near MIC/MBC of 1.2 µM was inflated by the disrupted outer membrane, through which the capsule is excreted and LPS is bound. Since wild-type and capsule-deficient strains of K. pneumoniae AJ218 have equivalent MICs and MBCs, the presence of the capsule appeared to confer no protection against colistin in AJ218. A spontaneously arising colistin mutant showed a tenfold increase in resistance to colistin; genetic analysis identified a single amino acid substitution (Q95P) in the PmrB sensor kinase in this colistin-resistant K. pneumoniae AJ218. Modification of the lipid A component of the LPS could result in a reduction of the net-negative charge of the outer membrane, which could hinder binding of colistin to the outer membrane and displacement of the divalent cations that bridge adjacent LPS molecules throughout the capsular polysaccharide network. Retention of the cross-linking divalent cations may explain why measurements of capsule thickness did not change significantly in the colistin-resistant strain after colistin exposure. These results contrast with those for other K. pneumoniae strains that suggest that the capsule confers colistin resistance.

Emergence of endemic MLST non-typeable vancomycin-resistant Enterococcus faecium.

BACKGROUND: Enterococcus faecium is a major nosocomial pathogen causing significant morbidity and mortality worldwide. Assessment of E. faecium using MLST to understand the spread of this organism is an important component of hospital infection control measures. Recent studies, however, suggest that MLST might be inadequate for E. faecium surveillance. OBJECTIVES: To use WGS to characterize recently identified vancomycin-resistant E. faecium (VREfm) isolates non-typeable by MLST that appear to be causing a multi-jurisdictional outbreak in Australia. METHODS: Illumina NextSeq and Pacific Biosciences SMRT sequencing platforms were used to determine the genome sequences of 66 non-typeable E. faecium (NTEfm) isolates. Phylogenetic and bioinformatics analyses were subsequently performed using a number of in silico tools. RESULTS: Sixty-six E. faecium isolates were identified by WGS from multiple health jurisdictions in Australia that could not be typed by MLST due to a missing pstS allele. SMRT sequencing and complete genome assembly revealed a large chromosomal rearrangement in representative strain DMG1500801, which likely facilitated the deletion of the pstS region. Phylogenomic analysis of this population suggests that deletion of pstS within E. faecium has arisen independently on at least three occasions. Importantly, the majority of these isolates displayed a vancomycin-resistant genotype. CONCLUSIONS: We have identified NTEfm isolates that appear to be causing a multi-jurisdictional outbreak in Australia. Identification of these isolates has important implications for MLST-based typing activities designed to monitor the spread of VREfm and provides further evidence supporting the use of WGS for hospital surveillance of E. faecium.

Prospective Whole-Genome Sequencing Enhances National Surveillance of Listeria monocytogenes.

Whole-genome sequencing (WGS) has emerged as a powerful tool for comparing bacterial isolates in outbreak detection and investigation. Here we demonstrate that WGS performed prospectively for national epidemiologic surveillance of Listeria monocytogenes has the capacity to be superior to our current approaches using pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), multilocus variable-number tandem-repeat analysis (MLVA), binary typing, and serotyping. Initially 423 L. monocytogenes isolates underwent WGS, and comparisons uncovered a diverse genetic population structure derived from three distinct lineages. MLST, binary typing, and serotyping results inferred in silico from the WGS data were highly concordant (>99%) with laboratory typing performed in parallel. However, WGS was able to identify distinct nested clusters within groups of isolates that were otherwise indistinguishable using our current typing methods. Routine WGS was then used for prospective epidemiologic surveillance on a further 97 L. monocytogenes isolates over a 12-month period, which provided a greater level of discrimination than that of conventional typing for inferring linkage to point source outbreaks. A risk-based alert system based on WGS similarity was used to inform epidemiologists required to act on the data. Our experience shows that WGS can be adopted for prospective L. monocytogenes surveillance and investigated for other pathogens relevant to public health.

Intercontinental dissemination of azithromycin-resistant shigellosis through sexual transmission: a cross-sectional study.

BACKGROUND: Shigellosis is an acute, severe bacterial colitis that, in high-income countries, is typically associated with travel to high-risk regions (Africa, Asia, and Latin America). Since the 1970s, shigellosis has also been reported as a sexually transmitted infection in men who have sex with men (MSM), in whom transmission is an important component of shigellosis epidemiology in high-income nations. We aimed to use sophisticated subtyping and international sampling to determine factors driving shigellosis emergence in MSM linked to an outbreak in the UK. METHODS: We did a large-scale, cross-sectional genomic epidemiological study of shigellosis cases collected from 29 countries between December, 1995, and June 8, 2014. Focusing on an ongoing epidemic in the UK, we collected and whole-genome sequenced clinical isolates of Shigella flexneri serotype 3a from high-risk and low-risk regions, including cases associated with travel and sex between men. We examined relationships between geographical, demographic, and clinical patient data with the isolate antimicrobial susceptibility, genetic data, and inferred evolutionary relationships. FINDINGS: We obtained 331 clinical isolates of S flexneri serotype 3a, including 275 from low-risk regions (44 from individuals who travelled to high-risk regions), 52 from high-risk regions, and four outgroup samples (ie, closely related, but genetically distinct isolates used to determine the root of the phylogenetic tree). We identified a recently emerged lineage of S flexneri 3a that has spread intercontinentally in less than 20 years throughout regions traditionally at low risk for shigellosis via sexual transmission in MSM. The lineage had acquired multiple antimicrobial resistance determinants, and prevailing sublineages were strongly associated with resistance to the macrolide azithromycin. Eight (4%) of 206 isolates from the MSM-associated lineage were obtained from patients who had previously provided an isolate; these serial isolations indicated atypical infection patterns (eg, reinfection). INTERPRETATION: We identified transmission-facilitating behaviours and atypical course(s) of infection as precipitating factors in shigellosis-affected MSM. The intercontinental spread of antimicrobial-resistant shigella through established transmission routes emphasises the need for new approaches to tackle the public health challenge of sexually transmitted infections in MSM. FUNDING: Wellcome Trust (grant number 098051).

SRST2: Rapid genomic surveillance for public health and hospital microbiology labs.

Rapid molecular typing of bacterial pathogens is critical for public health epidemiology, surveillance and infection control, yet routine use of whole genome sequencing (WGS) for these purposes poses significant challenges. Here we present SRST2, a read mapping-based tool for fast and accurate detection of genes, alleles and multi-locus sequence types (MLST) from WGS data. Using >900 genomes from common pathogens, we show SRST2 is highly accurate and outperforms assembly-based methods in terms of both gene detection and allele assignment. We include validation of SRST2 within a public health laboratory, and demonstrate its use for microbial genome surveillance in the hospital setting. In the face of rising threats of antimicrobial resistance and emerging virulence among bacterial pathogens, SRST2 represents a powerful tool for rapidly extracting clinically useful information from raw WGS data. Source code is available from http://katholt.github.io/srst2/.

Clinicopathological significance of expression of nestin, a neural stem/progenitor cell marker, in human glioma tissue.

The purpose of the study was to investigate the pathological and clinical significance of the expression of nestin, a type-VI intermediate filament transiently expressed during brain development, in glioma tissue. This study was conducted in 70 patients with newly diagnosed adult supratentorial gliomas who underwent multimodality treatment in our department, including surgery. The pathological diagnosis was grade II in 6 patients, grade III in 21 patients, and grade IV in 43 patients. Two specimen sections, one from the bulk of the removed tumor and one from the border between the tumor and normal brain tissue, were subjected to immunostaining with a mouse anti-human nestin monoclonal antibody. Analyses were performed to investigate possible correlation with pathological features, the relationship between nestin expression and the continuity of tumor with the subventricular zone (SVZ), correlation with the therapeutic prognosis, etc. Nestin was expressed specifically in astrocytoma lineage cells. In oligodendroglial tumors, nestin was expressed only in less-differentiated cells and cells suggestive of the presence of astrocytoma. In astrocytic tumors, the rate and level of nestin expression increased as the degree of malignancy increased. There was no significant correlation between the expression level of nestin and the continuity of tumor with the SVZ in the contrast-enhanced imaging before surgery. In addition, no correlation with the therapeutic prognosis was observed. Nestin, a neural stem cell marker, was specifically expressed in astrocytoma lineage cells in glioma tissue. A positive correlation was observed between the degree of malignancy and the level of nestin expression. However, the level of nestin expression was not related to the tumor localization in the SVZ and was not correlated with the therapeutic prognosis.

Identification of Streptococcus uberis multilocus sequence types highly associated with mastitis.

Multilocus sequence typing analysis of Streptococcus uberis has identified a cluster of isolates associated with clinical and subclinical mastitis and a cluster associated with cows with low somatic cell counts in their milk. Specific groups of genotypes (global clonal complex [GCC] sequence type 5s [ST5s] and GCC ST143s) were highly associated (P = 0.006) with clinical and subclinical mastitis and may represent a lineage of virulent isolates, whereas isolates belonging to GCC ST86 were associated with low-cell-count cows. This study has, for the first time, demonstrated the occurrence of identical sequence types (ST60 and ST184) between different continents (Australasia and Europe) and different countries (Australia and New Zealand). The standardized index of association and the empirical estimation of the rate of recombination showed substantial recombination within the S. uberis population in Australia, consistent with previous multilocus sequence type analyses.