An efficient error correction and accurate assembly tool for noisy long reads

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.

The advent of DNA and RNA sequencing has revolutionized the study of genomics and molecular biology. Next generation sequencing (NGS) technologies like Illumina, Ion Torrent, SOLiD sequencing etc. have brought about a quick and cheap way to sequence genomes. Recently, third generation sequencing (TGS) technologies like PacBio and Oxford Nanopore Technology (ONT) have also been developed. Different technologies use different underlying methods for sequencing and are prone to different error rates. Though many tools exist for error correction of sequencing data from NGS and TGS methods, no standard method is available yet to evaluate the accuracy and effectiveness of these error-correction tools. In this study, we present a Software Package for Error Correction Tool Assessment on nuCLEic acid sequences (SPECTACLE) providing comprehensive algorithms to evaluate error-correction methods for DNA and RNA sequencing, for NGS and TGS platforms. We also present a compilation of sequencing datasets for Illumina, PacBio and ONT platforms that present challenging scenarios for error-correction tools. Using these datasets and SPECTACLE, we evaluate the performance of 23 different error-correction tools and present unique and helpful insights into their strengths and weaknesses. We hope that our methodology will standardize the evaluation of DNA and RNA error-correction tools in the future.

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.

ABSTRACTMultiple displacement amplification (MDA) outperforms conventional PCR in long fragment and whole‐genome amplification, making it attractive to couple MDA with long‐read sequencing of samples with limited quantities of DNA to obtain improved genome assemblies. Here, we explore the efficacy and limits of MDA for efficient low‐cost genome sequence assembly using Oxford Nanopore Technologies (ONTs) rapid library preparations and minION sequencing. We successfully generated almost complete genome sequences for all organisms examined, including Gram‐positive (Staphylococcus aureus, Enterococcus faecium) and Gram‐negative (Escherichia coli) prokaryotes and one challenging eukaryotic pathogen (Cryptosporidium spp) representing a broad spectrum of critical infectious disease pathogens. High‐quality data from those samples were generated starting with only 0.025 ng of total DNA. Controlled sheared DNA samples exhibited a distinct pattern of size increase after MDA, which may be associated with the amplification of long, low‐abundance fragments present in the assay, as well as generating concatemeric sequences during amplification. To address concatemers, we developed a computational pipeline (CADECT: Concatemer Detection Tool) to identify and remove putative concatemeric sequences. This study highlights the efficacy of MDA in generating high‐quality genome assemblies from limited amounts of input DNA. Also, the CADECT pipeline effectively mitigated the impact of concatemeric sequences, enabling the assembly of contiguous sequences even in cases where the input genomic DNA was degraded. These results have significant implications for the study of organisms that are challenging to culture in vitro, such as Cryptosporidium, and for expediting critical results in clinical settings with limited quantities of available genomic DNA.

Globally, commercialized plum cultivars are mostly diploid Chinese plums (Prunus salicina Lindl.), also known as Japanese plums, and are one of the most abundant and variable fruit tree species. To advance Prunus genomic research, we present a chromosome-scale P. salicina genome assembly, constructed using an integrated strategy that combines Illumina, Oxford Nanopore, and high-throughput chromosome conformation capture (Hi-C) sequencing. The high-quality genome assembly consists of a 318.6-Mb sequence (contig N50 length of 2.3 Mb) with eight pseudo-chromosomes. The expansion of the P. salicina genome is led by recent segmental duplications and a long terminal repeat burst of approximately 0.2 Mya. This resulted in a significant expansion of gene families associated with flavonoid metabolism and plant resistance, which impacted fruit flavor and increased species adaptability. Population structure and domestication history suggest that Chinese plum may have originated from South China and provides a domestication route with accompanying genomic variations. Selection sweep and genetic diversity analysis enabled the identification of several critical genes associated with flowering time, stress tolerance, and flavonoid metabolism, demonstrating the essential roles of related pathways during domestication. Furthermore, we reconstructed and exploited flavonoid-anthocyanin metabolism using multi-omics analysis in Chinese plum and proposed a complete metabolic pathway. Collectively, our results will facilitate further candidate gene discovery for important agronomic traits in Chinese plum and provide insights into future functional genomic studies and DNA-informed breeding.

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.

Graph states are the main computational building blocks of measurement-based computation and a useful tool for error correction in the gate model architecture. The graph states form a class of quantum states which are eigenvectors for the abelian group of stabilizer operators. They own topological properties, arising from their graph structure, including the presence of highly connected nodes, called hubs. Starting from hub nodes, we demonstrate how to efficiently express a graph state through the generators of the stabilizer group. We provide examples by expressing the ring and the star topology, for which the number of stabilizers reduces from <tex xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">$n$</tex> to <tex xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">$\lceil\frac{n}{2}\rceil$</tex> , and from <tex xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">$n$</tex> to 1, respectively. We demonstrate that the graph states can be generated by a subgroup of the stabilizer group. Therefore, we provide an algebraic framework to manipulate the graph states with a reduced number of stabilizers.

Root-knot nematodes (Meloidogyne spp.) are among the most destructive agricultural pests that cause significant yield losses across a wide range of crops. Meloidogyne hapla is a valuable model for studying root-knot nematodes due to its parasitic diversity, small diploid genome, and a reproductive strategy that facilitates genetic analysis. Here, we report the most contiguous genome assembly to date for any plant-parasitic nematode built using PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C sequencing. Genetic linkage analysis of F2 populations derived from crosses between M. hapla strains validated the assembly but also revealed anomalies indicating chromosome structure differences between parental isolates such as fissions, fusions, and rearrangements. Strikingly, we identified sharply delimited zones with extraordinarily high recombination on most chromosomes. Notably, several of these high recombination zones were significantly enriched for genes encoding secreted proteins, many of which contribute to parasitism. These findings suggest that meiotic recombination facilitates effector diversification and offer insight into how these parasites diversify their effector protein repertoire to change or expand their extraordinary host range. We further report the discovery of a novel 16-nucleotide tandem repeat and lack of canonical telomere repeats at chromosome ends. The localization of this 16-nt repeat at chromosome ends highlights a potentially divergent mechanism of chromosome-end maintenance in this nematode group. Overall, our study integrates high-resolution structural genomics, genetic mapping, and functional inference to uncover links between genome architecture, recombination landscapes, and host-parasite interactions.

DNA sequencing technology is undergoing a revolution with the commercialization of next genera-tion technologies. Over the past eight to ten years massively parallel DNA sequencing platforms have become widely available with reducing the cost of DNA sequencing. Next generation platform (NGS) includes Helicos Heliscope™, Pacific Biosciences SMRT, Ion Torrent, Oxford Nanopore, etc. These plat-forms have the potential to dramatically accelerate biological research, by enabling the comprehensive analysis of genomes, transcriptomes and interactomes to become inexpensive, routine and widespread. Variant discovery by re-sequencing targeted regions of interest or whole genomes, de novo assemblies of prokaryotic and eukaryotic genomes. Cataloguing the transcriptomes of cells, genome-wide profiling of epigenetic marks and chromatin structure is using other seq-based methods and species classification and gene discovery by metagenomics studies.

Genomic resources have yielded unprecedented insights into ecological and evolutionary processes, not to mention their importance in economic and conservation management of specific organisms. However, the field of macroalgal genomics is hampered by difficulties in the isolation of suitable DNA. Even when DNA that appears high quality by standard metrics has been isolated, such samples may not perform well during the sequencing process. We here have compared Oxford Nanopore long-read sequencing results for three species of macroalgae to those of nonmacroalgal species and determined that when using macroalgal samples, sequencing activity declined rapidly, resulting in reduced sequencing yield. Chemical analysis of macroalgal DNA that would be considered suitable for sequencing revealed that DNA derived from dried macroalgae was enriched for polyphenol-DNA adducts (DNA with large polyphenols chemically attached to it), which may have led to sequencing inhibition. Of note, we observed the strongest evidence of sequencing inhibition and reduced sequence output when using samples dried using silica gel-suggesting that such storage approaches may not be appropriate for samples destined for Oxford Nanopore sequencing. Our findings have wide-ranging implications for the generation of genomic resources from macroalgae and suggest a need to develop new storage methods that are more amenable to Oxford Nanopore sequencing or to use fresh flash-frozen tissue wherever possible for genome sequencing.

For benchmarking, development and evaluation of tools we have sequenced the Yoruban reference genome NA19240 to 80x coverage using Oxford Nanopore PromethION. Additionally, we are developing a workflow for efficient detection and annotation of structural variation from long read sequencing data, integrating tools such as Sniffles and NanoSV.

A high quality genome assembly is a vital first step for the study of an organism. Recent advances in technology have made the creation of high quality chromosome scale assemblies feasible and low cost. However, the amount of input DNA needed for an assembly project can be a limiting factor for small organisms or precious samples. Here we demonstrate the feasibility of creating a chromosome scale assembly using a hybrid method for a low input sample, a single outbred Drosophila melanogaster. Our approach combines an Illumina shotgun library, Oxford nanopore long reads, and chromosome conformation capture for long range scaffolding. This single fly genome assembly has a N50 of 26 Mb, a length that encompasses entire chromosome arms, contains 95% of expected single copy orthologs, and a nearly complete assembly of this individual's Wolbachia endosymbiont. The methods described here enable the accurate and complete assembly of genomes from small, field collected organisms as well as precious clinical samples.

In recent years, both long-read sequencing and metagenomic analysis have been significantly advanced. Although long-read sequencing technologies have been primarily used for de novo genome assembly, they are rapidly maturing for widespread use in other applications. In particular, long reads could potentially lead to more precise taxonomic identification which has sparked an interest in using them for metagenomic analysis. Here we present a benchmark of several tools for metagenomic taxonomic classification, tested on in-silico datasets that were constructed using real long reads from isolate sequencing. We compared tools that were either newly developed for or modified to work with long reads, including Kraken, Centrifuge, CLARK, MetaMaps and MEGAN-LR. The test datasets were constructed with varying numbers of bacterial and eukaryotic genomes, to simulate different metagenomic applications. The tools were tested on their ability to accurately detect species and precisely estimate species abundances in the samples. Our analysis showed that all tested classifiers provide useful results, and that accuracy was strongly influenced by the comprehensiveness of the default database used. Using the same database for all tools provided comparable results across methods except for MetaMaps which had slightly better performance, but was slower than k-mer based tools.

Whole-genome sequencing has become a preferred method for studying bacterial plasmids, as it is generally assumed to capture the entire genome. However, long-read genome assemblers have been shown to sometimes miss plasmid sequences - an issue that has been associated with plasmid size. The purpose of this study was to investigate the relationship between plasmid size and plasmid recovery by the long-read-only assemblers, Flye, Raven, Miniasm, and Canu. This was accomplished by determining the number of times each assembler successfully recovered 33 plasmids, ranging from 1919 to 194 062 bp in size and belonging to 14 bacterial isolates from six bacterial genera, using Oxford Nanopore long reads. These results were additionally compared to plasmid recovery rates by the short-read-first assembler, Unicycler, using both Oxford Nanopore long reads and Illumina short reads. Results from this study indicate that Canu, Flye, Miniasm, and Raven are prone to missing plasmid sequences, whereas Unicycler was successful at recovering 100 % of plasmid sequences. Excluding Canu, most plasmid loss by long-read-only assemblers was due to failure to recover plasmids smaller than 10 kb. As such, it is recommended that Unicycler be used to increase the likelihood of plasmid recovery during bacterial genome assembly.

Piper nigrum and Piper longum are the most popular and economically essential spice crops globally valued for their aromatic alkaloids, especially Piperine. However, Piperine synthesis pathway mechanisms are not yet well known. This work was aimed to generate the full-length comparative berry transcriptome analysis dataset of P. nigrum and P. longum by Illumina and Nanopore sequencing platforms. While short-read sequencing technology is widely using to capture transcriptome profiles, there are still some limitations due to the read length. We used Oxford Nanopore technology for long reads and the Illumina sequencing platform for short reads to generate a hybrid transcriptome assembly from half matured and fully matured berries of P. nigrum and P. longum. From P. nigrum and P. longum 37.3 million and 38.1 million raw reads were generated respectively. A total of 308369 contigs from P. nigrum and 267715 contigs from P. longum were obtained and successfully annotated. The transcriptome data revealed gene families involved in piperine and other secondary metabolite biosynthetic pathways. The raw data were uploaded to NCBI database. This dataset shed light on the further exploration of the piperine biosynthetic pathway, its transcriptomic changes, and evolution. Data generated has been submitted to SRA of NCBI with Bio samples accession: (SAMN13981803, SAMN22826456).