scholarly journals Evaluation of hybrid and non-hybrid methods for de novo assembly of nanopore reads

2015 ◽  
Author(s):  
Ivan Sovic ◽  
Kresimir Krizanovic ◽  
Karolj Skala ◽  
Mile Sikic
Keyword(s):  
De Novo Assembly ◽  
De Novo ◽  
Hybrid Methods ◽  
Bacterial Genome ◽  
Error Rates ◽  
Sequencing Data ◽  
E Coli ◽  
Recent Emergence ◽  
K 12

Recent emergence of nanopore sequencing technology set a challenge for the established assembly methods not optimized for the combination of read lengths and high error rates of nanopore reads. In this work we assessed how existing de novo assembly methods perform on these reads. We benchmarked three non-hybrid (in terms of both error correction and scaffolding) assembly pipelines as well as two hybrid assemblers which use third generation sequencing data to scaffold Illumina assemblies. Tests were performed on several publicly available MinION and Illumina datasets of E. coli K-12, using several sequencing coverages of nanopore data (20x, 30x, 40x and 50x). We attempted to assess the quality of assembly at each of these coverages, to estimate the requirements for closed bacterial genome assembly. Results show that hybrid methods are highly dependent on the quality of NGS data, but much less on the quality and coverage of nanopore data and perform relatively well on lower nanopore coverages. Furthermore, when coverage is above 40x, all non-hybrid methods correctly assemble the E. coli genome, even a non-hybrid method tailored for Pacific Bioscience reads. While it requires higher coverage compared to a method designed particularly for nanopore reads, its running time is significantly lower.

scholarly journals When Less is More: "Slicing" Sequencing Data Improves Read Decoding Accuracy and De Novo Assembly Quality

2015 ◽  
Author(s):  
Stefano Lonardi ◽  
Hamid Mirebrahim ◽  
Steve Wanamaker ◽  
Matthew Alpert ◽  
Gianfranco Ciardo ◽  
...  
Keyword(s):  
Deep Sequencing ◽  
De Novo Assembly ◽  
De Novo ◽  
Optimal Size ◽  
Sequencing Data ◽  
Less Is More ◽  
Bac Clones ◽  
Deep Sequencing Data ◽  
First Time

Since the invention of DNA sequencing in the seventies, computational biologists have had to deal with the problem de novo genome assembly with limited (or insufficient) depth of sequencing. In this work, for the first time we investigate the opposite problem, that is, the challenge of dealing with excessive depth of sequencing. Specifically, we explore the effect of ultra-deep sequencing data in two domains: (i) the problem of decoding reads to BAC clones (in the context of the combinatorial pooling design proposed by our group), and (ii) the problem of de novo assembly of BAC clones. Using real ultra-deep sequencing data, we show that when the depth of sequencing increases over a certain threshold, sequencing errors make these two problems harder and harder (instead of easier, as one would expect with error-free data), and as a consequence the quality of the solution degrades with more and more data. For the first problem, we propose an effective solution based on "divide and conquer": we "slice" a large dataset into smaller samples of optimal size, decode each slice independently, then merge the results. Experimental results on over 15,000 barley BACs and over 4,000 cowpea BACs demonstrate a significant improvement in the quality of the decoding and the final assembly. For the second problem, we show for the first time that modern de novo assemblers cannot take advantage of ultra-deep sequencing data.

scholarly journals A complete bacterial genome assembled de novo using only nanopore sequencing data

2015 ◽  
Author(s):  
Nicholas James Loman ◽  
Joshua Quick ◽  
Jared T Simpson
Keyword(s):  
De Novo ◽  
Bacterial Genome ◽  
Reference Sequence ◽  
Nucleotide Identity ◽  
Sequencing Data ◽  
Celera Assembler ◽  
Level Data ◽  
Complete Bacterial Genome ◽  
Sequencing Platforms ◽  
K 12

A method for de novo assembly of data from the Oxford Nanopore MinION instrument is presented which is able to reconstruct the sequence of an entire bacterial chromosome in a single contig. Initially, overlaps between nanopore reads are detected. Reads are then subjected to one or more rounds of error correction by a multiple alignment process employing partial order graphs. After correction, reads are assembled using the Celera assembler. Finally, the assembly is polished using signal-level data from the nanopore employing a novel hidden Markov model. We show that this method is able to assemble nanopore reads from Escherichia coli K-12 MG1655 into a single contig of length 4.6Mb permitting a full reconstruction of gene order. The resulting draft assembly has 98.4% nucleotide identity compared to the finished reference genome. After polishing the assembly with our signal-level HMM, the nucleotide identity is improved to 99.4%. We show that MinION sequencing data can be used to reconstruct genomes without the need for a reference sequence or data from other sequencing platforms.

scholarly journals Comparison of long-read sequencing technologies in interrogating bacteria and fly genomes

2021 ◽  
Author(s):  
Eric S Tvedte ◽  
Mark Gasser ◽  
Benjamin C Sparklin ◽  
Jane Michalski ◽  
Carl E Hjelmen ◽  
...  
Keyword(s):  
Bacterial Genome ◽  
Hybrid Approach ◽  
Cost Effective ◽  
Fruit Fly ◽  
Drosophila Ananassae ◽  
Whole Genome Sequencing Data ◽  
Sequencing Data ◽  
E Coli ◽  
Hybrid Approaches ◽  
Long Read

Abstract The newest generation of DNA sequencing technology is highlighted by the ability to generate sequence reads hundreds of kilobases in length. Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) have pioneered competitive long read platforms, with more recent work focused on improving sequencing throughput and per-base accuracy. We used whole-genome sequencing data produced by three PacBio protocols (Sequel II CLR, Sequel II HiFi, RS II) and two ONT protocols (Rapid Sequencing and Ligation Sequencing) to compare assemblies of the bacteria Escherichia coli and the fruit fly Drosophila ananassae. In both organisms tested, Sequel II assemblies had the highest consensus accuracy, even after accounting for differences in sequencing throughput. ONT and PacBio CLR had the longest reads sequenced compared to PacBio RS II and HiFi, and genome contiguity was highest when assembling these datasets. ONT Rapid Sequencing libraries had the fewest chimeric reads in addition to superior quantification of E. coli plasmids versus ligation-based libraries. The quality of assemblies can be enhanced by adopting hybrid approaches using Illumina libraries for bacterial genome assembly or polishing eukaryotic genome assemblies, and an ONT-Illumina hybrid approach would be more cost-effective for many users. Genome-wide DNA methylation could be detected using both technologies, however ONT libraries enabled the identification of a broader range of known E. coli methyltransferase recognition motifs in addition to undocumented D. ananassae motifs. The ideal choice of long read technology may depend on several factors including the question or hypothesis under examination. No single technology outperformed others in all metrics examined.

scholarly journals Facile, High Quality Sequencing of Bacterial Genomes from Small Amounts of DNA

2014 ◽  
Vol 2014 ◽  
pp. 1-8
Author(s):  
Momchilo Vuyisich ◽  
Ayesha Arefin ◽  
Karen Davenport ◽  
Shihai Feng ◽  
Cheryl Gleasner ◽  
...  
Keyword(s):  
Genomic Dna ◽  
De Novo ◽  
Gc Content ◽  
Library Preparation ◽  
Sequencing Data ◽  
Bacterial Genomes ◽  
Dna Amount ◽  
High Quality ◽  
Preparation Methods

Sequencing bacterial genomes has traditionally required large amounts of genomic DNA (~1 μg). There have been few studies to determine the effects of the input DNA amount or library preparation method on the quality of sequencing data. Several new commercially available library preparation methods enable shotgun sequencing from as little as 1 ng of input DNA. In this study, we evaluated the NEBNext Ultra library preparation reagents for sequencing bacterial genomes. We have evaluated the utility of NEBNext Ultra for resequencing andde novoassembly of four bacterial genomes and compared its performance with the TruSeq library preparation kit. The NEBNext Ultra reagents enable high quality resequencing andde novoassembly of a variety of bacterial genomes when using 100 ng of input genomic DNA. For the two most challenging genomes (Burkholderiaspp.), which have the highest GC content and are the longest, we also show that the quality of both resequencing andde novoassembly is not decreased when only 10 ng of input genomic DNA is used.

scholarly journals Identification and Function of the pdxY Gene, Which Encodes a Novel Pyridoxal Kinase Involved in the Salvage Pathway of Pyridoxal 5′-Phosphate Biosynthesis in Escherichia coli K-12

1998 ◽  
Vol 180 (7) ◽  
pp. 1814-1821 ◽  
Author(s):  
Yong Yang ◽  
Ho-Ching Tiffany Tsui ◽  
Tsz-Kwong Man ◽  
Malcolm E. Winkler
Keyword(s):  
Escherichia Coli ◽  
De Novo ◽  
Specific Activity ◽  
Salvage Pathway ◽  
Pyridoxal Kinase ◽  
Triple Mutant ◽  
Reading Frame ◽  
E Coli ◽  
Specific Activities ◽  
K 12

ABSTRACT pdxK encodes a pyridoxine (PN)/pyridoxal (PL)/pyridoxamine (PM) kinase thought to function in the salvage pathway of pyridoxal 5′-phosphate (PLP) coenzyme biosynthesis. The observation that pdxK null mutants still contain PL kinase activity led to the hypothesis that Escherichia coli K-12 contains at least one other B6-vitamer kinase. Here we support this hypothesis by identifying the pdxY gene (formally, open reading frame f287b) at 36.92 min, which encodes a novel PL kinase. PdxY was first identified by its homology to PdxK in searches of the complete E. coli genome. Minimal clones of pdxY + overexpressed PL kinase specific activity about 10-fold. We inserted an omega cassette intopdxY and crossed the resultingpdxY::ΩKanr mutation into the bacterial chromosome of a pdxB mutant, in which de novo PLP biosynthesis is blocked. We then determined the growth characteristics and PL and PN kinase specific activities in extracts ofpdxK and pdxY single and double mutants. Significantly, the requirement of the pdxB pdxK pdxY triple mutant for PLP was not satisfied by PL and PN, and the triple mutant had negligible PL and PN kinase specific activities. Our combined results suggest that the PL kinase PdxY and the PN/PL/PM kinase PdxK are the only physiologically important B6vitamer kinases in E. coli and that their function is confined to the PLP salvage pathway. Last, we show thatpdxY is located downstream from pdxH (encoding PNP/PMP oxidase) and essential tyrS (encoding aminoacyl-tRNATyr synthetase) in a multifunctional operon.pdxY is completely cotranscribed with tyrS, but about 92% of tyrS transcripts terminate at a putative Rho-factor-dependent attenuator located in thetyrS-pdxY intercistronic region.

scholarly journals TrancriptomeReconstructoR, A Data-Driven Annotation of Complex Transcriptomes

2020 ◽  
Author(s):  
Maxim Ivanov ◽  
Albin Sandelin ◽  
Sebastian Marquardt
Keyword(s):  
De Novo ◽  
Gene Annotation ◽  
R Package ◽  
Sequence Information ◽  
Rna Seq ◽  
Sequencing Data ◽  
Gene Model ◽  
Preparation Methods ◽  
Downstream Analysis

Abstract Background: The quality of gene annotation determines the interpretation of results obtained in transcriptomic studies. The growing number of genome sequence information calls for experimental and computational pipelines for de novo transcriptome annotation. Ideally, gene and transcript models should be called from a limited set of key experimental data. Results: We developed TranscriptomeReconstructoR, an R package which implements a pipeline for automated transcriptome annotation. It relies on integrating features from independent and complementary datasets: i) full-length RNA-seq for detection of splicing patterns and ii) high-throughput 5' and 3' tag sequencing data for accurate definition of gene borders. The pipeline can also take a nascent RNA-seq dataset to supplement the called gene model with transient transcripts.We reconstructed de novo the transcriptional landscape of wild type Arabidopsis thaliana seedlings as a proof-of-principle. A comparison to the existing transcriptome annotations revealed that our gene model is more accurate and comprehensive than the two most commonly used community gene models, TAIR10 and Araport11. In particular, we identify thousands of transient transcripts missing from the existing annotations. Our new annotation promises to improve the quality of A.thaliana genome research.Conclusions: Our proof-of-concept data suggest a cost-efficient strategy for rapid and accurate annotation of complex eukaryotic transcriptomes. We combine the choice of library preparation methods and sequencing platforms with the dedicated computational pipeline implemented in the TranscriptomeReconstructoR package. The pipeline only requires prior knowledge on the reference genomic DNA sequence, but not the transcriptome. The package seamlessly integrates with Bioconductor packages for downstream analysis.

scholarly journals Optimizing de novo genome assembly from PCR-amplified metagenomes

2018 ◽  
Author(s):  
Simon Roux ◽  
Gareth Trubl ◽  
Danielle Goudeau ◽  
Nandita Nath ◽  
Estelle Couradeau ◽  
...  
Keyword(s):  
Genome Assembly ◽  
De Novo Assembly ◽  
De Novo ◽  
Pcr Amplification ◽  
Error Rates ◽  
De Novo Genome Assembly ◽  
Low Input ◽  
Assembly Algorithm ◽  
Coverage Bias ◽  
Assembly Pipeline

Background. Metagenomics has transformed our understanding of microbial diversity across ecosystems, with recent advances enabling de novo assembly of genomes from metagenomes. These metagenome-assembled genomes are critical to provide ecological, evolutionary, and metabolic context for all the microbes and viruses yet to be cultivated. Metagenomes can now be generated from nanogram to subnanogram amounts of DNA. However, these libraries require several rounds of PCR amplification before sequencing, and recent data suggest these typically yield smaller and more fragmented assemblies than regular metagenomes. Methods. Here we evaluate de novo assembly methods of 169 PCR-amplified metagenomes, including 25 for which an unamplified counterpart is available, to optimize specific assembly approaches for PCR-amplified libraries. We first evaluated coverage bias by mapping reads from PCR-amplified metagenomes onto reference contigs obtained from unamplified metagenomes of the same samples. Then, we compared different assembly pipelines in terms of assembly size (number of bp in contigs ≥ 10kb) and error rates to evaluate which are the best suited for PCR-amplified metagenomes. Results. Read mapping analyses revealed that the depth of coverage within individual genomes is significantly more uneven in PCR-amplified datasets versus unamplified metagenomes, with regions of high depth of coverage enriched in short inserts. This enrichment scales with the number of PCR cycles performed, and is presumably due to preferential amplification of short inserts. Standard assembly pipelines are confounded by this type of coverage unevenness, so we evaluated other assembly options to mitigate these issues. We found that a pipeline combining read deduplication and an assembly algorithm originally designed to recover genomes from libraries generated after whole genome amplification (single-cell SPAdes) frequently improved assembly of contigs ≥ 10kb by 10 to 100-fold for low input metagenomes. Conclusions. PCR-amplified metagenomes have enabled scientists to explore communities traditionally challenging to describe, including some with extremely low biomass or from which DNA is particularly difficult to extract. Here we show that a modified assembly pipeline can lead to an improved de novo genome assembly from PCR-amplified datasets, and enables a better genome recovery from low input metagenomes.

scholarly journals Fully phased human genome assembly without parental data using single-cell strand sequencing and long reads

2020 ◽  
Author(s):  
David Porubsky ◽  
◽  
Peter Ebert ◽  
Peter A. Audano ◽  
Mitchell R. Vollger ◽  
...  
Keyword(s):  
Single Cell ◽  
Genome Assembly ◽  
De Novo ◽  
Error Rates ◽  
Sequencing Data ◽  
Single Nucleotide Variants ◽  
De Novo Genome Assembly ◽  
Parental Data ◽  
Human Genome Assembly ◽  
Long Read

AbstractHuman genomes are typically assembled as consensus sequences that lack information on parental haplotypes. Here we describe a reference-free workflow for diploid de novo genome assembly that combines the chromosome-wide phasing and scaffolding capabilities of single-cell strand sequencing1,2 with continuous long-read or high-fidelity3 sequencing data. Employing this strategy, we produced a completely phased de novo genome assembly for each haplotype of an individual of Puerto Rican descent (HG00733) in the absence of parental data. The assemblies are accurate (quality value > 40) and highly contiguous (contig N50 > 23 Mbp) with low switch error rates (0.17%), providing fully phased single-nucleotide variants, indels and structural variants. A comparison of Oxford Nanopore Technologies and Pacific Biosciences phased assemblies identified 154 regions that are preferential sites of contig breaks, irrespective of sequencing technology or phasing algorithms.

Sign in / Sign up

Export Citation Format

Share Document