scholarly journals Assembling Large Genomes with Single-Molecule Sequencing and Locality Sensitive Hashing

2014 ◽  
Author(s):  
Konstantin Berlin ◽  
Sergey Koren ◽  
Chen-Shan Chin ◽  
James Drake ◽  
Jane M Landolin ◽  
...  
Keyword(s):  
Single Molecule ◽  
De Novo ◽  
Locality Sensitive Hashing ◽  
Model Organisms ◽  
Smrt Sequencing ◽  
High Coverage ◽  
Celera Assembler ◽  
Single Molecule Sequencing ◽  
Long Reads ◽  
Long Read

We report reference-grade de novo assemblies of four model organisms and the human genome from single-molecule, real-time (SMRT) sequencing. Long-read SMRT sequencing is routinely used to finish microbial genomes, but the available assembly methods have not scaled well to larger genomes. Here we introduce the MinHash Alignment Process (MHAP) for efficient overlapping of noisy, long reads using probabilistic, locality-sensitive hashing. Together with Celera Assembler, MHAP was used to reconstruct the genomes of Escherichia coli, Saccharomyces cerevisiae, Arabidopsis thaliana, Drosophila melanogaster, and human from high-coverage SMRT sequencing. The resulting assemblies include fully resolved chromosome arms and close persistent gaps in these important reference genomes, including heterochromatic and telomeric transition sequences. For D. melanogaster, MHAP achieved a 600-fold speedup relative to prior methods and a cloud computing cost of a few hundred dollars. These results demonstrate that single-molecule sequencing alone can produce near-complete eukaryotic genomes at modest cost.

scholarly journals Single-molecule sequencing of theDrosophila serratagenome

2016 ◽  
Author(s):  
Scott L. Allen ◽  
Emily K. Delaney ◽  
Artyom Kopp ◽  
Stephen F. Chenoweth
Keyword(s):  
Single Molecule ◽  
De Novo ◽  
Sequence Data ◽  
Error Rates ◽  
Model Species ◽  
High Coverage ◽  
Single Molecule Sequencing ◽  
Long Read ◽  
Drosophila Serrata ◽  
High Degree

ABSTRACTLong read sequencing technology promises to greatly enhancede novoassembly of genomes for non-model species. While error rates have been a large stumbling block, sequencing at high coverage allows reads to be self-corrected. Here we sequence andde novoassemble the genome ofDrosophila serrata, a non-model species from themontiumsubgroup that has been well studied for clines and sexual selection. Using 11 PacBio SMRT cells, we generated 12 Gbp of raw sequence data comprising approximately 65x whole genome coverage. Read lengths averaged 8,940 bp (NRead50 12,200) with the longest read at 53 Kbp. We self-corrected reads using the PBDagCon algorithm and assembled the genome using the MHAP algorithm within the PBcR assembler. Total genome length was 198 Mbp with an N50 just under 1 Mbp. Contigs displayed a high degree of arm-level conservation withD. melanogaster. We also provide an initial annotation for this genome usingin silicogene predictions that were supported by RNA-seq data.

scholarly journals lra: the Long Read Aligner for Sequences and Contigs

2020 ◽  
Author(s):  
Jingwen Ren ◽  
Mark JP Chaisson
Keyword(s):  
Single Molecule ◽  
De Novo ◽  
Data Types ◽  
Single Molecule Sequencing ◽  
Detection Algorithms ◽  
Link Type ◽  
Long Reads ◽  
Oxford Nanopore ◽  
Long Read ◽  
Linear Cost

AbstractMotivationIt is computationally challenging to detect variation by aligning long reads from single-molecule sequencing (SMS) instruments, or megabase-scale contigs from SMS assemblies. One approach to efficiently align long sequences is sparse dynamic programming (SDP), where exact matches are found between the sequence and the genome, and optimal chains of matches are found representing a rough alignment. Sequence variation is more accurately modeled when alignments are scored with a gap penalty that is a convex function of the gap length. Because previous implementations of SDP used a linear-cost gap function that does not accurately model variation, and implementations of alignment that have a convex gap penalty are either inefficient or use heuristics, we developed a method, lra, that uses SDP with a convex-cost gap penalty. We use lra to align long-read sequences from PacBio and Oxford Nanopore (ONT) instruments as well as de novo assembly contigs.ResultsAcross all data types, the runtime of lra is between 52-168% of the state of the art aligner minimap2 when generating SAM alignment, and 9-15% of an alternative method, ngmlr. This alignment approach may be used to provide additional evidence of SV calls in PacBio datasets, and an increase in sensitivity and specificity on ONT data with current SV detection algorithms. The number of calls discovered using pbsv with lra alignments are within 98.3-98.6% of calls made from minimap2 alignments on the same data, and give a nominal 0.2-0.4% increase in F1 score by Truvari analysis. On ONT data with SV called using Sniffles, the number of calls made from lra alignments is 3% greater than minimap2-based calls, and 30% greater than ngmlr based calls, with a 4.6-5.5% increase in Truvari F1 score. When applied to calling variation from de novo assembly contigs, there is a 5.8% increase in SV calls compared to minimap2+paftools, with a 4.3% increase in Truvari F1 score.Availability and implementationAvailable in bioconda: https://anaconda.org/bioconda/lra and github: https://github.com/ChaissonLab/[email protected], [email protected]

scholarly journals A whole genome atlas of 81 Psilocybe genomes as a resource for psilocybin production.

F1000Research ◽  
2021 ◽  
Vol 10 ◽  
pp. 961
Author(s):  
Kevin McKernan ◽  
Liam Kane ◽  
Yvonne Helbert ◽  
Lei Zhang ◽  
Nathan Houde ◽  
...  
Keyword(s):  
Gene Cluster ◽  
Single Molecule ◽  
Reference Genome ◽  
De Novo ◽  
Genomic Diversity ◽  
Sequence Coverage ◽  
Single Molecule Sequencing ◽  
Contiguous Gene ◽  
Long Read ◽  
Interesting Variation

The Psilocybe genus is well known for the synthesis of valuable psychoactive compounds such as Psilocybin, Psilocin, Baeocystin and Aeruginascin. The ubiquity of Psilocybin synthesis in Psilocybe has been attributed to a horizontal gene transfer mechanism of a ~20Kb gene cluster. A recently published highly contiguous reference genome derived from long read single molecule sequencing has underscored interesting variation in this Psilocybin synthesis gene cluster. This reference genome has also enabled the shotgun sequencing of spores from many Psilocybe strains to better catalog the genomic diversity in the Psilocybin synthesis pathway. Here we present the de novo assembly of 81 Psilocybe genomes compared to the P.envy reference genome. Surprisingly, the genomes of Psilocybe galindoi, Psilocybe tampanensis and Psilocybe azurescens lack sequence coverage over the previously described Psilocybin synthesis pathway but do demonstrate amino acid sequence homology to a less contiguous gene cluster and may illuminate the previously proposed evolution of psilocybin synthesis.

scholarly journals SMRT sequencing reveals differential patterns of methylation in two O111:H- Shiga toxigenicEscherichia coliisolates from a historic hemolytic uremic syndrome outbreak in Australia

2017 ◽  
Author(s):  
Brian M. Forde ◽  
Lauren J. McAllister ◽  
James C. Paton ◽  
Adrienne W. Paton ◽  
Scott A. Beatson
Keyword(s):  
Single Molecule ◽  
Genomic Analysis ◽  
Smrt Sequencing ◽  
Promoter Regions ◽  
Food Borne ◽  
Long Reads ◽  
Long Read ◽  
Shiga Toxin 2 ◽  
Uremic Syndrome ◽  
Potential Promoter

AbstractShiga toxigenicEscherichia coli(STEC) are important food-borne pathogens and a major cause of haemorrhagic colitis and haemolytic-uremic syndrome (HUS) worldwide. In 1995 a severe HUS outbreak in Adelaide occurred. A recent genomic analysis of STEC O111:H-strains 95JB1 and 95NR1 from this outbreak found that the more virulent isolate, 95NR1, harboured two additional copies of the Shiga toxin 2 (Stx2) genes although the structure of the Stx2-converting prophages could not be fully resolved due to the fragmented assembly. In this study we have used Pacific Biosciences (PacBio) single molecule real-time (SMRT) long read sequencing to characterise the complete epigenome (genome and methylome) of 95JB1 and 95NR1. Using long reads we completely resolved the structure of two, tandemly inserted, stx2-converting phage in 95NR1. Our analysis of the methylome of 95NR1 and 95JB1 identified hemi-methylation of a novel motif (5’-CTGCm6AG-3’) in more than 4000 sites in the 95NR1 genome. These sites were entirely unmethalyted in the 95JB1, including at least 180 potential promoter regions that could explain regulatory differences between the strains. We identified a Type IIG methyltransferase encoded in both genomes in association with three additional genes in an operon-like arrangement. IS1203mediated disruption of this operon in 95JB1 is the likely cause of the observed differential patterns of methylation between 95NR1 and 95JB1. This study demonstrates the enormous potential of PacBio SMRT sequencing to resolve complex prophage regions and reveal the genetic and epigenetic heterogeneity within a clonal population of bacteria.

scholarly journals The evaluation of RNA-Seq de novo assembly by PacBio long read sequencing

2019 ◽  
Author(s):  
Yifan Yang ◽  
Michael Gribskov
Keyword(s):  
Real Time ◽  
De Novo ◽  
Critical Issue ◽  
Evaluation Methods ◽  
Model Organisms ◽  
Rna Seq ◽  
Long Reads ◽  
Long Read ◽  
Set Up ◽  
Downstream Analysis

AbstractRNA-Seq de novo assembly is an important method to generate transcriptomes for non-model organisms before any downstream analysis. Given many great de novo assembly methods developed by now, one critical issue is that there is no consensus on the evaluation of de novo assembly methods yet. Therefore, to set up a benchmark for evaluating the quality of de novo assemblies is very critical. Addressing this challenge will help us deepen the insights on the properties of different de novo assemblers and their evaluation methods, and provide hints on choosing the best assembly sets as transcriptomes of non-model organisms for the further functional analysis. In this article, we generate a “real time” transcriptome using PacBio long reads as a benchmark for evaluating five de novo assemblers and two model-based de novo assembly evaluation methods. By comparing the de novo assmblies generated by RNA-Seq short reads with the “real time” transcriptome from the same biological sample, we find that Trinity is best at the completeness by generating more assemblies than the alternative assemblers, but less continuous and having more misassemblies; Oases is best at the continuity and specificity, but less complete; The performance of SOAPdenovo-Trans, Trans-AByss and IDBA-Tran are in between of five assemblers. For evaluation methods, DETONATE leverages multiple aspects of the assembly set and ranks the assembly set with an average performance as the best, meanwhile the contig score can serve as a good metric to select assemblies with high completeness, specificity, continuity but not sensitive to misassemblies; TransRate contig score is useful for removing misassemblies, yet often the assemblies in the optimal set is too few to be used as a transcriptome.

scholarly journals AsmMix: A pipeline for high quality diploid de novo assembly

2021 ◽  
Author(s):  
Pei Wu ◽  
Chao Liu ◽  
Ou Wang ◽  
Xia Zhao ◽  
Fang Chen ◽  
...  
Keyword(s):  
Single Molecule ◽  
De Novo ◽  
Variant Calling ◽  
The Other ◽  
Second Step ◽  
Small Scale ◽  
Mixing Process ◽  
High Quality ◽  
Single Molecule Sequencing ◽  
Long Read

AbstractIn this paper, we report a pipeline, AsmMix, which is capable of producing both contiguous and high-quality diploid genomes. The pipeline consists of two steps. In the first step, two sets of assemblies are generated: one is based on co-barcoded reads, which are highly accurate and haplotype-resolved but contain many gaps, the other assembly is based on single-molecule sequencing reads, which is contiguous but error-prone. In the second step, those two sets of assemblies are compared and integrated into a haplotype-resolved assembly with fewer errors. We test our pipeline using a dataset of human genome NA24385, perform variant calling from those assemblies and then compare against GIAB Benchmark. We show that AsmMix pipeline could produce highly contiguous, accurate, and haplotype-resolved assemblies. Especially the assembly mixing process could effectively reduce small-scale errors in the long read assembly.

scholarly journals Resolving the Complexity of Human Skin Metagenomes Using Single-Molecule Sequencing

mBio ◽  
2016 ◽  
Vol 7 (1) ◽  
Author(s):  
Yu-Chih Tsai ◽  
Sean Conlan ◽  
Clayton Deming ◽  
Julia A. Segre ◽  
Heidi H. Kong ◽  
...  
Keyword(s):  
Microbial Community ◽  
Human Skin ◽  
Single Molecule ◽  
Smrt Sequencing ◽  
High Quality ◽  
Single Nucleotide ◽  
Single Molecule Sequencing ◽  
Short Read ◽  
Hybrid Approaches ◽  
Long Read

ABSTRACT Deep metagenomic shotgun sequencing has emerged as a powerful tool to interrogate composition and function of complex microbial communities. Computational approaches to assemble genome fragments have been demonstrated to be an effective tool for de novo reconstruction of genomes from these communities. However, the resultant “genomes” are typically fragmented and incomplete due to the limited ability of short-read sequence data to assemble complex or low-coverage regions. Here, we use single-molecule, real-time (SMRT) sequencing to reconstruct a high-quality, closed genome of a previously uncharacterized Corynebacterium simulans and its companion bacteriophage from a skin metagenomic sample. Considerable improvement in assembly quality occurs in hybrid approaches incorporating short-read data, with even relatively small amounts of long-read data being sufficient to improve metagenome reconstruction. Using short-read data to evaluate strain variation of this C. simulans in its skin community at single-nucleotide resolution, we observed a dominant C. simulans strain with moderate allelic heterozygosity throughout the population. We demonstrate the utility of SMRT sequencing and hybrid approaches in metagenome quantitation, reconstruction, and annotation. IMPORTANCE The species comprising a microbial community are often difficult to deconvolute due to technical limitations inherent to most short-read sequencing technologies. Here, we leverage new advances in sequencing technology, single-molecule sequencing, to significantly improve reconstruction of a complex human skin microbial community. With this long-read technology, we were able to reconstruct and annotate a closed, high-quality genome of a previously uncharacterized skin species. We demonstrate that hybrid approaches with short-read technology are sufficiently powerful to reconstruct even single-nucleotide polymorphism level variation of species in this a community.

scholarly journals LongStitch: High-quality genome assembly correction and scaffolding using long reads

2021 ◽  
Author(s):  
Lauren Coombe ◽  
Janet X Li ◽  
Theodora Lo ◽  
Johnathan Wong ◽  
Vladimir Nikolic ◽  
...  
Keyword(s):  
Genome Assembly ◽  
De Novo ◽  
Draft Genome ◽  
Model Organisms ◽  
High Quality ◽  
De Novo Genome Assembly ◽  
Long Reads ◽  
Long Read ◽  
Genomic Regions ◽  
Genome Assemblies

Background Generating high-quality de novo genome assemblies is foundational to the genomics study of model and non-model organisms. In recent years, long-read sequencing has greatly benefited genome assembly and scaffolding, a process by which assembled sequences are ordered and oriented through the use of long-range information. Long reads are better able to span repetitive genomic regions compared to short reads, and thus have tremendous utility for resolving problematic regions and helping generate more complete draft assemblies. Here, we present LongStitch, a scalable pipeline that corrects and scaffolds draft genome assemblies exclusively using long reads. Results LongStitch incorporates multiple tools developed by our group and runs in up to three stages, which includes initial assembly correction (Tigmint-long), followed by two incremental scaffolding stages (ntLink and ARKS-long). Tigmint-long and ARKS-long are misassembly correction and scaffolding utilities, respectively, previously developed for linked reads, that we adapted for long reads. Here, we describe the LongStitch pipeline and introduce our new long-read scaffolder, ntLink, which utilizes lightweight minimizer mappings to join contigs. LongStitch was tested on short and long-read assemblies of three different human individuals using corresponding nanopore long-read data, and improves the contiguity of each assembly from 2.0-fold up to 304.6-fold (as measured by NGA50 length). Furthermore, LongStitch generates more contiguous and correct assemblies compared to state-of-the-art long-read scaffolder LRScaf in most tests, and consistently runs in under five hours using less than 23GB of RAM. Conclusions Due to its effectiveness and efficiency in improving draft assemblies using long reads, we expect LongStitch to benefit a wide variety of de novo genome assembly projects. The LongStitch pipeline is freely available at https://github.com/bcgsc/longstitch.

scholarly journals Genetic Findings of Sanger and Nanopore Single-Molecule Sequencing in Patients with X-Linked Hearing Loss and Incomplete Partition Type III

2021 ◽  
Author(s):  
Ying Chen ◽  
Jiajun Qiu ◽  
Yingwei Wu ◽  
Huan Jia ◽  
Yi Jiang ◽  
...  
Keyword(s):  
Single Molecule ◽  
Hearing Aids ◽  
Sanger Sequencing ◽  
Clinical Characteristics ◽  
De Novo ◽  
Genetic Mutations ◽  
Single Molecule Sequencing ◽  
Long Read ◽  
Incomplete Partition ◽  
The Mean

Abstract BackgroundPOU3F4 is the causative gene for X-linked deafness-2 (DFNX2), characterized by incomplete partition type III (IP-III) malformation of the inner ear. The aim of this study was to investigate the clinical characteristics and molecular findings by Sanger or Nanopore single-molecule sequencing in IP-III patients. MethodsDiagnosis of IP-III was mainly based on clinical characteristics including radiological and audiological findings. Sanger sequencing of POU3F4 were carried out for these IP-III patients. For those patients with negative results for POU3F4 Sanger sequencing, Nanopore long-read single-molecule sequencing was used to identify the possible pathogenic variants. Hearing intervention outcomes of hearing aids fitting and cochlear implantation were also analyzed. Grouped by different locations of POU3F4 variants, aided PTA was further compared between patients in whom the variants located in the exon region or in the upstream region.ResultsIn total, 18 male patients from 14 unrelated families were diagnosed with IP-III. 10 variants were identified in POU3F4 by Sanger sequencing and 9 of these were novel (p.Val321Gly, p.Gln181*, p.Cys233*, p.Val215Gly, p.Arg282Gln, p.Trp57*, p.Gln316*, c.903_912 delins TGCCA and p.Arg205del). Four different deletions (DELs) that varied from 80 to 486 kb were identified 876-1503 kb upstream of POU3F4 by Nanopore long-read single-molecule sequencing. Of them, de novo genetic mutations occurred in 21.4% (3/14) of patients with POU3F4 mutations. Of these 18 patients, 7 had bilateral hearing aids (HAs) and 10 patients received unilateral cochlear implantation (CI). The mean aided pure tone average (PTA) for HAs and CI users were 41.1±5.18 and 40.3±7.59 dB HL respectively. The mean PTAs for whom the variants located in the exon and upstream regions were 39.6±6.31 vs 43.0±7.10 dB HL, which presented no significant difference (p=0.342).ConclusionsAmong IP-III patients, 28.6% (4/14) had no definite mutation in exon region of POU3F4, however, possible pathogenic deletions were identified in upstream region of this gen. De novo genetic mutations occurred in 21.4% (3/14) of patients with POU3F4 mutation. Hearing intervention outcomes of IP-III patients presented no difference regardless of the variants locations on exon or upstream regions.

scholarly journals Highly-accurate long-read sequencing improves variant detection and assembly of a human genome

2019 ◽  
Author(s):  
Aaron M. Wenger ◽  
Paul Peluso ◽  
William J. Rowell ◽  
Pi-Chuan Chang ◽  
Richard J. Hall ◽  
...  
Keyword(s):  
Single Molecule ◽  
De Novo ◽  
Structural Variants ◽  
Short Reads ◽  
Sequencing Technologies ◽  
Long Reads ◽  
Long Read ◽  
Variant Detection ◽  
High Quality Genome ◽  
Circular Consensus Sequencing

AbstractThe major DNA sequencing technologies in use today produce either highly-accurate short reads or noisy long reads. We developed a protocol based on single-molecule, circular consensus sequencing (CCS) to generate highly-accurate (99.8%) long reads averaging 13.5 kb and applied it to sequence the well-characterized human HG002/NA24385. We optimized existing tools to comprehensively detect variants, achieving precision and recall above 99.91% for SNVs, 95.98% for indels, and 95.99% for structural variants. We estimate that 2,434 discordances are correctable mistakes in the high-quality Genome in a Bottle benchmark. Nearly all (99.64%) variants are phased into haplotypes, which further improves variant detection. De novo assembly produces a highly contiguous and accurate genome with contig N50 above 15 Mb and concordance of 99.998%. CCS reads match short reads for small variant detection, while enabling structural variant detection and de novo assembly at similar contiguity and markedly higher concordance than noisy long reads.

Sign in / Sign up

Export Citation Format

Share Document