Darrin Lemmer

Genomic Analysis of the Emergence and Rapid Global Dissemination of the Clonal Group 258 Klebsiella pneumoniae Pandemic

Multidrug-resistant Klebsiella pneumoniae producing the KPC carbapenemase have rapidly spread throughout the world, causing severe healthcare-associated infections with limited antimicrobial treatment options. Dissemination of KPC-producing K. pneumoniae is largely attributed to expansion of a single dominant strain, ST258. In this study, we explore phylogenetic relationships and evolution within ST258 and its clonal group, CG258, using whole genome sequence analysis of 167 isolates from 20 countries collected over 17 years. Our results show a common ST258 ancestor emerged from its diverse parental clonal group around 1995 and likely acquired bla KPC prior to dissemination. Over the past two decades, ST258 has remained highly clonal despite diversity in accessory elements and divergence in the capsule polysaccharide synthesis locus. Apart from the large recombination event that gave rise to ST258, few mutations set it apart from its clonal group. However, one mutation occurs in a global transcription regulator. Characterization of outer membrane protein sequences revealed a profile in ST258 that includes a truncated OmpK35 and modified OmpK37. Our work illuminates potential genomic contributors to the pathogenic success of ST258, helps us better understand the global dissemination of this strain, and identifies genetic markers unique to ST258.

Nasp: An accurate, rapid method for the identification of snps in wgs datasets that supports flexible input and output formats

Whole-genome sequencing (WGS) of bacterial isolates has become standard practice in many laboratories. Applications for WGS analysis include phylogeography and molecular epidemiology, using single nucleotide polymorphisms (SNPs) as the unit of evolution. NASP was developed as a reproducible method that scales well with the hundreds to thousands of WGS data typically used in comparative genomics applications. In this study, we demonstrate how NASP compares with other tools in the analysis of two real bacterial genomics datasets and one simulated dataset. Our results demonstrate that NASP produces similar, and often better, results in comparison with other pipelines, but is much more flexible in terms of data input types, job management systems, diversity of supported tools and output formats. We also demonstrate differences in results based on the choice of the reference genome and choice of inferring phylogenies from concatenated SNPs or alignments including monomorphic positions. NASP represents a source-available, version-controlled, unit-tested method and can be obtained from tgennorth.github.io/NASP.

The Northern Arizona SNP Pipeline (NASP): accurate, flexible, and rapid identification of SNPs in WGS datasets

Whole genome sequencing (WGS) of bacteria is becoming standard practice in many laboratories. Applications for WGS analysis include phylogeography and molecular epidemiology, using single nucleotide polymorphisms (SNPs) as the unit of evolution. The Northern Arizona SNP Pipeline (NASP) was developed as a reproducible pipeline that scales well with the large amount of WGS data typically used in comparative genomics applications. In this study, we demonstrate how NASP compares to other tools in the analysis of two real bacterial genomics datasets and one simulated dataset. Our results demonstrate that NASP produces comparable, and often better, results to other pipelines, but is much more flexible in terms of data input types, job management systems, diversity of supported tools, and output formats. We also demonstrate differences in results based on the choice of the reference genome and choice of inferring phylogenies from concatenated SNPs or alignments including monomorphic positions. NASP represents a source-available, version-controlled, unit-tested method and can be obtained from tgennorth.github.io/NASP.

Rapid Drug Susceptibility Testing of Drug-Resistant Mycobacterium tuberculosis Isolates Directly from Clinical Samples by Use of Amplicon Sequencing: a Proof-of-Concept Study

ABSTRACT Increasingly complex drug-resistant tuberculosis (DR-TB) is a major global health concern and one of the primary reasons why TB is now the leading infectious cause of death worldwide. Rapid characterization of a DR-TB patients complete drug resistance profile would facilitate individualized treatment in place of empirical treatment, improve treatment outcomes, prevent amplification of resistance, and reduce the transmission of DR-TB. The use of targeted next-generation sequencing (NGS) to obtain drug resistance profiles directly from patient sputum samples has the potential to enable comprehensive evidence-based treatment plans to be implemented quickly, rather than in weeks to months, which is currently needed for phenotypic drug susceptibility testing (DST) results. In this pilot study, we evaluated the performance of amplicon sequencing of Mycobacterium tuberculosis DNA from patient sputum samples using a tabletop NGS technology and automated data analysis to provide a rapid DST solution (the Next Gen-RDST assay). One hundred sixty-six out of 176 (94.3%) sputum samples from the Republic of Moldova yielded complete Next Gen-RDST assay profiles for 7 drugs of interest. We found a high level of concordance of our Next Gen-RDST assay results with phenotypic DST (97.0%) and pyrosequencing (97.8%) results from the same clinical samples. Our Next Gen-RDST assay was also able to estimate the proportion of resistant-to-wild-type alleles down to mixtures of ≤1%, which demonstrates the ability to detect very low levels of resistant variants not detected by pyrosequencing and possibly below the threshold for phenotypic growth methods. The assay as described here could be used as a clinical or surveillance tool.

Genome sequence of Staphylococcus aureus strain CA-347, a USA600 methicillin-resistant isolate

ABSTRACT The Staphylococcus aureus clonal lineage CC45 is a predominant colonizer of healthy individuals in northern Europe and constitutes a highly basal cluster of the S. aureus population. Here, we report the complete genome sequence of S. aureus strain CA-347 (NRS648), a representative of the methicillin-resistant USA600 clone predominantly found in the United States.

KlebSeq, a Diagnostic Tool for Surveillance, Detection, and Monitoring of Klebsiella pneumoniae

ABSTRACT Health care-acquired infections (HAIs) kill tens of thousands of people each year and add significantly to health care costs. Multidrug-resistant and epidemic strains are a large proportion of HAI agents, and multidrug-resistant strains of Klebsiella pneumoniae, a leading HAI agent, have caused an urgent public health crisis. In the health care environment, patient colonization by K. pneumoniae precedes infection, and transmission via colonization leads to outbreaks. Periodic patient screening for K. pneumoniae colonization has the potential to curb the number of HAIs. In this report, we describe the design and validation of KlebSeq, a highly informative screening tool that detects Klebsiella species and identifies clinically important strains and characteristics by using highly multiplexed amplicon sequencing without a live-culturing step. We demonstrate the utility of this tool on several complex specimen types, including urine, wound swabs and tissue, and several types of respiratory and fecal specimens, showing K. pneumoniae species and clonal group identification and antimicrobial resistance and virulence profiling, including capsule typing. Use of this amplicon sequencing tool to screen patients for Klebsiella carriage could inform health care staff of the risk of infection and outbreak potential. KlebSeq also serves as a model for next-generation molecular tools for public health and health care, as expansion of this tool can be used for several other HAI agents or applications.

Cryptic Microheteroresistance Explains Mycobacterium tuberculosis Phenotypic Resistance

Rationale: Minority drug‐resistant Mycobacterium tuberculosis subpopulations can be associated with phenotypic resistance but are poorly detected by Sanger sequencing or commercial molecular diagnostic assays. Objectives: To determine the role of targeted next‐generation sequencing in resolving these minor variant subpopulations. Methods: We used single molecule overlapping reads (SMOR), a targeted next‐generation sequencing approach that dramatically reduces sequencing error, to analyze primary cultured isolates phenotypically resistant to rifampin, fluoroquinolones, or aminoglycosides, but for which Sanger sequencing found no resistance‐associated variants (RAVs) within respective resistance‐determining regions (study group). Isolates also underwent single‐colony selection on antibiotic‐containing agar, blinded to sequencing results. As a positive control, isolates with multiple colocalizing chromatogram peaks were also analyzed (control group). Measurements and Main Results: Among 61 primary culture isolates (25 study group and 36 control group), SMOR described 66 (49%) and 45 (33%) of 135 total heteroresistant RAVs at frequencies less than 5% and less than 1% of the total mycobacterial population, respectively. In the study group, SMOR detected minor resistant variant subpopulations in 80% (n = 20/25) of isolates with no Sanger‐identified RAVs (median subpopulation size, 1.0%; interquartile range, 0.2‐3.9%). Single‐colony selection on drug‐containing media corroborated SMOR results for 90% (n = 18/20) of RAV‐containing specimens, and the absence of RAVs in 60% (n = 3/5) of isolates. Overall, Sanger sequencing was concordant with SMOR for 77% (n = 53/69) of macroheteroresistant (5‐95% total population), but only 5% of microheteroresistant (<5%) subpopulations (n = 3/66) across both groups. Conclusions: Cryptic minor variant mycobacterial subpopulations exist below the resolving capability of current drug susceptibility testing methodologies, and may explain an important proportion of false‐negative resistance determinations.

Whole-Genome Sequencing to Determine Origin of Multinational Outbreak of Sarocladium kiliense Bloodstream Infections.

Next-generation technologies and bioinformatics enabled source attribution and implementation of effective control strategies.

Improved Subtyping of Staphylococcus aureus Clonal Complex 8 Strains Based on Whole-Genome Phylogenetic Analysis

Staphylococcus aureus is a major human pathogen worldwide in both community and health care settings. Surveillance for S. aureus strains is important to our understanding of their spread and to informing infection prevention and control. Confusion surrounding the strain nomenclature of one of the most prevalent lineages of S. aureus, clonal complex 8 (CC8), and the imprecision of current tools for typing S. aureus make surveillance and source tracing difficult and sometimes misleading. In this study, we clarify the CC8 strain designations and propose a new typing scheme for CC8 isolates that is rapid and easy to use. This typing scheme is based on relatively stable genomic markers, and we demonstrate its superiority over traditional typing techniques. This scheme has the potential to greatly improve epidemiological investigations of S. aureus. ABSTRACT Strains of Staphylococcus aureus in clonal complex 8 (CC8), including USA300, USA500, and the Iberian clone, are prevalent pathogens in the United States, both inside and outside health care settings. Methods for typing CC8 strains are becoming obsolete as the strains evolve and diversify, and whole-genome sequencing has shown that some strain types fall into multiple sublineages within CC8. In this study, we attempt to clarify the strain nomenclature of CC8, classifying the major strain types based on whole-genome sequence phylogenetics using both methicillin-resistant S. aureus (MRSA) and methicillin-susceptible S. aureus (MSSA) genomes. We show that isolates of the Archaic and Iberian clones from decades ago make up the most basal clade of the main CC8 lineages and that at least one successful lineage of CC8, made up mostly of MSSA, diverged before the other well-known strain types USA500 and USA300. We also show that the USA500 type includes two clades separated by the previously described “Canadian epidemic MRSA” strain CMRSA9, that one clade containing USA500 also contains the USA300 clade, and that the USA300-0114 strain type is not a monophyletic group. Additionally, we present a rapid, simple CC8 strain-typing scheme using real-time PCR assays that target single nucleotide polymorphisms (SNPs) derived from our CC8 phylogeny and show the significant benefit of using more stable genomic markers based on evolutionary lineages over traditional S. aureus typing techniques. This more accurate and accessible S. aureus typing system may improve surveillance and better inform the epidemiology of this very important pathogen. IMPORTANCE Staphylococcus aureus is a major human pathogen worldwide in both community and health care settings. Surveillance for S. aureus strains is important to our understanding of their spread and to informing infection prevention and control. Confusion surrounding the strain nomenclature of one of the most prevalent lineages of S. aureus, clonal complex 8 (CC8), and the imprecision of current tools for typing S. aureus make surveillance and source tracing difficult and sometimes misleading. In this study, we clarify the CC8 strain designations and propose a new typing scheme for CC8 isolates that is rapid and easy to use. This typing scheme is based on relatively stable genomic markers, and we demonstrate its superiority over traditional typing techniques. This scheme has the potential to greatly improve epidemiological investigations of S. aureus.

Direct detection of Coccidioides from Arizona soils using CocciENV, a highly sensitive and specific real-time PCR assay

Coccidioides immitis and Coccidioides posadasii are soil fungi endemic to desert regions of the southwestern United States, and the causative agents of valley fever, or coccidioidomycosis. Studies have shown that the distribution of Coccidioides in soils is sporadic and cannot be explained by soil characteristics alone, suggesting that biotic and other abiotic factors should be examined. However, tools to reliably and robustly screen the large number of soils needed to investigate these potential associations have not been available. Thus, we developed a real-time polymerase chain reaction (PCR) assay for testing environmental samples by modifying CocciDx, an assay validated for testing clinical specimens to facilitate coccidioidomycosis diagnosis. For this study, we collected soil samples from previously established locations of C. posadasii in Arizona and new locations in fall 2013 and spring 2014, and screened the extracted DNA with the new assay known as CocciEnv. To verify the presence of Coccidioides in soil using an alternate method, we employed next generation amplicon sequencing targeting the ITS2 region. Results show our modified assay, CocciEnv, is a rapid and robust method for detecting Coccidioides DNA in complex environmental samples. The ability to test a large number of soils for the presence of Coccidioides is a much-needed tool in the understanding of the ecology of the organism and epidemiology of the disease and will greatly improve our understanding of this human pathogen.