Long read genome assemblers struggle with small plasmids

Complete genomes of 568 diverse Klebsiella pneumoniae species complex isolates from humans, animals, and marine sources in Norway from 2001 to 2020

We report 578 hybrid genome assemblies (568 complete) of Klebsiella pneumoniae species complex isolates from human, animal, and marine sources in Norway collected from 2001 to 2020, belonging to five phylogroups including K. pneumoniae (n = 492) and K. variicola (n = 69) and 364 unique sequence types.

Trends and characteristics of multidrug-resistant MRSA in Norway 2008-2020

Infections caused by multidrug-resistant (MDR) bacteria are recognized as a critical One Health concern which poses a significant threat to public health, leading to increased morbidity and mortality across both high- and low-income countries. In this study, we investigated the epidemiology and molecular mechanisms of multidrug-resistant methicillin-resistant Staphylococcus aureus (MDR-MRSA) strains identified in Norway from 2008 to 2020, in order to gain a better understanding of the evolution and dissemination of multidrug resistance in S. aureus. A total of 452 MDR-MRSA strains isolated from 429 individuals were analyzed from a dataset of 23,412 MRSA strains. Methods included epidemiological characterization, antimicrobial susceptibility testing (AST), and genetic analysis of a selection of strains using nanopore sequencing to identify antimicrobial resistance (AMR) genes and mutations, as well as their location on plasmids, SCCmec and other mobile genetic elements (MGEs). The study revealed an overall increasing trend in MDR-MRSA strains, with healthcare-associated strains being more prevalent among MDR-MRSA compared to the overall MRSA population. Significant heterogeneity in spa-types and clonal complexes exhibiting multidrug resistance was observed, with high resistance rates against multiple antibiotic groups, particularly erythromycin, ciprofloxacin/norfloxacin, tetracycline, gentamicin, and clindamycin in addition to cefoxitin. The predominant MDR-MRSA clones included t1476/CC8, t127/CC1, t189/CC188, and t030, t037/CC239. Among these, MRSA t1476/CC8 showed an upward trend toward the conclusion of the study period, indicating the emergence of a MDR-MRSA clone. A broad range of AMR genes and mutations were detected, linked to a wide variety of MGEs, highlighting the complex mechanisms of resistance development and dissemination within the MRSA population. This study highlights the rising challenge posed by MDR-MRSA strains, and reveals the multifactorial nature of AMR in S. aureus, thus emphasizing the importance of continued surveillance, antibiotic stewardship and infection control measures, as well as global cooperation, in order to combat the spread of these multidrug-resistant pathogens.

High intra-laboratory reproducibility of nanopore sequencing in bacterial species underscores advances in its accuracy

Nanopore sequencing is a third-generation technology known for its portability, real-time analysis and ability to generate long reads. It has great potential for use in clinical diagnostics, but thorough validation is required to address accuracy concerns and ensure reliable and reproducible results. In this study, we automated an open-source workflow (freely available at https://gitlab.com/FLI_Bioinfo/nanobacta) for the assembly of Oxford Nanopore sequencing data and used it to investigate the reproducibility of assembly results under consistent conditions. We used a benchmark dataset of five bacterial reference strains and generated eight technical sequencing replicates of the same DNA using the Ligation and Rapid Barcoding kits together with the Flongle and MinION flow cells. We assessed reproducibility by measuring discrepancies such as substitution and insertion/deletion errors, analysing plasmid recovery results and examining genetic markers and clustering information. We compared the results of genome assemblies with and without short-read polishing. Our results show an average reproducibility accuracy of 99.999955% for nanopore-only assemblies and 99.999996% when the short reads were used for polishing. The genomic analysis results were highly reproducible for the nanopore-only assemblies without short read in the following areas: identification of genetic markers for antimicrobial resistance and virulence, classical MLST, taxonomic classification, genome completeness and contamination analysis. Interestingly, the clustering information results from the core genome SNP and core genome MLST analyses were also highly reproducible for the nanopore-only assemblies, with pairwise differences of up to two allele differences in core genome MLST and two SNPs in core genome SNP across replicates. After polishing the assemblies with short reads, the pairwise differences for cgMLST were 0 and for cgSNP were 0-1 SNP across replicates. These results highlight the advances in sequencing accuracy of nanopore data without the use of short reads.

Variability in cadmium tolerance of closely related Listeria monocytogenes isolates originating from dairy processing environments

Increased tolerance to cadmium in Listeria monocytogenes has been suggested to contribute to their persistence in natural and food production environments. This study investigated the phenotypic cadmium response of L. monocytogenes strains with efflux pump cadAC (variants 1-4) and related strains with cadA1C1. Growth of cadAC variant strains (n = 5) in 0 µM-120 µM cadmium salts (CdCl2, CdSO4) in Mueller-Hinton broth (MHB) was evaluated. Additionally, 88 L. monocytogenes strains from dairy processing facilities were exposed to 43.8 µM CdCl2 in MHB, and their lag phase duration (LPD) was measured. Strains with cadA1 through cadA3 showed similar growth trends in the presence of cadmium, while the cadA4 variant (Scott A) had the highest CdCl2 minimum inhibitory concentration (175 µM). Growth varied between the two salts, with CdSO4 significantly increasing LPD (P < 0.05) compared to CdCl2. In 43.8 µM CdCl2, cadA1 strains displayed LPDs ranging from 0.99 ± 0.14 h to 6.44 ± 0.08 h, with no clear genomic differences explaining this variability. Strains without cadA did not grow at 43.8 µM CdCl2 but exhibited low tolerance (10.9 µM CdCl2), potentially due to non-specific soft metal ATPases (626 aa; 737 aa) and soft metal resistance proteins encoded by czc genes (289 aa; 291 aa; 303 aa) within their chromosomes. These findings enhance our understanding of L. monocytogenes cadmium tolerance and underscore the need for further research to explore the genetic and physiological factors underlying these trends.

IMPORTANCE

Mobile genetic elements in Listeria monocytogenes contribute to its survival in natural and food processing environments. This study focused on how different genetic variants of the efflux pump gene cadAC and group of closely related cadA1C1 strains respond to cadmium exposure. When exposed to two cadmium salts, cadmium chloride and cadmium sulfate, we observed varying growth patterns, with a significantly longer lag phase in cadmium sulfate compared to cadmium chloride. Strains with cadA1 to cadA3 had similar growth trends, whereas a strain with the cadA4 variant had the highest minimum inhibitory concentration value. Among 88 strains from dairy processing facilities, significant phenotypic differences were observed despite core genome similarities, indicating other underlying genetic and physiological factors contribute to cadmium tolerance. Since cadmium tolerance studies in L. monocytogenes are limited, with rare phenotypic comparisons between closely related strains, our study makes an important observation and contribution to understanding of L. monocytogenes tolerance to cadmium by providing phenotypic comparisons between numerous strains within the same clonal group (<16 single nucleotide polymorphisms).

Complex exchanges among plasmids and clonal expansion of lineages shape the population structure and virulence of Borrelia burgdorferi

Background

In the United States, Borrelia burgdorferi (Bb) is the principal etiologic agent of Lyme disease. The complex structure of Bb genomes has posed challenges for genomic studies because homology among the bacterium's many plasmids, which account for ~40% of the genome by length, has made them difficult to sequence and assemble.

Results

We used long-read sequencing to generate near-complete assemblies of 62 isolates of human-derived Bb and collected public genomes with plasmid sequences. We characterized genetic diversity and population structure in the resulting set of 82 plasmid-complete Borrelia burgdorferi sensu stricto genomes. The Bb core genome is encoded by a chromosome and the conserved plasmids cp26, lp54, and lp17; the accessory genome is encoded by all other plasmids and the distal arm of the chromosome. Near-complete genomes reveal that the most granular Bb genotypes are clonal expansions of complex rearrangements among accessory genome elements. Ribosomal spacer types (RST) represent multiple collections of such genotypes, whereas OspC types are usually clonal. Structural rearrangements are non-randomly distributed throughout the genome, with cp32 plasmids undergoing dense exchanges and most linear plasmids, except lp54, sharing blocks among themselves and with the distal arm of the chromosome. OspC type A strains, known to possess greater virulence in humans, are distinguished by the presence of lp28-1 and lp56. Rearrangements among plasmids tended to preserve gene content, suggesting functional constraints among gene networks. Using k-partite graph decompositions, we identified gene sets with correlation patterns suggestive of conserved functional modules.

Conclusions

Long-read assemblies reveal that Bb population genetic structure results from clonal expansion of lineages that have undergone complex rearrangements among plasmid-encoded accessory genome elements. Genetic structure is preserved among genes even when plasmid rearrangements occur, suggesting that selection among epistatic loci maintains functional genetic networks. The analysis of near-complete genomes assembled using long-read sequencing methods advances our understanding of Bb biology and Lyme disease pathogenesis by providing the first detailed view of population variation in previously inaccessible areas of the Bb genome.

Sequencing Strategy to Ensure Accurate Plasmid Assembly

Despite the wide use of plasmids in research and clinical production, the need to verify plasmid sequences is a bottleneck that is too often underestimated in the manufacturing process. Although sequencing platforms continue to improve, the method and assembly pipeline chosen still influence the final plasmid assembly sequence. Furthermore, few dedicated tools exist for plasmid assembly, especially for de novo assembly. Here, we evaluated short-read, long-read, and hybrid (both short and long reads) de novo assembly pipelines across three replicates of a 24-plasmid library. Consistent with previous characterizations of each sequencing technology, short-read assemblies had issues resolving GC-rich regions, and long-read assemblies commonly had small insertions and deletions, especially in repetitive regions. The hybrid approach facilitated the most accurate, consistent assembly generation and identified mutations relative to the reference sequence. Although Sanger sequencing can be used to verify specific regions, some GC-rich and repetitive regions were difficult to resolve using any method, suggesting that easily sequenced genetic parts should be prioritized in the design of new genetic constructs.

Plaseval: a framework for comparing and evaluating plasmid detection tools

BACKGROUND

Plasmids play a major role in the transfer of antimicrobial resistance (AMR) genes among bacteria via horizontal gene transfer. The identification of plasmids in short-read assemblies is a challenging problem and a very active research area. Plasmid binning aims at detecting, in a draft genome assembly, groups (bins) of contigs likely to originate from the same plasmid. Several methods for plasmid binning have been developed recently, such as PlasBin-flow, HyAsP, gplas, MOB-suite, and plasmidSPAdes. This motivates the problem of evaluating the performances of plasmid binning methods, either against a given ground truth or between them.

RESULTS

We describe PlasEval, a novel method aimed at comparing the results of plasmid binning tools. PlasEval computes a dissimilarity measure between two sets of plasmid bins, that can originate either from two plasmid binning tools, or from a plasmid binning tool and a ground truth set of plasmid bins. The PlasEval dissimilarity accounts for the contig content of plasmid bins, the length of contigs and is repeat-aware. Moreover, the dissimilarity score computed by PlasEval is broken down into several parts, that allows to understand qualitative differences between the compared sets of plasmid bins. We illustrate the use of PlasEval by benchmarking four recently developed plasmid binning tools-PlasBin-flow, HyAsP, gplas, and MOB-recon-on a data set of 53 E. coli bacterial genomes.

CONCLUSION

Analysis of the results of plasmid binning methods using PlasEval shows that their behaviour varies significantly. PlasEval can be used to decide which specific plasmid binning method should be used for a specific dataset. The disagreement between different methods also suggests that the problem of plasmid binning on short-read contigs requires further research. We believe that PlasEval can prove to be an effective tool in this regard. PlasEval is publicly available at https://github.com/acme92/PlasEval.

Comparative evaluation of commercial DNA isolation approaches for nanopore-only bacterial genome assembly and plasmid recovery

The advent of Oxford Nanopore Technologies has undergone significant improvements in terms of sequencing costs, accuracy, and sequencing read lengths, making it a cost-effective, and readily accessible approach for analyzing microbial genomes. A major challenge for bacterial whole genome sequencing by Nanopore technology is the requirement for a higher quality and quantity of high molecular weight DNA compared to short-read sequencing platforms. In this study, using eight pathogenic bacteria, we evaluated the quality, quantity, and fragmented size distribution of extracted DNA obtained from three different commercial DNA extraction kits, and one automated robotic platform. Our results demonstrated significant variation in DNA yield and purity among the extraction kits. The ZymoBIOMICS DNA Miniprep Kit (ZM) provided a higher purity of DNA compared to other kit-based extractions. All kit-based DNA extractions were successfully performed on all twenty-four samples using a single MinION flow cell, with the Nanobind CBB Big DNA kit (NB) yielding the longest raw reads. The Fire Monkey HMW-DNA Extraction Kit (FM) and the automated Roche MagNaPure 96 platform (RO) outperformed in genome assembly, particularly in gram-negative bacteria. Based on our finding, we recommend a minimum read coverage and raw read N50, obtained from the appropriate DNA extraction kit for each bacterial species, to optimize genome assembly and plasmid recovery. This approach will assist end-users in selecting the most effective kit-based extraction method for bacterial whole-genome assembly using only long-read nanopore sequences.

Extraction of high-molecular-weight DNA from Streptococcus spp. for nanopore sequencing in resource-limited settings

The long-read sequencing platform MinION, developed by Oxford Nanopore Technologies, enables the sequencing of bacterial genomes in resource-limited settings, such as field conditions or low- and middle-income countries. For this purpose, protocols for extracting high-molecular-weight DNA using nonhazardous, inexpensive reagents and equipment are needed, and some methods have been developed for gram-negative bacteria. However, we found that without modification, these protocols are unsuitable for gram-positive Streptococcus spp., a major threat to fish farming and food security in low- and middle-income countries. Multiple approaches were evaluated, and the most effective was an extraction method using lysozyme, sodium dodecyl sulfate, and proteinase K for lysis of bacterial cells and magnetic beads for DNA recovery. We optimized the method to consistently achieve sufficient yields of pure high-molecular-weight DNA with minimal reagents and time and developed a version of the protocol which can be performed without a centrifuge or electrical power. The suitability of the method was verified by MinION sequencing and assembly of 12 genomes of epidemiologically diverse fish-pathogenic Streptococcus iniae and Streptococcus agalactiae isolates. The combination of effective high-molecular-weight DNA extraction and MinION sequencing enabled the discovery of a naturally occurring 15 kb low-copy number mobilizable plasmid in S. iniae, which we name pSI1. We expect that our resource-limited settings-adapted protocol for high-molecular-weight DNA extraction could be implemented successfully for similarly recalcitrant-to-lysis gram-positive bacteria, and it represents a method of choice for MinION-based disease diagnostics in low- and middle-income countries.

Sequencing Strategy to Ensure Accurate Plasmid Assembly

Despite the wide use of plasmids in research and clinical production, the need to verify plasmid sequences is a bottleneck that is too often underestimated in the manufacturing process. Although sequencing platforms continue to improve, the method and assembly pipeline chosen still influence the final plasmid assembly sequence. Furthermore, few dedicated tools exist for plasmid assembly, especially for de novo assembly. Here, we evaluated short-read, long-read, and hybrid (both short and long reads) de novo assembly pipelines across three replicates of a 24-plasmid library. Consistent with previous characterizations of each sequencing technology, short-read assemblies had issues resolving GC-rich regions, and long-read assemblies commonly had small insertions and deletions, especially in repetitive regions. The hybrid approach facilitated the most accurate, consistent assembly generation and identified mutations relative to the reference sequence. Although Sanger sequencing can be used to verify specific regions, some GC-rich and repetitive regions were difficult to resolve using any method, suggesting that easily sequenced genetic parts should be prioritized in the design of new genetic constructs.

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.

Nanopore adaptive sampling effectively enriches bacterial plasmids

Bacterial plasmids play a major role in the spread of antibiotic resistance genes. However, their characterization via DNA sequencing suffers from the low abundance of plasmid DNA in those samples. Although sample preparation methods can enrich the proportion of plasmid DNA before sequencing, these methods are expensive and laborious, and they might introduce a bias by enriching only for specific plasmid DNA sequences. Nanopore adaptive sampling could overcome these issues by rejecting uninteresting DNA molecules during the sequencing process. In this study, we assess the application of adaptive sampling for the enrichment of low-abundant plasmids in known bacterial isolates using two different adaptive sampling tools. We show that a significant enrichment can be achieved even on expired flow cells. By applying adaptive sampling, we also improve the quality of de novo plasmid assemblies and reduce the sequencing time. However, our experiments also highlight issues with adaptive sampling if target and non-target sequences span similar regions.

IMPORTANCE

Antimicrobial resistance causes millions of deaths every year. Mobile genetic elements like bacterial plasmids are key drivers for the dissemination of antimicrobial resistance genes. This makes the characterization of plasmids via DNA sequencing an important tool for clinical microbiologists. Since plasmids are often underrepresented in bacterial samples, plasmid sequencing can be challenging and laborious. To accelerate the sequencing process, we evaluate nanopore adaptive sampling as an in silico method for the enrichment of low-abundant plasmids. Our results show the potential of this cost-efficient method for future plasmid research but also indicate issues that arise from using reference sequences.

Genomic mining of Vibrio parahaemolyticus highlights prevalence of antimicrobial resistance genes and new genetic markers associated with AHPND and tdh + /trh + genotypes

BACKGROUND

Acute Hepatopancreatic Necrosis Disease (AHPND) causes significant mortality in shrimp aquaculture. The infection is primarily instigated by Vibrio parahaemolyticus (Vp) strains carrying a plasmid encoding the binary toxin PirAB. Yet, comprehension of supplementary virulence factors associated with this relatively recent disease remains limited. Furthermore, the same holds for gastroenteritis in humans caused by other Vp genotypes. Additionally, given the prevalent use of antibiotics to combat bacterial infections, it becomes imperative to illuminate the presence of antimicrobial resistance genes within these bacteria.

RESULTS

A subsampled number of 1,036 Vp genomes was screened for the presence of antimicrobial resistance genes, revealing an average prevalence of 5 ± 2 (SD) genes. Additional phenotypic antimicrobial susceptibility testing of three Vp strains (M0904, TW01, and PV1) sequenced in this study demonstrated resistance to ampicillin by all tested strains. Additionally, Vp M0904 showed multidrug resistance (against ampicillin, tetracycline, and trimethoprim-sulfamethoxazole). With a focus on AHPND, a screening of all Vibrio spp. for the presence of pirA and/or pirB indicates an estimated prevalence of 0.6%, including four V. campbellii, four V. owensii, and a Vibrio sp. next to Vp. Their pirAB-encoding plasmids exhibited a highly conserved backbone, with variations primarily in the region of the Tn3 family transposase. Furthermore, an assessment of the subsampled Vp genomes for the presence of known virulence factors showed a correlation between the presence of the Type 3 Secretion System 2 and tdh, while the presence of the Type 6 Secretion System 1 was clade dependent. Furthermore, a genome-wide association study (GWAS) unveiled (new) genes associated with pirA, pirB, tdh, and trh genotypes. Notable associations with the pirAB genotype included outer membrane proteins, immunoglobulin-like domain containing proteins, and toxin-antitoxin systems. For the tdh + /trh + genotypes (containing tdh, trh, or both genes), associations were found with T3SS2 genes, urease-related genes and nickel-transport system genes, and genes involved in a 'minimal' type I-F CRISPR mechanism.

CONCLUSIONS

This study highlights the prevalence of antimicrobial resistance and virulence genes in Vp, identifying novel genetic markers associated with AHPND and tdh + /trh + genotypes. These findings contribute valuable insights into the genomic basis of these genotypes, with implications for shrimp aquaculture and food safety.

A practical guide and Galaxy workflow to avoid inter-plasmidic repeat collapse and false gene loss in Unicycler's hybrid assemblies

Generating complete, high-quality genome assemblies is key for any downstream analysis, such as comparative genomics. For bacterial genome assembly, various algorithms and fully automated pipelines exist, which are free-of-charge and easily accessible. However, these assembly tools often cannot unambiguously resolve a bacterial genome, for example due to the presence of sequence repeat structures on the chromosome or on plasmids. Then, a more sophisticated approach and/or manual curation is needed. Such modifications can be challenging, especially for non-bioinformaticians, because they are generally not considered as a straightforward process. In this study, we propose a standardized approach for manual genome completion focusing on the popular hybrid assembly pipeline Unicycler. The provided Galaxy workflow addresses two weaknesses in Unicycler's hybrid assemblies: (i) collapse of inter-plasmidic repeats and (ii) false loss of single-copy sequences. To demonstrate and validate how to detect and resolve these assembly errors, we use two genomes from the Bacillus cereus group. By applying the proposed pipeline following an automated assembly, the genome sequence quality can be significantly improved.

Short- and Long-Read Sequencing Reveals the Presence and Evolution of an IncF Plasmid Harboring blaCTX-M-15 and blaCTX-M-27 Genes in Escherichia coli ST131

Escherichia coli sequence type 131 (ST131) has contributed to the spread of extended-spectrum beta-lactamase (ESBL) and has emerged as the dominant cause of hospital- and community-acquired urinary tract infections. Here, we report for the first time an in-depth analysis of whole-genome sequencing (WGS) of 4 ESBL-producing E. coli ST131 isolates recovered from patients in two hospitals in Armenia using Illumina short-read sequencing for accurate base calling to determine their genotype and to infer their phylogeny and using Oxford Nanopore Technologies long-read sequencing to resolve plasmid and chromosomal genetic elements. Genotypically, the four Armenian isolates were identified as part of the H30Rx/clade C2 (n = 2) and H41/clade A (n = 2) lineages and were phylogenetically closely related to isolates from the European Nucleotide Archive (ENA) database previously recovered from patients in the United States, Australia, and New Zealand. The Armenian isolates recovered in this study had chromosomal integration of the blaCTX-M-15 gene in the H30Rx isolates and a high number of virulence genes found in the H41 isolates associated with the carriage of a rare genomic island (in the context of E. coli ST131) containing the S fimbrial, salmochelin siderophore, and microcin H47 virulence genes. Furthermore, our data show the evolution of the IncF[2:A2:B20] plasmid harboring both blaCTX-M-15 and blaCTX-M-27 genes, derived from the recombination of genes from an IncF[F2:A-:B-] blaCTX-M-15-associated plasmid into the IncF[F1:A2:B20] blaCTX-M-27-associated plasmid backbone seen in two genetically closely related H41 Armenian isolates. IMPORTANCE Combining short and long reads from whole-genome sequencing analysis provided a genetic context for uncommon genes of clinical importance to better understand transmission and evolutionary features of ESBL-producing uropathogenic E. coli (UPEC) ST131 isolates recovered in Armenia. Using hybrid genome assembly in countries lacking genomic surveillance studies can inform us about new lineages not seen in other countries with genes encoding high virulence and antibiotic resistance harbored on mobile genetic elements.

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.