Plassembler: an automated bacterial plasmid assembly tool.

SUMMARY: With recent advances in sequencing technologies, it is now possible to obtain near-perfect complete bacterial chromosome assemblies cheaply and efficiently by combining a long-read-first assembly approach with short-read polishing. However, existing methods for assembling bacterial plasmids from long-read-first assemblies often misassemble or even miss bacterial plasmids entirely and accordingly require manual curation. Plassembler was developed to provide a tool that automatically assembles and outputs bacterial plasmids using a hybrid assembly approach. It achieves increased accuracy and computational efficiency compared to the existing gold standard tool Unicycler by removing chromosomal reads from the input read sets using a mapping approach. AVAILABILITY AND IMPLEMENTATION: Plassembler is implemented in Python and is installable as a bioconda package using 'conda install -c bioconda plassembler'. The source code is available on GitHub at https://github.com/gbouras13/plassembler. The full benchmarking pipeline can be found at https://github.com/gbouras13/plassembler_simulation_benchmarking, while the benchmarking input FASTQ and output files can be found at https://doi.org/10.5281/zenodo.7996690.

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

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.

A Candida auris Outbreak and Its Control in an Intensive Care Setting.

BACKGROUND: Candida auris is an emerging and multidrug-resistant pathogen. Here we report the epidemiology of a hospital outbreak of C. auris colonization and infection. METHODS: After identification of a cluster of C. auris infections in the neurosciences intensive care unit (ICU) of the Oxford University Hospitals, United Kingdom, we instituted an intensive patient and environmental screening program and package of interventions. Multivariable logistic regression was used to identify predictors of C. auris colonization and infection. Isolates from patients and from the environment were analyzed by whole-genome sequencing. RESULTS: A total of 70 patients were identified as being colonized or infected with C. auris between February 2, 2015, and August 31, 2017; of these patients, 66 (94%) had been admitted to the neurosciences ICU before diagnosis. Invasive C. auris infections developed in 7 patients. When length of stay in the neurosciences ICU and patient vital signs and laboratory results were controlled for, the predictors of C. auris colonization or infection included the use of reusable skin-surface axillary temperature probes (multivariable odds ratio, 6.80; 95% confidence interval [CI], 2.96 to 15.63; P<0.001) and systemic fluconazole exposure (multivariable odds ratio, 10.34; 95% CI, 1.64 to 65.18; P=0.01). C. auris was rarely detected in the general environment. However, it was detected in isolates from reusable equipment, including multiple axillary skin-surface temperature probes. Despite a bundle of infection-control interventions, the incidence of new cases was reduced only after removal of the temperature probes. All outbreak sequences formed a single genetic cluster within the C. auris South African clade. The sequenced isolates from reusable equipment were genetically related to isolates from the patients. CONCLUSIONS: The transmission of C. auris in this hospital outbreak was found to be linked to reusable axillary temperature probes, indicating that this emerging pathogen can persist in the environment and be transmitted in health care settings. (Funded by the National Institute for Health Research Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance at Oxford University and others.).

Reconciling the Potentially Irreconcilable? Genotypic and Phenotypic Amoxicillin-Clavulanate Resistance in Escherichia coli.

Resistance to amoxicillin-clavulanate, a widely used beta-lactam/beta-lactamase inhibitor combination antibiotic, is rising globally, and yet susceptibility testing remains challenging. To test whether whole-genome sequencing (WGS) could provide a more reliable assessment of susceptibility than traditional methods, we predicted resistance from WGS for 976 Escherichia coli bloodstream infection isolates from Oxfordshire, United Kingdom, comparing against phenotypes from the BD Phoenix (calibrated against EUCAST guidelines). A total of 339/976 (35%) isolates were amoxicillin-clavulanate resistant. Predictions based solely on beta-lactamase presence/absence performed poorly (sensitivity, 23% [78/339]) but improved when genetic features associated with penicillinase hyperproduction (e.g., promoter mutations and copy number estimates) were considered (sensitivity, 82% [277/339]; P < 0.0001). Most discrepancies occurred in isolates with MICs within ±1 doubling dilution of the breakpoint. We investigated two potential causes: the phenotypic reference and the binary resistant/susceptible classification. We performed reference standard, replicated phenotyping in a random stratified subsample of 261/976 (27%) isolates using agar dilution, following both EUCAST and CLSI guidelines, which use different clavulanate concentrations. As well as disagreeing with each other, neither agar dilution phenotype aligned perfectly with genetic features. A random-effects model investigating associations between genetic features and MICs showed that some genetic features had small, variable and additive effects, resulting in variable resistance classification. Using model fixed-effects to predict MICs for the non-agar dilution isolates, predicted MICs were in essential agreement (±1 doubling dilution) with observed (BD Phoenix) MICs for 691/715 (97%) isolates. This suggests amoxicillin-clavulanate resistance in E. coli is quantitative, rather than qualitative, explaining the poorly reproducible binary (resistant/susceptible) phenotypes and suboptimal concordance between different phenotypic methods and with WGS-based predictions.

Klebsiella quasipneumoniae Provides a Window into Carbapenemase Gene Transfer, Plasmid Rearrangements, and Patient Interactions with the Hospital Environment.

Several emerging pathogens have arisen as a result of selection pressures exerted by modern health care. Klebsiella quasipneumoniae was recently defined as a new species, yet its prevalence, niche, and propensity to acquire antimicrobial resistance genes are not fully described. We have been tracking inter- and intraspecies transmission of the Klebsiella pneumoniae carbapenemase (KPC) gene, blaKPC, between bacteria isolated from a single institution. We applied a combination of Illumina and PacBio whole-genome sequencing to identify and compare K. quasipneumoniae from patients and the hospital environment over 10- and 5-year periods, respectively. There were 32 blaKPC-positive K. quasipneumoniae isolates, all of which were identified as K. pneumoniae in the clinical microbiology laboratory, from 8 patients and 11 sink drains, with evidence for seven separate blaKPC plasmid acquisitions. Analysis of a single subclade of K. quasipneumoniae subsp. quasipneumoniae (n = 23 isolates) from three patients and six rooms demonstrated seeding of a sink by a patient, subsequent persistence of the strain in the hospital environment, and then possible transmission to another patient. Longitudinal analysis of this strain demonstrated the acquisition of two unique blaKPC plasmids and then subsequent within-strain genetic rearrangement through transposition and homologous recombination. Our analysis highlights the apparent molecular propensity of K. quasipneumoniae to persist in the environment as well as acquire carbapenemase plasmids from other species and enabled an assessment of the genetic rearrangements which may facilitate horizontal transmission of carbapenemases.

The Role of fosA in Challenges with Fosfomycin Susceptibility Testing of Multispecies Klebsiella pneumoniae Carbapenemase-Producing Clinical Isolates.

With multidrug-resistant (MDR) Enterobacterales on the rise, a nontoxic antimicrobial agent with a unique mechanism of action such as fosfomycin seems attractive. However, establishing accurate fosfomycin susceptibility testing for non-Escherichia coli isolates in a clinical microbiology laboratory remains problematic. We evaluated fosfomycin susceptibility by multiple methods with 96 KPC-producing clinical isolates of multiple strains and species collected at a single center between 2008 and 2016. In addition, we assessed the presence of fosfomycin resistance genes from whole-genome sequencing (WGS) data using NCBI's AMRFinder and custom HMM search. Susceptibility testing was performed using a glucose-6-phosphate-supplemented fosfomycin Etest and Kirby-Bauer disk diffusion (DD) assays, and the results were compared to those obtained by agar dilution. Clinical Laboratory and Standards Institute (CLSI) breakpoints for E. coli were applied for interpretation. Overall, 63% (60/96) of isolates were susceptible by Etest, 70% (67/96) by DD, and 88% (84/96) by agar dilution. fosA was detected in 80% (70/88) of previously sequenced isolates, with species-specific associations and alleles, and fosA-positive isolates were associated with higher MIC distributions. Disk potentiation testing was performed using sodium phosphonoformate to inhibit fosA and showed significant increases in the zone diameter of DD testing for isolates that were fosA positive compared to those that were fosA negative. The addition of sodium phosphonoformate (PPF) corrected 10/14 (71%) major errors in categorical agreement with agar dilution. Our results indicate that fosA influences the inaccuracy of susceptibility testing by methods readily available in a clinical laboratory compared to agar dilution. Further research is needed to determine the impact of fosA on clinical outcomes.

Don't overlook the little guy: An evaluation of the frequency of small plasmids co-conjugating with larger carbapenemase gene containing plasmids.

As the spread of antimicrobial resistance (AMR) genes becomes an increasing global threat, improved understanding of mobile genetic elements which contribute to the spread of antimicrobial resistance genes, becomes more critical. We created transconjugants from the mating of three chromosomally isogenic Klebsiella pneumoniae carbapenemase (blaKPC) positive Citrobacter freundii isolates with a laboratory strain of Escherichia coli and evaluated the movement of small cryptic plasmids (SCPs), p3223 and p1916, when larger blaKPC-plasmids were transferred. In all of the 143 transconjugants, multiple plasmids, both large and small, transferred with each mating. When two blaKPC-plasmids were present in the host, frequently (87%; 98/113) both would be transferred during mating. p3223 is found in a wide range of bacterial hosts that harbor AMR genes; p1916 has been identified in only a limited number of publicly available sequences to date. From our evaluation, there is still much to learn about SCPs, and the high rate of co-transfer of multiple plasmids from real-world carbapenemase-producing Enterobacteriales.

High Rates of Human Fecal Carriage of mcr-1-Positive Multidrug-Resistant Enterobacteriaceae Emerge in China in Association With Successful Plasmid Families.

Background: mcr-1-mediated colistin resistance in Enterobacteriaceae is concerning, as colistin is used in treating multidrug-resistant Enterobacteriaceae infections. We identified trends in human fecal mcr-1-positivity rates and colonization with mcr-1-positive, third-generation cephalosporin-resistant (3GC-R) Enterobacteriaceae in Guangzhou, China, and investigated the genetic contexts of mcr-1 in mcr-1-positive 3GC-R strains. Methods: Fecal samples were collected from in-/out-patients submitting specimens to 3 hospitals (2011-2016). mcr-1 carriage trends were assessed using iterative sequential regression. A subset of mcr-1-positive isolates was sequenced (whole-genome sequencing [WGS], Illumina), and genetic contexts (flanking regions, plasmids) of mcr-1 were characterized. Results: Of 8022 fecal samples collected, 497 (6.2%) were mcr-1 positive, and 182 (2.3%) harbored mcr-1-positive 3GC-R Enterobacteriaceae. We observed marked increases in mcr-1 (0% [April 2011] to 31% [March 2016]) and more recent (since January 2014; 0% [April 2011] to 15% [March 2016]) increases in human colonization with mcr-1-positive 3GC-R Enterobacteriaceae (P < .001). mcr-1-positive 3GC-R isolates were commonly multidrug resistant. WGS of mcr-1-positive 3GC-R isolates (70 Escherichia coli, 3 Klebsiella pneumoniae) demonstrated bacterial strain diversity; mcr-1 in association with common plasmid backbones (IncI, IncHI2/HI2A, IncX4) and sometimes in multiple plasmids; frequent mcr-1 chromosomal integration; and high mobility of the mcr-1-associated insertion sequence ISApl1. Sequence data were consistent with plasmid spread among animal/human reservoirs. Conclusions: The high prevalence of mcr-1 in multidrug-resistant E. coli colonizing humans is a clinical threat; diverse genetic mechanisms (strains/plasmids/insertion sequences) have contributed to the dissemination of mcr-1, and will facilitate its persistence.

Host-Pathogen Coevolution: The Selective Advantage of Bacillus thuringiensis Virulence and Its Cry Toxin Genes.

Reciprocal coevolution between host and pathogen is widely seen as a major driver of evolution and biological innovation. Yet, to date, the underlying genetic mechanisms and associated trait functions that are unique to rapid coevolutionary change are generally unknown. We here combined experimental evolution of the bacterial biocontrol agent Bacillus thuringiensis and its nematode host Caenorhabditis elegans with large-scale phenotyping, whole genome analysis, and functional genetics to demonstrate the selective benefit of pathogen virulence and the underlying toxin genes during the adaptation process. We show that: (i) high virulence was specifically favoured during pathogen-host coevolution rather than pathogen one-sided adaptation to a nonchanging host or to an environment without host; (ii) the pathogen genotype BT-679 with known nematocidal toxin genes and high virulence specifically swept to fixation in all of the independent replicate populations under coevolution but only some under one-sided adaptation; (iii) high virulence in the BT-679-dominated populations correlated with elevated copy numbers of the plasmid containing the nematocidal toxin genes; (iv) loss of virulence in a toxin-plasmid lacking BT-679 isolate was reconstituted by genetic reintroduction or external addition of the toxins. We conclude that sustained coevolution is distinct from unidirectional selection in shaping the pathogen's genome and life history characteristics. To our knowledge, this study is the first to characterize the pathogen genes involved in coevolutionary adaptation in an animal host-pathogen interaction system.

Chromosomal Integration of the Klebsiella pneumoniae Carbapenemase Gene, blaKPC, in Klebsiella Species Is Elusive but Not Rare.

Carbapenemase genes in Enterobacteriaceae are mostly described as being plasmid associated. However, the genetic context of carbapenemase genes is not always confirmed in epidemiological surveys, and the frequency of their chromosomal integration therefore is unknown. A previously sequenced collection of blaKPC-positive Enterobacteriaceae from a single U.S. institution (2007 to 2012; n = 281 isolates from 182 patients) was analyzed to identify chromosomal insertions of Tn4401, the transposon most frequently harboring blaKPC Using a combination of short- and long-read sequencing, we confirmed five independent chromosomal integration events from 6/182 (3%) patients, corresponding to 15/281 (5%) isolates. Three patients had isolates identified by perirectal screening, and three had infections which were all successfully treated. When a single copy of blaKPC was in the chromosome, one or both of the phenotypic carbapenemase tests were negative. All chromosomally integrated blaKPC genes were from Klebsiella spp., predominantly K. pneumoniae clonal group 258 (CG258), even though these represented only a small proportion of the isolates. Integration occurred via IS15-ΔI-mediated transposition of a larger, composite region encompassing Tn4401 at one locus of chromosomal integration, seen in the same strain (K. pneumoniae ST340) in two patients. In summary, we identified five independent chromosomal integrations of blaKPC in a large outbreak, demonstrating that this is not a rare event. blaKPC was more frequently integrated into the chromosome of epidemic CG258 K. pneumoniae lineages (ST11, ST258, and ST340) and was more difficult to detect by routine phenotypic methods in this context. The presence of chromosomally integrated blaKPC within successful, globally disseminated K. pneumoniae strains therefore is likely underestimated.

Enhanced Klebsiella pneumoniae Carbapenemase Expression from a Novel Tn4401 Deletion.

The Klebsiella pneumoniae carbapenemase gene (blaKPC) is typically located within mobile transposon Tn4401 Enhanced KPC expression has been associated with deletions in the putative promoter region upstream of blaKPC Illumina sequences from blaKPC-positive clinical isolates from a single institution were mapped to a Tn4401b reference sequence, which carries no deletions. The novel isoform Tn4401h (188-bp deletion [between istB and blaKPC]) was present in 14% (39/281) of clinical isolates. MICs showed that Escherichia coli strains containing plasmids with Tn4401a and Tn4401h were more resistant to meropenem (≥16 and ≥16, respectively), ertapenem (≥8 and 4, respectively), and cefepime (≥64 and 4, respectively) than E. coli strains with Tn4401b (0.5, ≤0.5, and ≤1, respectively). Quantitative real-time PCR (qRT-PCR) demonstrated that Tn4401a had a 16-fold increase and Tn4401h a 4-fold increase in blaKPC mRNA levels compared to the reference Tn4401b. A lacZ reporter plasmid was used to test the activity of the promoter regions from the different variants, and the results showed that the Tn4401a and Tn4401h promoter sequences generated higher ß-galactosidase activity than the corresponding Tn4401b sequence. Further dissection of the promoter region demonstrated that putative promoter P1 was not functional. The activity of the isolated P2 promoter was greatly enhanced by inclusion of the P1-P2 intervening sequence. These studies indicated that gene expression could be an important consideration in understanding resistance phenotypes predicted by genetic signatures in the context of sequencing-based rapid diagnostics.

Ordering the mob: Insights into replicon and MOB typing schemes from analysis of a curated dataset of publicly available plasmids.

Plasmid typing can provide insights into the epidemiology and transmission of plasmid-mediated antibiotic resistance. The principal plasmid typing schemes are replicon typing and MOB typing, which utilize variation in replication loci and relaxase proteins respectively. Previous studies investigating the proportion of plasmids assigned a type by these schemes ('typeability') have yielded conflicting results; moreover, thousands of plasmid sequences have been added to NCBI in recent years, without consistent annotation to indicate which sequences represent complete plasmids. Here, a curated dataset of complete Enterobacteriaceae plasmids from NCBI was compiled, and used to assess the typeability and concordance of in silico replicon and MOB typing schemes. Concordance was assessed at hierarchical replicon type resolutions, from replicon family-level to plasmid multilocus sequence type (pMLST)-level, where available. We found that 85% and 65% of the curated plasmids could be replicon and MOB typed, respectively. Overall, plasmid size and the number of resistance genes were significant independent predictors of replicon and MOB typing success. We found some degree of non-concordance between replicon families and MOB types, which was only partly resolved when partitioning plasmids into finer-resolution groups (replicon and pMLST types). In some cases, non-concordance was attributed to ambiguous boundaries between MOBP and MOBQ types; in other cases, backbone mosaicism was considered a more plausible explanation. ß-lactamase resistance genes tended not to show fidelity to a particular plasmid type, though some previously reported associations were supported. Overall, replicon and MOB typing schemes are likely to continue playing an important role in plasmid analysis, but their performance is constrained by the diverse and dynamic nature of plasmid genomes.

Complete Sequencing of Plasmids Containing blaOXA-163 and blaOXA-48 in Escherichia coli Sequence Type 131.

OXA-48-like enzymes have emerged as important extended-spectrum ß-lactamases/carbapenemases in Escherichia coli sequence type 131 (ST131). We report the structures of the first fully sequenced blaOXA-163 plasmid and of two other blaOXA-48 plasmids in this lineage. blaOXA-163 was located on a 71-kb IncN plasmid with other resistance genes. blaOXA-48 was present on IncL/M plasmids, genetically similar to other blaOXA-48 plasmid sequences, and consistent with interspecies/interlineage spread. The presence of blaOXA-48-like genes on epidemic plasmids in ST131 is of concern.

First Report of blaIMP-14 on a Plasmid Harboring Multiple Drug Resistance Genes in Escherichia coli Sequence Type 131.

The blaIMP-14 carbapenem resistance gene has largely previously been observed in Pseudomonas aeruginosa and Acinetobacter spp. As part of global surveillance and sequencing of carbapenem-resistant Escherichia coli, we identified a sequence type 131 strain harboring blaIMP-14 within a class 1 integron, itself nested within an ∼54-kb multidrug resistance region on an epidemic IncA/C2 plasmid. The emergence of blaIMP-14 in this context in the ST131 lineage is of potential clinical concern.

Nested Russian Doll-Like Genetic Mobility Drives Rapid Dissemination of the Carbapenem Resistance Gene blaKPC.

The recent widespread emergence of carbapenem resistance in Enterobacteriaceae is a major public health concern, as carbapenems are a therapy of last resort against this family of common bacterial pathogens. Resistance genes can mobilize via various mechanisms, including conjugation and transposition; however, the importance of this mobility in short-term evolution, such as within nosocomial outbreaks, is unknown. Using a combination of short- and long-read whole-genome sequencing of 281 blaKPC-positive Enterobacteriaceae isolates from a single hospital over 5 years, we demonstrate rapid dissemination of this carbapenem resistance gene to multiple species, strains, and plasmids. Mobility of blaKPC occurs at multiple nested genetic levels, with transmission of blaKPC strains between individuals, frequent transfer of blaKPC plasmids between strains/species, and frequent transposition of blaKPC transposon Tn4401 between plasmids. We also identify a common insertion site for Tn4401 within various Tn2-like elements, suggesting that homologous recombination between Tn2-like elements has enhanced the spread of Tn4401 between different plasmid vectors. Furthermore, while short-read sequencing has known limitations for plasmid assembly, various studies have attempted to overcome this by the use of reference-based methods. We also demonstrate that, as a consequence of the genetic mobility observed in this study, plasmid structures can be extremely dynamic, and therefore these reference-based methods, as well as traditional partial typing methods, can produce very misleading conclusions. Overall, our findings demonstrate that nonclonal resistance gene dissemination can be extremely rapid, presenting significant challenges for public health surveillance and achieving effective control of antibiotic resistance.

Capsular Typing Method for Streptococcus agalactiae Using Whole-Genome Sequence Data.

Group B streptococcus (GBS) capsular serotypes are major determinants of virulence and affect potential vaccine coverage. Here we report a whole-genome-sequencing-based method for GBS serotype assignment. This method shows strong agreement (kappa of 0.92) with conventional methods and increased serotype assignment (100%) to all 10 capsular types.

High instability of a nematicidal Cry toxin plasmid in Bacillus thuringiensis.

In bacterial pathogens, virulence factors are often carried on plasmids and other mobile genetic elements, and as such, plasmid evolution is central in understanding pathogenicity. Bacillus thuringiensis is an invertebrate pathogen that uses plasmid-encoded crystal (Cry) toxins to establish infections inside the host. Our study aimed to quantify stability of two Cry toxin-encoding plasmids, BTI_23p and BTI_16p, under standard laboratory culturing conditions. These two plasmids are part of the genome of the B. thuringiensis strain MYBT18679, which is of particular interest because of its high pathogenicity towards nematodes. One of the plasmids, BTI_23p, was found to be highly unstable, with substantial loss occurring within a single growth cycle. Nevertheless, longer term experimental evolution in the absence of a host revealed maintenance of the plasmid at low levels in the bacterial populations. BTI_23p encodes two nematicidal Cry toxins, Cry21Aa2 and Cry14Aa1. Consistent with previous findings, loss of the plasmid abolished pathogenicity towards the nematode Caenorhabditis elegans, which could be rescued by addition of Cry21Aa2-expressing Escherichia coli. These results implicate BTI_23p as a plasmid that is required for successful infection, yet unstable when present at high frequency in the population, consistent with the role of Cry toxins as public goods.

Klebsiella pneumoniae carbapenemase (KPC)-producing K. pneumoniae at a single institution: insights into endemicity from whole-genome sequencing.

The global emergence of Klebsiella pneumoniae carbapenemase-producing K. pneumoniae (KPC-Kp) multilocus sequence type ST258 is widely recognized. Less is known about the molecular and epidemiological details of non-ST258 K. pneumoniae in the setting of an outbreak mediated by an endemic plasmid. We describe the interplay of blaKPC plasmids and K. pneumoniae strains and their relationship to the location of acquisition in a U.S. health care institution. Whole-genome sequencing (WGS) analysis was applied to KPC-Kp clinical isolates collected from a single institution over 5 years following the introduction of blaKPC in August 2007, as well as two plasmid transformants. KPC-Kp from 37 patients yielded 16 distinct sequence types (STs). Two novel conjugative blaKPC plasmids (pKPC_UVA01 and pKPC_UVA02), carried by the hospital index case, accounted for the presence of blaKPC in 21/37 (57%) subsequent cases. Thirteen (35%) isolates represented an emergent lineage, ST941, which contained pKPC_UVA01 in 5/13 (38%) and pKPC_UVA02 in 6/13 (46%) cases. Seven (19%) isolates were the epidemic KPC-Kp strain, ST258, mostly imported from elsewhere and not carrying pKPC_UVA01 or pKPC_UVA02. Using WGS-based analysis of clinical isolates and plasmid transformants, we demonstrate the unexpected dispersal of blaKPC to many non-ST258 lineages in a hospital through spread of at least two novel blaKPC plasmids. In contrast, ST258 KPC-Kp was imported into the institution on numerous occasions, with other blaKPC plasmid vectors and without sustained transmission. Instead, a newly recognized KPC-Kp strain, ST941, became associated with both novel blaKPC plasmids and spread locally, making it a future candidate for clinical persistence and dissemination.

Hybracter: enabling scalable, automated, complete and accurate bacterial genome assemblies.

Improvements in the accuracy and availability of long-read sequencing mean that complete bacterial genomes are now routinely reconstructed using hybrid (i.e. short- and long-reads) assembly approaches. Complete genomes allow a deeper understanding of bacterial evolution and genomic variation beyond single nucleotide variants. They are also crucial for identifying plasmids, which often carry medically significant antimicrobial resistance genes. However, small plasmids are often missed or misassembled by long-read assembly algorithms. Here, we present Hybracter which allows for the fast, automatic and scalable recovery of near-perfect complete bacterial genomes using a long-read first assembly approach. Hybracter can be run either as a hybrid assembler or as a long-read only assembler. We compared Hybracter to existing automated hybrid and long-read only assembly tools using a diverse panel of samples of varying levels of long-read accuracy with manually curated ground truth reference genomes. We demonstrate that Hybracter as a hybrid assembler is more accurate and faster than the existing gold standard automated hybrid assembler Unicycler. We also show that Hybracter with long-reads only is the most accurate long-read only assembler and is comparable to hybrid methods in accurately recovering small plasmids.

Hybracter: Enabling Scalable, Automated, Complete and Accurate Bacterial Genome Assemblies.

Improvements in the accuracy and availability of long-read sequencing mean that complete bacterial genomes are now routinely reconstructed using hybrid (i.e. short- and long-reads) assembly approaches. Complete genomes allow a deeper understanding of bacterial evolution and genomic variation beyond single nucleotide variants (SNVs). They are also crucial for identifying plasmids, which often carry medically significant antimicrobial resistance (AMR) genes. However, small plasmids are often missed or misassembled by long-read assembly algorithms. Here, we present Hybracter which allows for the fast, automatic, and scalable recovery of near-perfect complete bacterial genomes using a long-read first assembly approach. Hybracter can be run either as a hybrid assembler or as a long-read only assembler. We compared Hybracter to existing automated hybrid and long-read only assembly tools using a diverse panel of samples of varying levels of long-read accuracy with manually curated ground truth reference genomes. We demonstrate that Hybracter as a hybrid assembler is more accurate and faster than the existing gold standard automated hybrid assembler Unicycler. We also show that Hybracter with long-reads only is the most accurate long-read only assembler and is comparable to hybrid methods in accurately recovering small plasmids.