Genomic analysis of carbapenemase-encoding plasmids from Klebsiella pneumoniae across Europe highlights three major patterns of dissemination

AbstractThe incidence of Klebsiella pneumoniae infections that are resistant to carbapenems, a last-line class of antibiotics, has been rapidly increasing. The primary mechanism of carbapenem resistance is production of carbapenemase enzymes, which are most frequently encoded on plasmids by blaOXA-48-like, blaVIM, blaNDM and blaKPC genes. Using short-read sequence data, we previously analysed genomes of 1717 isolates from the K. pneumoniae species complex submitted during the European survey of carbapenemase-producing Enterobacteriaceae (EuSCAPE). Here, we investigated the diversity, prevalence and transmission dynamics of carbapenemase-encoding plasmids using long-read sequencing of representative isolates (n=79) from this collection in combination with short-read data from all isolates. We highlight three major patterns by which carbapenemase genes have disseminated via plasmids. First, blaOXA-48-like genes have spread across diverse lineages primarily via a highly conserved, epidemic pOXA-48-like plasmid. Second, blaVIM and blaNDM genes have spread via transient associations of diverse plasmids with numerous lineages. Third, blaKPC genes have transmitted predominantly by stable association with one clonal lineage (ST258/512) despite frequent mobilisation between pre-existing yet diverse plasmids within the lineage. Despite contrasts in these three modes of carbapenemase gene spread, which can be summarised as using one plasmid/multiple lineages, multiple plasmids/multiple lineages, and multiple plasmids/one lineage, all are underpinned by significant propagation along high-risk clonal lineages.

Download Full-text

Integrated chromosomal and plasmid sequence analyses reveal diverse modes of carbapenemase gene spread among Klebsiella pneumoniae

Proceedings of the National Academy of Sciences ◽

10.1073/pnas.2003407117 ◽

2020 ◽

Vol 117 (40) ◽

pp. 25043-25054 ◽

Cited By ~ 2

Author(s):

Sophia David ◽

Victoria Cohen ◽

Sandra Reuter ◽

Anna E. Sheppard ◽

Tommaso Giani ◽

...

Keyword(s):

Antibiotic Resistance ◽

Klebsiella Pneumoniae ◽

Resistance Genes ◽

Sequence Data ◽

Antibiotic Resistance Genes ◽

Clinical Settings ◽

Surveillance Systems ◽

Carbapenemase Gene ◽

Short Read Sequence ◽

Stable Association

Molecular and genomic surveillance systems for bacterial pathogens currently rely on tracking clonally evolving lineages. By contrast, plasmids are usually excluded or analyzed with low-resolution techniques, despite being the primary vectors of antibiotic resistance genes across many key pathogens. Here, we used a combination of long- and short-read sequence data of Klebsiella pneumoniae isolates (n = 1,717) from a European survey to perform an integrated, continent-wide study of chromosomal and plasmid diversity. This revealed three contrasting modes of dissemination used by carbapenemase genes, which confer resistance to last-line carbapenems. First, blaOXA-48-like genes have spread primarily via the single epidemic pOXA-48–like plasmid, which emerged recently in clinical settings and spread rapidly to numerous lineages. Second, blaVIM and blaNDM genes have spread via transient associations of many diverse plasmids with numerous lineages. Third, blaKPC genes have transmitted predominantly by stable association with one successful clonal lineage (ST258/512) yet have been mobilized among diverse plasmids within this lineage. We show that these plasmids, which include pKpQIL-like and IncX3 plasmids, have a long association (and are coevolving) with the lineage, although frequent recombination and rearrangement events between them have led to a complex array of mosaic plasmids carrying blaKPC. Taken altogether, these results reveal the diverse trajectories of antibiotic resistance genes in clinical settings, summarized as using one plasmid/multiple lineages, multiple plasmids/multiple lineages, and multiple plasmids/one lineage. Our study provides a framework for the much needed incorporation of plasmid data into genomic surveillance systems, an essential step toward a more comprehensive understanding of resistance spread.

Download Full-text

Using Short Read Sequencing to Characterise Balanced Reciprocal Translocations in Pigs

10.21203/rs.3.rs-28830/v1 ◽

2020 ◽

Author(s):

Aniek Cornelia Bouwman ◽

Martijn F.L. Derks ◽

Marleen L.W.J. Broekhuijse ◽

Barbara Harlizius ◽

Roel F. Veerkamp

Keyword(s):

Sequence Data ◽

Variant Calling ◽

Reciprocal Translocations ◽

Short Read ◽

Short Read Sequencing ◽

Long Read ◽

Short Read Sequence ◽

Staining Techniques ◽

Chromosome Staining ◽

Paired End Sequencing

Abstract Background A balanced constitutional reciprocal translocation (RT) is a mutual exchange of terminal segments of two non-homologous chromosomes without any loss or gain of DNA in germline cells. Carriers of balanced RTs are viable individuals with no apparent phenotypical consequences. These animals produce, however, unbalanced gametes and show therefore reduced fertility and offspring with congenital abnormalities. This cytogenetic abnormality is usually detected using chromosome staining techniques. The aim of this study was to test the possibilities of using paired end short read sequencing for detection of balanced RTs in boars and investigate their breakpoints and junctions.Results Balanced RTs were recovered in a blinded analysis, using structural variant calling software DELLY, in 6 of the 7 carriers with 30 fold short read paired end sequencing. In 15 non-carriers we did not detect any RTs. Reducing the coverage to 20 fold, 15 fold and 10 fold showed that at least 20 fold coverage is required to obtain good results. One RT was not detected using the blind screening, however, a highly likely RT was discovered after unblinding. This RT was located in a repetitive region, showing the limitations of short read sequence data. The detailed analysis of the breakpoints and junctions suggested three junctions showing microhomology, three junctions with blunt-end ligation, and three micro-insertions at the breakpoint junctions. The RTs detected also showed to disrupt genes.Conclusions We conclude that paired end short read sequence data can be used to detect and characterize balanced reciprocal translocations, if sequencing depth is at least 20 fold coverage. However, translocations in repetitive areas may require large fragments or even long read sequence data.

Download Full-text

Paragraph: a graph-based structural variant genotyper for short-read sequence data

Genome Biology ◽

10.1186/s13059-019-1909-7 ◽

2019 ◽

Vol 20 (1) ◽

Cited By ~ 19

Author(s):

Sai Chen ◽

Peter Krusche ◽

Egor Dolzhenko ◽

Rachel M. Sherman ◽

Roman Petrovski ◽

...

Keyword(s):

Sequence Data ◽

Whole Genome Sequence ◽

Whole Genome ◽

Structural Variations ◽

Short Read ◽

Three Samples ◽

Genomics Research ◽

Long Read ◽

Short Read Sequence ◽

Population Scale

AbstractAccurate detection and genotyping of structural variations (SVs) from short-read data is a long-standing area of development in genomics research and clinical sequencing pipelines. We introduce Paragraph, an accurate genotyper that models SVs using sequence graphs and SV annotations. We demonstrate the accuracy of Paragraph on whole-genome sequence data from three samples using long-read SV calls as the truth set, and then apply Paragraph at scale to a cohort of 100 short-read sequenced samples of diverse ancestry. Our analysis shows that Paragraph has better accuracy than other existing genotypers and can be applied to population-scale studies.

Download Full-text

Using short read sequencing to characterise balanced reciprocal translocations in pigs

10.21203/rs.3.rs-28830/v3 ◽

2020 ◽

Author(s):

Aniek C. Bouwman ◽

Martijn F.L. Derks ◽

Marleen L.W.J. Broekhuijse ◽

Barbara Harlizius ◽

Roel F. Veerkamp

Keyword(s):

Sequence Data ◽

Variant Calling ◽

Reciprocal Translocations ◽

Short Read ◽

Short Read Sequencing ◽

Long Read ◽

Short Read Sequence ◽

Staining Techniques ◽

Chromosome Staining ◽

Paired End Sequencing

Abstract Background A balanced constitutional reciprocal translocation (RT) is a mutual exchange of terminal segments of two non-homologous chromosomes without any loss or gain of DNA in germline cells. Carriers of balanced RTs are viable individuals with no apparent phenotypical consequences. These animals produce, however, unbalanced gametes and show therefore reduced fertility and offspring with congenital abnormalities. This cytogenetic abnormality is usually detected using chromosome staining techniques. The aim of this study was to test the possibilities of using paired end short read sequencing for detection of balanced RTs in boars and investigate their breakpoints and junctions. Results Balanced RTs were recovered in a blinded analysis, using structural variant calling software DELLY, in 6 of the 7 carriers with 30 fold short read paired end sequencing. In 15 non-carriers we did not detect any RTs. Reducing the coverage to 20 fold, 15 fold and 10 fold showed that at least 20 fold coverage is required to obtain good results. One RT was not detected using the blind screening, however, a highly likely RT was discovered after unblinding. This RT was located in a repetitive region, showing the limitations of short read sequence data. The detailed analysis of the breakpoints and junctions suggested three junctions showing microhomology, three junctions with blunt-end ligation, and three micro-insertions at the breakpoint junctions. The RTs detected also showed to disrupt genes. Conclusions We conclude that paired end short read sequence data can be used to detect and characterize balanced reciprocal translocations, if sequencing depth is at least 20 fold coverage. However, translocations in repetitive areas may require large fragments or even long read sequence data.

Download Full-text

Paragraph: A graph-based structural variant genotyper for short-read sequence data

10.1101/635011 ◽

2019 ◽

Cited By ~ 5

Author(s):

Sai Chen ◽

Peter Krusche ◽

Egor Dolzhenko ◽

Rachel M. Sherman ◽

Roman Petrovski ◽

...

Keyword(s):

Sequence Data ◽

Whole Genome Sequence ◽

Whole Genome ◽

Structural Variations ◽

Short Read ◽

Three Samples ◽

Genomics Research ◽

Long Read ◽

Short Read Sequence ◽

Population Scale

AbstractAccurate detection and genotyping of structural variations (SVs) from short-read data is a long-standing area of development in genomics research and clinical sequencing pipelines. We introduce Paragraph, an accurate genotyper that models SVs using sequence graphs and SV annotations. We demonstrate the accuracy of Paragraph on whole-genome sequence data from three samples using long read SV calls as the truth set, and then apply Paragraph at scale to a cohort of 100 short-read sequenced samples of diverse ancestry. Our analysis shows that Paragraph has better accuracy than other existing genotypers and can be applied to population-scale studies.

Download Full-text

Using short read sequencing to characterise balanced reciprocal translocations in pigs

10.21203/rs.3.rs-28830/v2 ◽

2020 ◽

Author(s):

Aniek C. Bouwman ◽

Martijn F.L. Derks ◽

Marleen L.W.J. Broekhuijse ◽

Barbara Harlizius ◽

Roel F. Veerkamp

Keyword(s):

Sequence Data ◽

Variant Calling ◽

Reciprocal Translocations ◽

Short Read ◽

Short Read Sequencing ◽

Long Read ◽

Short Read Sequence ◽

Staining Techniques ◽

Chromosome Staining ◽

Paired End Sequencing

Abstract Background A balanced constitutional reciprocal translocation (RT) is a mutual exchange of terminal segments of two non-homologous chromosomes without any loss or gain of DNA in germline cells. Carriers of balanced RTs are viable individuals with no apparent phenotypical consequences. These animals produce, however, unbalanced gametes and show therefore reduced fertility and offspring with congenital abnormalities. This cytogenetic abnormality is usually detected using chromosome staining techniques. The aim of this study was to test the possibilities of using paired end short read sequencing for detection of balanced RTs in boars and investigate their breakpoints and junctions. Results Balanced RTs were recovered in a blinded analysis, using structural variant calling software DELLY, in 6 of the 7 carriers with 30 fold short read paired end sequencing. In 15 non-carriers we did not detect any RTs. Reducing the coverage to 20 fold, 15 fold and 10 fold showed that at least 20 fold coverage is required to obtain good results. One RT was not detected using the blind screening, however, a highly likely RT was discovered after unblinding. This RT was located in a repetitive region, showing the limitations of short read sequence data. The detailed analysis of the breakpoints and junctions suggested three junctions showing microhomology, three junctions with blunt-end ligation, and three micro-insertions at the breakpoint junctions. The RTs detected also showed to disrupt genes. Conclusions We conclude that paired end short read sequence data can be used to detect and characterize balanced reciprocal translocations, if sequencing depth is at least 20 fold coverage. However, translocations in repetitive areas may require large fragments or even long read sequence data.

Download Full-text

Fatal Respiratory Diphtheria Caused by ß-Lactam–Resistant Corynebacterium diphtheriae

Clinical Infectious Diseases ◽

10.1093/cid/ciaa1147 ◽

2020 ◽

Cited By ~ 1

Author(s):

Brian M Forde ◽

Andrew Henderson ◽

Elliott G Playford ◽

David Looke ◽

Belinda C Henderson ◽

...

Keyword(s):

Single Molecule ◽

Genomic Analysis ◽

Carbapenem Resistance ◽

Corynebacterium Diphtheriae ◽

Penicillin Binding Protein ◽

Long Read ◽

Amoxicillin Clavulanic Acid ◽

Resistance To Penicillin ◽

High Level

Abstract Background Diphtheria is a potentially fatal respiratory disease caused by toxigenic Corynebacterium diphtheriae. Although resistance to erythromycin has been recognized, β-lactam resistance in toxigenic diphtheria has not been described. Here, we report a case of fatal respiratory diphtheria caused by toxigenic C. diphtheriae resistant to penicillin and all other β-lactam antibiotics, and describe a novel mechanism of inducible carbapenem resistance associated with the acquisition of a mobile resistance element. Methods Long-read whole-genome sequencing was performed using Pacific Biosciences Single Molecule Real-Time sequencing to determine the genome sequence of C. diphtheriae BQ11 and the mechanism of β-lactam resistance. To investigate the phenotypic inducibility of meropenem resistance, short-read sequencing was performed using an Illumina NextSeq500 sequencer on the strain both with and without exposure to meropenem. Results BQ11 demonstrated high-level resistance to penicillin (benzylpenicillin minimum inhibitory concentration [MIC] ≥ 256 μg/ml), β-lactam/β-lactamase inhibitors and cephalosporins (amoxicillin/clavulanic acid MIC ≥ 256 μg/mL; ceftriaxone MIC ≥ 8 μg/L). Genomic analysis of BQ11 identified acquisition of a novel transposon carrying the penicillin-binding protein (PBP) Pbp2c, responsible for resistance to penicillin and cephalosporins. When strain BQ11 was exposed to meropenem, selective pressure drove amplification of the transposon in a tandem array and led to a corresponding change from a low-level to a high-level meropenem-resistant phenotype. Conclusions We have identified a novel mechanism of inducible antibiotic resistance whereby isolates that appear to be carbapenem susceptible on initial testing can develop in vivo resistance to carbapenems with repeated exposure. This phenomenon could have significant implications for the treatment of C. diphtheriae infection, and may lead to clinical failure.

Download Full-text

1348. In vitro Activity of Plazomicin, a Next-Generation Aminoglycoside, Against Carbapenemase-Producing Klebsiella pneumoniae

Open Forum Infectious Diseases ◽

10.1093/ofid/ofy210.1179 ◽

2018 ◽

Vol 5 (suppl_1) ◽

pp. S412-S413

Author(s):

Michael R Jacobs ◽

Caryn E Good ◽

Ayman M Abdelhamed ◽

Daniel D Rhoads ◽

Kristine M Hujer ◽

...

Keyword(s):

Klebsiella Pneumoniae ◽

Carbapenem Resistance ◽

Vitro Activity ◽

In Vitro Activity ◽

Next Generation ◽

Aminoglycoside Resistance ◽

Carbapenem Resistant ◽

Carbapenemase Gene ◽

Research Grant

Abstract Background Plazomicin is a next-generation aminoglycoside with in vitro activity against multidrug-resistant Gram-negative species, including carbapenem-resistant isolates. The Consortium on Resistance Against Carbapenems in Klebsiella and other Enterobacteriaceae (CRACKLE) is a federally funded, prospective multicenter consortium of 20 hospitals from nine US healthcare systems to track carbapenem-resistant Enterobacteriaceae. Methods Minimum inhibitory concentrations (MICs) of plazomicin were determined by broth microdilution according to current CLSI guidelines against a collection of 697 carbapenem-resistant Klebsiella pneumoniae with defined carbapenem resistance mechanisms, including KPC and OXA carbapenemases. Isolates were submitted by participating CRACKLE centers. Results Carbapenemases present in study isolates included KPC-2 (n = 323), KPC-3 (n = 364), KPC-4 (n = 2), OXA-48 like (n = 7), and NDM (n = 1). Plazomicin MICs ranged from ≤0.12 to >32 mg/L, with MIC50 and MIC90 values of 0.25 and 1 mg/L, respectively (figure). MICs of 689 (98.8%) isolates were ≤4 mg/L, while MICs of the remaining eight isolates were >32 mg/L. Plazomicin MICs were related to specific carbapenemases present in isolates: of eight isolates with MICs >32 mg/L, seven contained OXA-48 like and one contained KPC-3, suggesting that these isolates possess an aminoglycoside-resistance mechanism on the same plasmid as their carbapenemase gene, such as a 16S ribosomal RNA methyltransferase, against which plazomicin is not active. Conclusion Plazomicin has good in vitro potency against a collection of carbapenemase-producing K. pneumoniae, with MIC90 value of 1 mg/L and MICs of ≤4 mg/L for 98.9% of isolates. Disclosures M. R. Jacobs, Achaogen: Investigator, Research grant. Shionogi: Investigator, Research grant. L. Connolly, Achaogen, Inc.: Consultant, Consulting fee. K. M. Krause, Achaogen: Employee, Salary. S. S. Richter, bioMerieux: Grant Investigator, Research grant. BD Diagnostics: Grant Investigator, Research grant. Roche: Grant Investigator, Research grant. Hologic: Grant Investigator, Research grant. Diasorin: Grant Investigator, Research grant. Accelerate: Grant Investigator, Research grant. Biofire: Grant Investigator, Research grant. D. Van Duin, achaogen: Scientific Advisor, Consulting fee. shionogi: Scientific Advisor, Consulting fee. Allergan: Scientific Advisor, Consulting fee. Astellas: Scientific Advisor, Consulting fee. Neumedicine: Scientific Advisor, Consulting fee. Roche: Scientific Advisor, Consulting fee. T2 Biosystems: Scientific Advisor, Consulting fee.

Download Full-text

Rapid Mycobacterium tuberculosis spoligotyping from uncorrected long reads using Galru

10.1101/2020.05.31.126490 ◽

2020 ◽

Author(s):

Andrew J. Page ◽

Nabil-Fareed Alikhan ◽

Michael Strinden ◽

Thanh Le Viet ◽

Timofey Skvortsov

Keyword(s):

Mycobacterium Tuberculosis ◽

State Of The Art ◽

Sequence Data ◽

Human Pathogen ◽

Sequencing Data ◽

Short Read ◽

Short Read Sequencing ◽

Long Reads ◽

Long Read

AbstractSpoligotyping of Mycobacterium tuberculosis provides a subspecies classification of this major human pathogen. Spoligotypes can be predicted from short read genome sequencing data; however, no methods exist for long read sequence data such as from Nanopore or PacBio. We present a novel software package Galru, which can rapidly detect the spoligotype of a Mycobacterium tuberculosis sample from as little as a single uncorrected long read. It allows for near real-time spoligotyping from long read data as it is being sequenced, giving rapid sample typing. We compare it to the existing state of the art software and find it performs identically to the results obtained from short read sequencing data. Galru is freely available from https://github.com/quadram-institute-bioscience/galru under the GPLv3 open source licence.

Download Full-text

Complete Genome Resequencing of Thermus thermophilus Strain TMY by Hybrid Assembly of Long- and Short-Read Sequencing Technologies

Microbiology Resource Announcements ◽

10.1128/mra.00979-21 ◽

2021 ◽

Vol 10 (46) ◽

Author(s):

Kentaro Miyazaki ◽

Natsuko Tokito

Keyword(s):

Complete Genome ◽

Thermus Thermophilus ◽

Genomic Analysis ◽

Comparative Genomic ◽

Hybrid Assembly ◽

Genome Resequencing ◽

Short Read ◽

Content Type ◽

Sequencing Technologies ◽

Long Read

Complete genome resequencing was conducted for Thermus thermophilus strain TMY by hybrid assembly of Oxford Nanopore Technologies long-read and MGI short-read data. Errors in the previously reported genome sequence determined by PacBio technology alone were corrected, allowing for high-quality comparative genomic analysis of closely related T. thermophilus genomes.

Download Full-text