Mark B. Schultz

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

Genomic analysis of diversity, population structure, virulence, and antimicrobial resistance in Klebsiella pneumoniae, an urgent threat to public health

Significance Klebsiella pneumoniae is rapidly becoming untreatable using last-line antibiotics. It is especially problematic in hospitals, where it causes a range of acute infections. To approach controlling such a bacterium, we first must define what it is and how it varies genetically. Here we have determined the DNA sequence of K. pneumoniae isolates from around the world and present a detailed analysis of these data. We show that there is a wide spectrum of diversity, including variation within shared sequences and gain and loss of whole genes. Using this detailed blueprint, we show that there is an unrecognized association between the possession of specific gene profiles associated with virulence and antibiotic resistance and the differing disease outcomes seen for K. pneumoniae. Klebsiella pneumoniae is now recognized as an urgent threat to human health because of the emergence of multidrug-resistant strains associated with hospital outbreaks and hypervirulent strains associated with severe community-acquired infections. K. pneumoniae is ubiquitous in the environment and can colonize and infect both plants and animals. However, little is known about the population structure of K. pneumoniae, so it is difficult to recognize or understand the emergence of clinically important clones within this highly genetically diverse species. Here we present a detailed genomic framework for K. pneumoniae based on whole-genome sequencing of more than 300 human and animal isolates spanning four continents. Our data provide genome-wide support for the splitting of K. pneumoniae into three distinct species, KpI (K. pneumoniae), KpII (K. quasipneumoniae), and KpIII (K. variicola). Further, for K. pneumoniae (KpI), the entity most frequently associated with human infection, we show the existence of >150 deeply branching lineages including numerous multidrug-resistant or hypervirulent clones. We show K. pneumoniae has a large accessory genome approaching 30,000 protein-coding genes, including a number of virulence functions that are significantly associated with invasive community-acquired disease in humans. In our dataset, antimicrobial resistance genes were common among human carriage isolates and hospital-acquired infections, which generally lacked the genes associated with invasive disease. The convergence of virulence and resistance genes potentially could lead to the emergence of untreatable invasive K. pneumoniae infections; our data provide the whole-genome framework against which to track the emergence of such threats.

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.

Integrated shotgun sequencing and bioinformatics pipeline allows ultra-fast mitogenome recovery and confirms substantial gene rearrangements in Australian freshwater crayfishes

BackgroundAlthough it is possible to recover the complete mitogenome directly from shotgun sequencing data, currently reported methods and pipelines are still relatively time consuming and costly. Using a sample of the Australian freshwater crayfish Engaeus lengana, we demonstrate that it is possible to achieve three-day turnaround time (four hours hands-on time) from tissue sample to NCBI-ready submission file through the integration of MiSeq sequencing platform, Nextera sample preparation protocol, MITObim assembly algorithm and MITOS annotation pipeline.ResultsThe complete mitochondrial genome of the parastacid freshwater crayfish, Engaeus lengana, was recovered by modest shotgun sequencing (1.2 giga bases) using the Illumina MiSeq benchtop sequencing platform. Genome assembly using the MITObim mitogenome assembler recovered the mitochondrial genome as a single contig with a 97-fold mean coverage (min. = 17; max. = 138). The mitogenome consists of 15,934 base pairs and contains the typical 37 mitochondrial genes and a non-coding AT-rich region. The genome arrangement is similar to the only other published parastacid mitogenome from the Australian genus Cherax.ConclusionsWe infer that the gene order arrangement found in Cherax destructor is common to Australian crayfish and may be a derived feature of the southern hemisphere family Parastacidae. Further, we report to our knowledge, the simplest and fastest protocol for the recovery and assembly of complete mitochondrial genomes using the MiSeq benchtop sequencer.

Convergent Adaptation in the Dominant Global Hospital Clone ST239 of Methicillin-Resistant Staphylococcus aureus

ABSTRACT Infections caused by highly successful clones of hospital-associated methicillin-resistant Staphylococcus aureus (HA-MRSA) are a major public health burden. The globally dominant sequence type 239 (ST239) HA-MRSA clone has persisted in the health care setting for decades, but the basis of its success has not been identified. Taking a collection of 123 ST239 isolates spanning 32 years, we have used population-based functional genomics to investigate the evolution of this highly persistent and successful clone. Phylogenetic reconstruction and population modeling uncovered a previously unrecognized distinct clade of ST239 that was introduced into Australia from Asia and has perpetuated the epidemic in this region. Functional analysis demonstrated attenuated virulence and enhanced resistance to last-line antimicrobials, the result of two different phenomena, adaptive evolution within the original Australian ST239 clade and the introduction of a new clade displaying shifts in both phenotypes. The genetic diversity between the clades allowed us to employ genome-wide association testing and identify mutations in other essential regulatory systems, including walKR, that significantly associate with and may explain these key phenotypes. The phenotypic convergence of two independently evolving ST239 clades highlights the very strong selective pressures acting on HA-MRSA, showing that hospital environments have favored the accumulation of mutations in essential MRSA genes that increase resistance to antimicrobials, attenuate virulence, and promote persistence in the health care environment. Combinations of comparative genomics and careful phenotypic measurements of longitudinal collections of clinical isolates are giving us the knowledge to intelligently address the impact of current and future antibiotic usage policies and practices on hospital pathogens globally. IMPORTANCE Methicillin-resistant Staphylococcus aureus (MRSA) is responsible for innumerable drug-resistant health care-associated infections globally. This study, the first to investigate the evolutionary response of hospital-associated MRSA (HA-MRSA) over many decades, demonstrates how MRSA can persist in a region through the reintroduction of a previously unrecognized distinct clade. This study also demonstrates the crucial adaptive responses of HA-MRSA to the highly selective environment of the health care system, the evolution of MRSA isolates to even higher levels of antibiotic resistance at the cost of attenuated virulence. However, in vivo persistence is maintained, resulting in a clone of HA-MRSA able to resist almost all antimicrobial agents and still cause invasive disease in the heavily compromised hosts found in modern health care settings. Methicillin-resistant Staphylococcus aureus (MRSA) is responsible for innumerable drug-resistant health care-associated infections globally. This study, the first to investigate the evolutionary response of hospital-associated MRSA (HA-MRSA) over many decades, demonstrates how MRSA can persist in a region through the reintroduction of a previously unrecognized distinct clade. This study also demonstrates the crucial adaptive responses of HA-MRSA to the highly selective environment of the health care system, the evolution of MRSA isolates to even higher levels of antibiotic resistance at the cost of attenuated virulence. However, in vivo persistence is maintained, resulting in a clone of HA-MRSA able to resist almost all antimicrobial agents and still cause invasive disease in the heavily compromised hosts found in modern health care settings.

DNA evidence of whale sharks (Rhincodon typus) feeding on red crab (Gecarcoidea natalis) larvae at Christmas Island, Australia

Whale sharks (Rhincodon typus) are thought to aggregate in nearshore waters around Christmas Island (105°37′E, 10o29′S) to consume the marine larvae of the endemic red land crab (Gecarcoidea natalis). However, there have been no direct observations of sharks feeding on crab larvae. Whale shark faeces were analysed using genetic testing to confirm the presence of crab larvae in their diet. Primers were designed for amplifying two Gecarcoidea natalis mitochondrial small-subunit (mtSSU) rDNA regions. Gel electrophoresis of polymerase chain reaction (PCR) products amplified from whale shark faecal DNA produced bands of the expected size for G. natalis templates. Specificity of both primer sets for G. natalis mtSSU rDNA was expected to be high from comparisons with mtSSU rDNA regions from closely related crabs and we confirmed their specificity empirically. The amplification of fragments from faecal DNA of the same size as those produced from G. natalis DNA indicates that the whale shark had been feeding on G. natalis and that enough of the crab DNA survived digestion to be detected by these PCRs. Our study provides further evidence that aggregations of whale sharks in coastal waters occur in response to ephemeral but predictable increases in planktonic prey.

Evolution underground: A molecular phylogenetic investigation of Australian burrowing freshwater crayfish (Decapoda: Parastacidae) with particular focus on Engaeus Erichson

Phylogenetic relationships and species boundaries of Australian burrowing freshwater crayfish belonging to the genera Engaeus, Engaewa, Geocharax, Gramastacus and Tenuibranchiurus are investigated using combined mitochondrial and nuclear DNA sequence data and Bayesian and Maximum Parsimony methods. Phylogenies are statistically compared to previously published hypotheses. Engaeus, Engaewa, Geocharax, Gramastacus and Tenuibranchiurus form a strongly supported monophyletic clade. This grouping is independently supported by morphology but unites geographically highly disjunct lineages. Our data show two cryptic species in Geocharax, one cryptic species in Gramastacus and two cryptic species within the highly divergent Engaeus lyelli lineage. Using a Bayesian relaxed molecular clock method, the 16S rDNA data show generic-level diversification coinciding with the transition from a wet to arid palaeoclimate near the mid Miocene.

Nuclear and Mitochondrial Patterns of Population Structure in North Pacific False Killer Whales (Pseudorca crassidens)

False killer whales (Pseudorca crassidens) are large delphinids typically found in deep water far offshore. However, in the Hawaiian Archipelago, there are 2 resident island-associated populations of false killer whales, one in the waters around the main Hawaiian Islands (MHI) and one in the waters around the Northwestern Hawaiian Islands (NWHI). We use mitochondrial DNA (mtDNA) control region sequences and genotypes from 16 nuclear DNA (nucDNA) microsatellite loci from 206 individuals to examine levels of differentiation among the 2 island-associated populations and offshore animals from the central and eastern North Pacific. Both mtDNA and nucDNA exhibit highly significant differentiation between populations, confirming limited gene flow in both sexes. The mtDNA haplotypes exhibit a strong pattern of phylogeographic concordance, with island-associated populations sharing 3 closely related haplotypes not found elsewhere in the Pacific. However, nucDNA data suggest that NWHI animals are at least as differentiated from MHI animals as they are from offshore animals. The patterns of differentiation revealed by the 2 marker types suggest that the island-associated false killer whale populations likely share a common colonization history, but have limited contemporary gene flow.

Five decades of genome evolution in the globally distributed, extensively antibiotic resistant Acinetobacter baumannii global clone 1

The majority of Acinetobacter baumannii isolates that are multiply, extensively and pan-antibiotic resistant belong to two globally disseminated clones, GC1 and GC2, that were first noticed in the 1970s. Here, we investigated microevolution and phylodynamics within GC1 via analysis of 45 whole-genome sequences, including 23 sequenced for this study. The most recent common ancestor of GC1 arose around 1960 and later diverged into two phylogenetically distinct lineages. In the 1970s, the main lineage acquired the AbaR resistance island, conferring resistance to older antibiotics, via a horizontal gene transfer event. We estimate a mutation rate of ∼5 SNPs genome− 1 year− 1 and detected extensive recombination within GC1 genomes, introducing nucleotide diversity into the population at >20 times the substitution rate (the ratio of SNPs introduced by recombination compared with mutation was 22). The recombination events were non-randomly distributed in the genome and created significant diversity within loci encoding outer surface molecules (including the capsular polysaccharide, the outer core lipooligosaccharide and the outer membrane protein CarO), and spread antimicrobial resistance-conferring mutations affecting the gyrA and parC genes and insertion sequence insertions activating the ampC gene. Both GC1 lineages accumulated resistance to newer antibiotics through various genetic mechanisms, including the acquisition of plasmids and transposons or mutations in chromosomal genes. Our data show that GC1 has diversified into multiple successful extensively antibiotic-resistant subclones that differ in their surface structures. This has important implications for all avenues of control, including epidemiological tracking, antimicrobial therapy and vaccination.

Cryptic diversity in Engaeus Erichson, Geocharax Clark and Gramastacus Riek (Decapoda:Parastacidae) revealed by mitochondrial 16S rDNA sequences

Nucleotide sequence data from the mitochondrial 16S rDNA region were utilised to investigate phylogenetic relationships and species boundaries among Australian freshwater crayfish belonging to the genera Engaeus Erichson, 1846, Geocharax Clark, 1936 and Gramastacus Riek, 1972. Geocharax and Gramastacus were found to be monophyletic genera but one species currently assigned to Engaeus may belong to another genus. Relationships between the three existing genera were not resolved. Analysis of species boundaries within Geocharax suggests that there are an additional two species in this genus, and our analysis of Gramastacus indicates that undescribed populations from central New South Wales may comprise a second species. The data provide at least one instance of a taxon crossing the Great Dividing Range and provide confirmation of previously proposed hypotheses seeking to explain trans-Bass Strait distributions of species.