Ryan R. Wick

Bandage: interactive visualization of de novo genome assemblies

Summary: Although de novo assembly graphs contain assembled contigs (nodes), the connections between those contigs (edges) are difficult for users to access. Bandage (a Bioinformatics Application for Navigating De novo Assembly Graphs Easily) is a tool for visualizing assembly graphs with connections. Users can zoom in to specific areas of the graph and interact with it by moving nodes, adding labels, changing colors and extracting sequences. BLAST searches can be performed within the Bandage graphical user interface and the hits are displayed as highlights in the graph. By displaying connections between contigs, Bandage presents new possibilities for analyzing de novo assemblies that are not possible through investigation of contigs alone. Availability and implementation: Source code and binaries are freely available at https://github.com/rrwick/Bandage. Bandage is implemented in C++ and supported on Linux, OS X and Windows. A full feature list and screenshots are available at http://rrwick.github.io/Bandage. Contact: rrwick@gmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads

The Illumina DNA sequencing platform generates accurate but short reads, which can be used to produce accurate but fragmented genome assemblies. Pacific Biosciences and Oxford Nanopore Technologies DNA sequencing platforms generate long reads that can produce complete genome assemblies, but the sequencing is more expensive and error-prone. There is significant interest in combining data from these complementary sequencing technologies to generate more accurate “hybrid” assemblies. However, few tools exist that truly leverage the benefits of both types of data, namely the accuracy of short reads and the structural resolving power of long reads. Here we present Unicycler, a new tool for assembling bacterial genomes from a combination of short and long reads, which produces assemblies that are accurate, complete and cost-effective. Unicycler builds an initial assembly graph from short reads using the de novo assembler SPAdes and then simplifies the graph using information from short and long reads. Unicycler uses a novel semi-global aligner to align long reads to the assembly graph. Tests on both synthetic and real reads show Unicycler can assemble larger contigs with fewer misassemblies than other hybrid assemblers, even when long-read depth and accuracy are low. Unicycler is open source (GPLv3) and available at github.com/rrwick/Unicycler.

Identification of Klebsiella capsule synthesis loci from whole genome data

Klebsiella pneumoniae is a growing cause of healthcare-associated infections for which multi-drug resistance is a concern. Its polysaccharide capsule is a major virulence determinant and epidemiological marker. However, little is known about capsule epidemiology since serological typing is not widely accessible and many isolates are serologically non-typeable. Molecular typing techniques provide useful insights, but existing methods fail to take full advantage of the information in whole genome sequences. We investigated the diversity of the capsule synthesis loci (K-loci) among 2503 K. pneumoniae genomes. We incorporated analyses of full-length K-locus nucleotide sequences and also clustered protein-encoding sequences to identify, annotate and compare K-locus structures. We propose a standardized nomenclature for K-loci and present a curated reference database. A total of 134 distinct K-loci were identified, including 31 novel types. Comparative analyses indicated 508 unique protein-encoding gene clusters that appear to reassort via homologous recombination. Extensive intra- and inter-locus nucleotide diversity was detected among the wzi and wzc genes, indicating that current molecular typing schemes based on these genes are inadequate. As a solution, we introduce Kaptive, a novel software tool that automates the process of identifying K-loci based on full locus information extracted from whole genome sequences (https://github.com/katholt/Kaptive). This work highlights the extensive diversity of Klebsiella K-loci and the proteins that they encode. The nomenclature, reference database and novel typing method presented here will become essential resources for genomic surveillance and epidemiological investigations of this pathogen.

Inducible colistin resistance via a disrupted plasmid-borne mcr-1 gene in a 2008 Vietnamese Shigella sonnei isolate

OBJECTIVES The objective of this study was to assess the presence of mcr-1 in Shigella sonnei isolated in Vietnam. METHODS WGS data were analysed for the presence of the mcr-1 gene sequence. The association of mcr-1 with a plasmid was assessed by PCR and by conjugation. RESULTS Through genome sequencing we identified a plasmid-associated inactive form of mcr-1 in a 2008 Vietnamese isolate of Shigella sonnei. The plasmid was conjugated into Escherichia coli and mcr-1 was activated upon exposure to colistin, resulting in highly colistin-resistant transconjugants. CONCLUSIONS This is the first description of the mcr-1 gene in Shigella, which is atypical given that colistin is not ordinarily used to treat diarrhoea. Our data suggest the mcr-1 gene has been circulating in human-restricted pathogens for some time but likely carries a selective fitness cost.

Gastrointestinal Carriage Is a Major Reservoir of Klebsiella pneumoniae Infection in Intensive Care Patients

Summary Klebsiella pneumoniae colonization is a significant risk factor for infection in ICU, with approximately half of K. pneumoniae infections resulting from patients’ own microbiota. Screening for colonization on admission could limit risk of infection in the colonized patient and others.

Completing bacterial genome assemblies with multiplex MinION sequencing

Illumina sequencing platforms have enabled widespread bacterial whole genome sequencing. While Illumina data is appropriate for many analyses, its short read length limits its ability to resolve genomic structure. This has major implications for tracking the spread of mobile genetic elements, including those which carry antimicrobial resistance determinants. Fully resolving a bacterial genome requires long-read sequencing such as those generated by Oxford Nanopore Technologies (ONT) platforms. Here we describe our use of the ONT MinION to sequence 12 isolates of Klebsiella pneumoniae on a single flow cell. We assembled each genome using a combination of ONT reads and previously available Illumina reads, and little to no manual intervention was needed to achieve fully resolved assemblies using the Unicycler hybrid assembler. Assembling only ONT reads with Canu was less effective, resulting in fewer resolved genomes and higher error rates even following error correction with Nanopolish. We demonstrate that multiplexed ONT sequencing is a valuable tool for high-throughput bacterial genome finishing. Specifically, we advocate the use of Illumina sequencing as a first analysis step, followed by ONT reads as needed to resolve genomic structure.

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.

Repeated local emergence of carbapenem-resistant Acinetobacter baumannii in a single hospital ward

We recently reported a dramatic increase in the prevalence of carbapenem-resistant Acinetobacter baumannii infections in the intensive care unit (ICU) of a Vietnamese hospital. This upsurge was associated with a specific oxa23-positive clone that was identified by multilocus VNTR analysis. Here, we used whole-genome sequence analysis to dissect the emergence of carbapenem-resistant A. baumannii causing ventilator-associated pneumonia (VAP) in the ICU during 2009–2012. To provide historical context and distinguish microevolution from strain introduction, we compared these genomes with those of A. baumannii asymptomatic carriage and VAP isolates from this same ICU collected during 2003–2007. We identified diverse lineages co-circulating over many years. Carbapenem resistance was associated with the presence of oxa23, oxa40, oxa58 and ndm1 genes in multiple lineages. The majority of resistant isolates were oxa23-positive global clone GC2; fine-scale phylogenomic analysis revealed five distinct GC2 sublineages within the ICU that had evolved locally via independent chromosomal insertions of oxa23 transposons. The increase in infections caused by carbapenem-resistant A. baumannii was associated with transposon-mediated transmission of a carbapenemase gene, rather than clonal expansion or spread of a carbapenemase-harbouring plasmid. Additionally, we found evidence of homologous recombination creating diversity within the local GC2 population, including several events resulting in replacement of the capsule locus. We identified likely donors of the imported capsule locus sequences amongst the A. baumannii isolated on the same ward, suggesting that diversification was largely facilitated via reassortment and sharing of genetic material within the localized A. baumannii population.

Frequent emergence of pathogenic lineages of Klebsiella pneumoniae via mobilisation of yersiniabactin and colibactin

Mobile genetic elements (MGEs) that frequently transfer within and between bacterial species play a critical role in bacterial evolution, and often carry key accessory genes that associate with a bacteria9s ability to cause disease. MGEs carrying antimicrobial resistance (AMR) and/or virulence determinants are common in opportunistic pathogen Klebsiella pneumoniae, which are a leading cause of highly drug-resistant infections in hospitals. Well-characterised virulence determinants in K. pneumoniae include the polyketide synthesis loci ybt and clb (also known as pks), encoding the iron-scavenging siderophore yersiniabactin and genotoxin colibactin respectively. These loci are located within an MGE called ICEKp, which is the most common virulence-associated MGE of K. pneumoniae, providing a mechanism for these virulence factors to spread within the population. Here we apply population genomics to investigate the prevalence, evolution and mobility of ybt and clb in K. pneumoniae populations through comparative analysis of 2,498 whole genome sequences. The ybt locus was detected in 40% of K. pneumoniae genomes, particularly amongst those associated with invasive infections. We identified 17 distinct ybt lineages and 3 clb lineages, each associated with one of 14 different structural variants of ICEKp. Comparison with the wider Enterobacteriaceae population showed occasional ICEKp acquisition by other members. The clb locus was present in 14% of all K. pneumoniae and 38.4% of ybt+ genomes. Hundreds of independent ICEKp integration events were detected affecting hundreds of phylogenetically distinct K. pneumoniae lineages, including ≥19 in the globally-disseminated carbapenem-resistant clone CG258. A novel plasmid-encoded form of ybt was also identified, representing a new mechanism for ybt dispersal in K. pneumoniae populations. These data show that MGEs carrying ybt and clb circulate freely in the K. pneumoniae population, including among multidrug-resistant strains, and should be considered a target for genomic surveillance along with AMR determinants.Klebsiella pneumoniae (Kp) is a commensal bacterium that causes opportunistic infections. Evidence is mounting that Kp strains carrying acquired siderophores (yersiniabactin, salmochelin and aerobactin) and/or the genotoxin colibactin are highly pathogenic and can cause invasive disease. Here we explored the diversity of the Kp integrative conjugative element (ICEKp), which mobilises the yersiniabactin locus ybt, by comparing 2499 diverse Kp genomes. We identified 17 distinct ybt lineages and 14 ICEKp structural variants (some of which carry colibactin (clb) or salmochelin synthesis loci). Hundreds of ICEKp transmission events were detected affecting hundreds of Kp lineages, including nearly >20 transfers into the globally-disseminated, carbapenem-resistant clonal group CG258. Additionally, we identify a plasmid-encoded lineage of ybt, representing a new mechanism for ybt dispersal in Kp populations. We introduce a novel sequence-based typing approach for identifying ybt and clb variants, to aid the identification of emerging pathogenic lineages and the convergence of antibiotic resistance and hypervirulence. SIGNIFICANCE Klebsiella pneumoniae infections are increasingly difficult to treat with antibiotics. Some K. pneumoniae carry extra genes that allow them to synthesise yersiniabactin, an iron-scavenging molecule, which enhances their ability to cause disease. These genes are located on a genetic element that can easily transfer between strains. Here, we screened 2499 K. pneumoniae genome sequences and found substantial diversity in the yersiniabactin genes and the associated genetic elements, including a novel mechanism of transfer, and detected hundreds of distinct yersiniabactin acquisition events between K. pneumoniae strains. We also developed tools to identify and type yersiniabactin genes, to help track the evolution and spread of yersiniabactin in global K. pneumoniae populations and to monitor for acquisition of yersiniabactin in antibiotic-resistant strains.

Genetic diversity, mobilisation and spread of the yersiniabactin-encoding mobile element ICEKp in Klebsiella pneumoniae populations.

Mobile genetic elements (MGEs) that frequently transfer within and between bacterial species play a critical role in bacterial evolution, and often carry key accessory genes that associate with a bacteria’s ability to cause disease. MGEs carrying antimicrobial resistance (AMR) and/or virulence determinants are common in the opportunistic pathogen Klebsiella pneumoniae, which is a leading cause of highly drug-resistant infections in hospitals. Well-characterised virulence determinants in K. pneumoniae include the polyketide synthesis loci ybt and clb (also known as pks), encoding the iron-scavenging siderophore yersiniabactin and genotoxin colibactin, respectively. These loci are located within an MGE called ICEKp, which is the most common virulence-associated MGE of K. pneumoniae, providing a mechanism for these virulence factors to spread within the population. Here we apply population genomics to investigate the prevalence, evolution and mobility of ybt and clb in K. pneumoniae populations through comparative analysis of 2498 whole-genome sequences. The ybt locus was detected in 40 % of K. pneumoniae genomes, particularly amongst those associated with invasive infections. We identified 17 distinct ybt lineages and 3 clb lineages, each associated with one of 14 different structural variants of ICEKp. Comparison with the wider population of the family Enterobacteriaceae revealed occasional ICEKp acquisition by other members. The clb locus was present in 14 % of all K. pneumoniae and 38.4 % of ybt+ genomes. Hundreds of independent ICEKp integration events were detected affecting hundreds of phylogenetically distinct K. pneumoniae lineages, including at least 19 in the globally-disseminated carbapenem-resistant clone CG258. A novel plasmid-encoded form of ybt was also identified, representing a new mechanism for ybt dispersal in K. pneumoniae populations. These data indicate that MGEs carrying ybt and clb circulate freely in the K. pneumoniae population, including among multidrug-resistant strains, and should be considered a target for genomic surveillance along with AMR determinants.