Long read genome assemblers struggle with small plasmids.

Whole bacterial genome sequence reconstruction using Oxford Nanopore Technologies ('Nanopore') long-read-only sequencing may offer a lower-cost, higher-throughput alternative for pathogen surveillance to 'hybrid' assembly with recent improvements in Nanopore sequencing accuracy. We evaluated the accuracy, including plasmid reconstruction, of Nanopore long-read-only genome assemblies of Enterobacterales. We sequenced 92 genomes from clinical Enterobacterales isolates, collected in England under a national surveillance programme, with long-read Nanopore (R10.4.1, Dorado v5.0.0 super-high-accuracy basecalled) and short-read Illumina (NovaSeq) sequencing approaches. Genomes were assembled using three long-read-only (Flye, Hybracter long and Autocycler) and three hybrid assemblers (Hybracter hybrid, Unicycler normal and bold). Three polishing modalities (Medaka v2 with subsampled or un-subsampled long-reads; Polypolish+Pypolca with short-reads) were investigated. Autocycler circularised the most chromosomes [87/92 (95%)]. Plasmid sequence reconstruction was comparable among all assemblers except Flye, all recovering 90-96% of plasmids, although the 'ground truth' was uncertain. Flye performed worse than other assemblers on almost all metrics. Autocycler+Medaka (un-subsampled long-reads) was the most accurate long-read-only assembler/polisher combination, comparable to hybrid assemblies [median 0 (IQR: 0-0) single nucleotide variants (SNVs) and 0 (IQR: 0-1) insertions/deletions (indels) per genome; median quality value/Q score 100 (IQR: 64-100)], with only 4/92 genome sequences having >10 SNVs/indels. Medaka polishing with un-subsampled long-reads resulted in small improvements in indels, but not SNVs for both Flye and Autocycler assemblies. Seven-locus multi-locus sequence type, antimicrobial resistance, virulence and stress gene annotation was equivalent across assembler/polisher combinations. Nanopore long-read-only bacterial genome assembly with Autocycler combined with Medaka polishing (using un-subsampled reads) is similarly accurate and possibly more complete than hybrid assemblies, representing a viable alternative for incorporating high-quality genomic data, including plasmids, into Enterobacterales surveillance.

Recent advances and lower costs in rapid high-throughput sequencing have engendered hope that whole genome sequencing (WGS) might afford complete resistome characterization in bacterial isolates. WGS is particularly useful for the clinical characterization of fastidious and slow-growing bacteria. Despite its potential, several challenges should be addressed before adopting WGS to detect antimicrobial resistance (AMR) genes in the clinical laboratory. Here, with three distinct ESKAPE bacteria (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter spp.), different approaches were compared to identify best practices for detecting AMR genes, including: total genomic DNA and plasmid DNA extractions, the solo assembly of Illumina short-reads and of Oxford Nanopore Technologies (ONT) long-reads, two hybrid assembly pipelines, and three in silico AMR databases. We also determined the susceptibility of each strain to 21 antimicrobials. We found that all AMR genes detected in pure plasmid DNA were also detectable in total genomic DNA, indicating that, at least in these three enterobacterial genera, the purification of plasmid DNA was not necessary to detect plasmid-borne AMR genes. Illumina short-reads used with ONT long-reads in either hybrid or polished assemblies of total genomic DNA enhanced the sensitivity and accuracy of AMR gene detection. Phenotypic susceptibility closely corresponded with genotypes identified by sequencing; however, the three AMR databases differed significantly in distinguishing mobile dedicated AMR genes from non-mobile chromosomal housekeeping genes in which rare spontaneous resistance mutations might occur. This study indicates that each method employed in a WGS workflow has an impact on the detection of AMR genes. A combination of short- and long-reads, followed by at least three different AMR databases, should be used for the consistent detection of such genes. Further, an additional step for plasmid DNA purification and sequencing may not be necessary. This study reveals the need for standardized biochemical and informatic procedures and database resources for consistent, reliable AMR genotyping to take full advantage of WGS in order to expedite patient treatment and track AMR genes within the hospital and community.

Over the last ten years, genome sequencing capabilities have expanded exponentially. There have been tremendous advances in sequencing technology, DNA sample preparation, genome assembly, and data analysis. This has led to advances in a number of facets of bacterial genomics, including metagenomics, clinical medicine, bacterial archaeology, and bacterial evolution. This review examines the strengths and weaknesses of techniques in bacterial genome sequencing, upcoming technologies, and assembly techniques, as well as highlighting recent studies that highlight new applications for bacterial genomics.

Zancudomyces culisetae is an obligate symbiotic fungus inhabiting the digestive tracts of aquatic insect larvae, including black flies, midges, and mosquitoes. With a global distribution and high prevalence in disease-transmitting insects, Z. culisetae serves as a model for studying insect gut fungi. A previous draft genome assembly using Illumina short reads provided insights into its genome composition, such as a low GC ratio and evidence of horizontal gene transfer. However, its fragmented nature has limited deeper exploration of the evolutionary mechanisms shaping these gut symbionts. To address this gap, we generated a wealth of genomic resources for Z. culisetae using multiple sequencing platforms, including Illumina, Oxford Nanopore, PacBio-CLR (Complete Long Reads), and PacBio-HiFi (High Fidelity). This also provides an opportunity to compare these popular sequencing methods to suggest the optimal approach for fungal genome assembly. Our results suggest that PacBio-HiFi produced the most complete assembly, yielding a 27.8 Mb genome size with 26 contigs, representing the highest-quality genome of insect gut fungi to date. Additionally, we generated transcriptomic data to support genome annotation, identifying 8,484 protein-coding genes. Despite the improved genome quality, Z. culisetae lacks ∼20% of Benchmarking Universal Single-Copy Orthologue commonly found in fungi, reflecting adaptations to its obligate symbiotic lifestyle. This study not only provides valuable genomic resources for insect gut fungal research but also evaluates the strengths and limitations of current genome sequencing and assembly approaches, offering best practices for fungal genome analysis and genetic research.

The inbred Babraham pig serves as a valuable biomedical model for research due to its high level of homozygosity, including in the major histocompatibility complex (MHC) loci and likely other important immune-related gene complexes, which are generally highly diverse in outbred populations. As the ability to control for this diversity using inbred organisms is of great utility, we sought to improve this resource by generating a long-read whole genome assembly and transcriptome atlas of a Babraham pig. The genome was de novo assembled using PacBio long reads and error-corrected using Illumina short reads. Assembled contigs were then mapped to the porcine reference assembly, Sscrofa11.1, to generate chromosome-level scaffolds. The resulting TPI_Babraham_pig_v1 assembly is nearly as contiguous as Sscrofa11.1 with a contig N50 of 34.95 Mb and contig L50 of 23. The remaining sequence gaps are generally the result of poor assembly across large and highly repetitive regions such as the centromeres and tandemly duplicated gene families, including immune-related gene complexes, that often vary in gene content between haplotypes. We also further confirm homozygosity across the Babraham MHC and characterize the allele content and tissue expression of several other immune-related gene complexes, including the antibody and T cell receptor loci, the natural killer complex, and the leukocyte receptor complex. The Babraham pig genome assembly provides an alternate highly contiguous porcine genome assembly as a resource for the livestock genomics community. The assembly will also aid biomedical and veterinary research that utilizes this animal model such as when controlling for genetic variation is critical.

StrainSelect: A novel microbiome reference database that disambiguates all bacterial strains, genome assemblies and extant cultures worldwide

Oxford Nanopore sequencing can be used to achieve complete bacterial genomes. However, the error rates of Oxford Nanopore long reads are greater compared to Illumina short reads. Long-read assemblers using a variety of assembly algorithms have been developed to overcome this deficiency, which have not been benchmarked for genomic analyses of bacterial pathogens using Oxford Nanopore long reads. In this study, long-read assemblers, namely Canu, Flye, Miniasm/Racon, Raven, Redbean, and Shasta, were thus benchmarked using Oxford Nanopore long reads of bacterial pathogens. Ten species were tested for mediocre- and low-quality simulated reads, and 10 species were tested for real reads. Raven was the most robust assembler, obtaining complete and accurate genomes. All Miniasm/Racon and Raven assemblies of mediocre-quality reads provided accurate antimicrobial resistance (AMR) profiles, while the Raven assembly of Klebsiella variicola with low-quality reads was the only assembly with an accurate AMR profile among all assemblers and species. All assemblers functioned well for predicting virulence genes using mediocre-quality and real reads, whereas only the Raven assemblies of low-quality reads had accurate numbers of virulence genes. Regarding multilocus sequence typing (MLST), Miniasm/Racon was the most effective assembler for mediocre-quality reads, while only the Raven assemblies of Escherichia coli O157:H7 and K. variicola with low-quality reads showed positive MLST results. Miniasm/Racon and Raven were the best performers for MLST using real reads. The Miniasm/Racon and Raven assemblies showed accurate phylogenetic inference. For the pan-genome analyses, Raven was the strongest assembler for simulated reads, whereas Miniasm/Racon and Raven performed the best for real reads. Overall, the most robust and accurate assembler was Raven, closely followed by Miniasm/Racon.

Halisarca dujardinii is a marine sponge known for its ability to completely regenerate via cell reaggregation. Here we present the first chromosome-level genome assembly of H. dujardinii, generated using Oxford Nanopore long reads, Illumina short reads, and Hi-C data. The final assembly spans 226.5 Mbp and is organized into 21 chromosome-scale scaffolds, representing the full nuclear genome, with an assembly N50 of 10.1 Mbp. In addition, we report the complete mitochondrial genome for H. dujardinii, assembled as a circular molecule and annotated to contain 14 protein-coding genes. We provide a comprehensive genome annotation comprising 14,565 nuclear protein-coding genes, of which 85.6% are functionally annotated, along with repetitive elements and non-coding RNAs. Transcript models were refined using extensive bulk and single-cell RNA sequencing data, enabling accurate annotation of 3' untranslated regions. This new genomic resource provides a valuable foundation for investigating the molecular basis of sponge regeneration, as well as for comparative studies of sponge evolution and marine biology.

Extended laboratory culture and antimicrobial susceptibility testing timelines hinder rapid species identification and susceptibility profiling of bacterial pathogens associated with bovine respiratory disease, the most prevalent cause of cattle mortality in the United States. Whole-genome sequencing offers a culture-independent alternative to current bacterial identification methods, but requires a library of bacterial reference genomes for comparison. To contribute new bacterial genome assemblies and evaluate genetic diversity and variation in antimicrobial resistance genotypes, whole-genome sequencing was performed on bovine respiratory disease–associated bacterial isolates (Histophilus somni, Mycoplasma bovis, Mannheimia haemolytica, and Pasteurella multocida) from dairy and beef cattle. One hundred genomically distinct assemblies were added to the NCBI database, doubling the available genomic sequences for these four species. Computer-based methods identified 11 predicted antimicrobial resistance genes in three species, with none being detected in M. bovis. While computer-based analysis can identify antibiotic resistance genes within whole-genome sequences (genotype), it may not predict the actual antimicrobial resistance observed in a living organism (phenotype). Antimicrobial susceptibility testing on 64 H. somni, M. haemolytica, and P. multocida isolates had an overall concordance rate between genotype and phenotypic resistance to the associated class of antimicrobials of 72.7% (P < 0.001), showing substantial discordance. Concordance rates varied greatly among different antimicrobial, antibiotic resistance gene, and bacterial species combinations. This suggests that antimicrobial susceptibility phenotypes are needed to complement genomically predicted antibiotic resistance gene genotypes to better understand how the presence of antibiotic resistance genes within a given bacterial species could potentially impact optimal bovine respiratory disease treatment and morbidity/mortality outcomes.

The genome sequence of multidrug-resistant Raoultella terrigena RT01-5M1 strain isolated from Canadian farmed salmon was determined using Oxford nanopore and Illumina MiSeq sequencers. The assembled chromosome was estimated at 5,699,993 bp in size, with two plasmids, 164,879 bp and 82,046 bp. The chromosome and smaller plasmid contained antimicrobial resistance genes.

BackgroundHuman cytomegalovirus (HCMV) has a double-stranded DNA genome of approximately 235 Kbp that is structurally complex including extended GC-rich repeated regions. Genomic recombination events are frequent in HCMV cultures but have also been observed in vivo. Thus, the assembly of HCMV whole genomes from technologies producing shorter than 500 bp sequences is technically challenging. Here we improved the reconstruction of HCMV full genomes by means of a hybrid, de novo genome-assembly bioinformatics pipeline upon data generated from the recently released MinION MkI B sequencer from Oxford Nanopore Technologies.ResultsThe MinION run of the HCMV (strain TB40/E) library resulted in ~ 47,000 reads from a single R9 flowcell and in ~ 100× average read depth across the virus genome. We developed a novel, self-correcting bioinformatics algorithm to assemble the pooled HCMV genomes in three stages. In the first stage of the bioinformatics algorithm, long contigs (N50 = 21,892) of lower accuracy were reconstructed. In the second stage, short contigs (N50 = 5686) of higher accuracy were assembled, while in the final stage the high quality contigs served as template for the correction of the longer contigs resulting in a high-accuracy, full genome assembly (N50 = 41,056). We were able to reconstruct a single representative haplotype without employing any scaffolding steps. The majority (98.8%) of the genomic features from the reference strain were accurately annotated on this full genome construct. Our method also allowed the detection of multiple alternative sub-genomic fragments and non-canonical structures suggesting rearrangement events between the unique (UL /US) and the repeated (T/IRL/S) genomic regions.ConclusionsThird generation high-throughput sequencing technologies can accurately reconstruct full-length HCMV genomes including their low-complexity and highly repetitive regions. Full-length HCMV genomes could prove crucial in understanding the genetic determinants and viral evolution underpinning drug resistance, virulence and pathogenesis.

Whole genome sequencing provides the highest resolution for characterizing pathogen evolution, epidemiology, and diagnostics. Genome assemblies contain information on the identity and potential phenotypes of a pathogen. Likewise, variant calling can inform on transmission patterns and evolutionary relationships. Recent improvements in Oxford Nanopore long-read sequencing have made its use attractive for genomic epidemiology. However, the accuracy and optimal strategy for analysis of Nanopore reads remains to be determined. We compared the use of Illumina short reads and Oxford Nanopore long reads for genome assembly and variant calling of phytopathogenic bacteria. We generated short- and long-read datasets for diverse phytopathogenic Agrobacterium strains. We then analyzed these data using multiple pipelines designed for either short or long reads and compared the results. We found that assemblies made from long reads were more complete than those made from short-read data and contained few sequence errors. Variant calling pipelines differed in their ability to accurately call variants and infer genotypes from long reads. Results suggest that computationally fragmenting long reads can improve the accuracy of variant calling in population-level studies. Using fragmented long reads, pipelines designed for short reads were more accurate at recovering genotypes than pipelines designed for long reads. Further, short- and long-read datasets can be analyzed together with the same pipelines. These findings show that Oxford Nanopore sequencing is accurate and can be sufficient for microbial pathogen genomics and epidemiology. Ultimately, this enhances the ability of researchers and clinicians to understand and mitigate the spread of pathogens.

Italian ryegrass (Lolium multiflorum Lam.) is an important forage grass, providing a major source of feed for ruminants in temperate regions. Due to its highly heterozygous and repeat‐rich genome, high‐quality chromosome‐level genome assemblies are scarce for Italian ryegrass. Here, we sequenced the genome of a genotype from the Italian ryegrass cultivar ‘Rabiosa’ (hereafter referred to as Rabiosa), and we obtained Oxford Nanopore Technologies long reads (∼60‐fold coverage), Illumina short reads (∼85‐fold coverage), and high‐throughput chromosome conformation capture data (∼60‐fold coverage). With Rabiosa as one of the parents, we constructed an F1 population consisting of 304 individuals, which were genotyped by reduced representation sequencing for linkage map construction and quantitative trait locus (QTL) analysis. Using whole‐genome sequencing data of Rabiosa and the genetic linkage map, we first generated a chromosome‐level unphased haploid assembly (scaffold N50 of 338.75 Mb, total Benchmarking Universal Single‐Copy Orthologs [BUSCO] score of 94.60%). Then, based on the unphased assembly and a reference‐based phasing approach, we generated a chromosome‐level haplotype‐resolved assembly containing both haplotypes (scaffold N50 of ∼250 Mb and total BUSCO score of ∼90% for each haplome). Between the two haplotypes of Rabiosa, we observed a highly collinear gene order at the chromosome level and a high sequence variation at the local level. With a graph‐based reference built from the unphased and the haplotype‐resolved assemblies of Rabiosa, we conducted a QTL analysis, and two QTL significantly associated with stem rust resistance were detected. The genome assemblies of Rabiosa will serve as invaluable genomic resources to facilitate genomic applications in forage grass research and breeding.

Edwardsiella piscicida is a pathogenic bacterium responsible for significant losses in important wild and cultured fish species. E. piscicida strain MS-18-199 recovered from a diseased hybrid catfish from East Mississippi and showed resistance to florfenicol, chloramphenicol, oxytetracycline, doxycycline, erythromycin, tetracycline, azitromycin, spectinomycin, sulfonamide, and bacitracin. To explore the mechanisms of resistance in E. piscicida strain MS-18-199, genomic DNA was extracted and subjected to whole genome sequencing (WGS) using a combination of long (Oxford Nanopore) and short (Illumina) reads. The genome of strain MS-18-199 revealed a novel plasmid named pEPMS-18199. The 117,448 bp plasmid contains several antimicrobial resistance (AMR) elements/genes, including florfenicol efflux pump (floR), tetracycline efflux pump (tetA), tetracycline repressor protein (tetR), sulfonamide resistance (sul2), aminoglycoside O-phosphotransferase aph(6)-Id (strB), and aminoglycoside O-phosphotransferase aph(3)-Ib (strA). Two genes, arsA and arsD, that encode protein components related to transport/resistance to arsenic were also found in pEPMS-18199. In addition, pEPMS-18199 carried twelve conjugative transfer genes (tra), eight transposases and insertion elements, two plasmid stability proteins, two replication proteins, and three partitioning proteins (par system). Results from mobilization and stability experiments revealed that pEPMS-18199 is highly stable in the host cell and could be transferred to Escherichia coli and Edwardsiella ictaluri by conjugation. To our knowledge, this is the first detection of a multidrug resistance (MDR) conjugative plasmid in E. piscicida in the United States. Careful tracking of this plasmid in the aquaculture system is warranted. Knowledge regarding the molecular mechanisms of AMR in aquaculture is important for antimicrobial stewardship.

High-quality genome assembly of the silkworm, Bombyx mori