Dataset Information

Effective machine-learning assembly for next-generation amplicon sequencing with very low coverage.

ABSTRACT: BACKGROUND:In short-read DNA sequencing experiments, the read coverage is a key parameter to successfully assemble the reads and reconstruct the sequence of the input DNA. When coverage is very low, the original sequence reconstruction from the reads can be difficult because of the occurrence of uncovered gaps. Reference guided assembly can then improve these assemblies. However, when the available reference is phylogenetically distant from the sequencing reads, the mapping rate of the reads can be extremely low. Some recent improvements in read mapping approaches aim at modifying the reference according to the reads dynamically. Such approaches can significantly improve the alignment rate of the reads onto distant references but the processing of insertions and deletions remains challenging. RESULTS:Here, we introduce a new algorithm to update the reference sequence according to previously aligned reads. Substitutions, insertions and deletions are performed in the reference sequence dynamically. We evaluate this approach to assemble a western-grey kangaroo mitochondrial amplicon. Our results show that more reads can be aligned and that this method produces assemblies of length comparable to the truth while limiting error rate when classic approaches fail to recover the correct length. Finally, we discuss how the core algorithm of this method could be improved and combined with other approaches to analyse larger genomic sequences. CONCLUSIONS:We introduced an algorithm to perform dynamic alignment of reads on a distant reference. We showed that such approach can improve the reconstruction of an amplicon compared to classically used bioinformatic pipelines. Although not portable to genomic scale in the current form, we suggested several improvements to be investigated to make this method more flexible and allow dynamic alignment to be used for large genome assemblies.

SUBMITTER: Ranjard L

PROVIDER: S-EPMC6907241 | biostudies-literature | 2019 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Effective machine-learning assembly for next-generation amplicon sequencing with very low coverage.

Ranjard Louis L Wong Thomas K F TKF Rodrigo Allen G AG

BMC bioinformatics 20191211 1

<h4>Background</h4>In short-read DNA sequencing experiments, the read coverage is a key parameter to successfully assemble the reads and reconstruct the sequence of the input DNA. When coverage is very low, the original sequence reconstruction from the reads can be difficult because of the occurrence of uncovered gaps. Reference guided assembly can then improve these assemblies. However, when the available reference is phylogenetically distant from the sequencing reads, the mapping rate of the r ...[more]

PMID: 31829137

Similar Datasets

Project description:Molecular technologies have allowed laboratories to detect and establish the profiles of human cancers by identifying a variety of somatic variants. In order to improve personalized patient care, we have established a next-generation sequencing (NGS) test to screen for somatic variants in primary or advanced cancers. In this study, we describe the laboratory quality management program for NGS testing, and also provide an overview of the somatic variants identified in over 1000 patient samples as well as their implications in clinical practice.Over the past one-and-a-half years, our laboratory received a total of 1028 formalin-fixed, paraffin-embedded (FFPE) tumor tissues, which consisted of non-small-cell lung carcinomas (NSCLCs), colon adenocarcinomas, glioma/glioblastomas, melanomas, breast carcinomas, and other tumor types. During this time period, we implemented a series of quality control (QC) checks that included (1) pre-DNA extraction, (2) DNA quantification, (3) DNA quality, (4) library quantification, (5) post-emulsification PCR, and (6) post-sequencing metrics. At least 10 ng of genomic DNA (gDNA) were used to prepare barcoded libraries using the AmpliSeq CHPv2. Samples were multiplexed and sequenced on Ion Torrent 318 chips using the Ion PGM System. Variants were identified using the Variant Caller Plugin, and annotation and functional predictions were performed using the Golden Helix SVS.A total of 1005 samples passed QC1-3, and following additional library preparation QC checkpoints, 877 samples were sequenced. Samples were classified into two categories: wild-type (127) and positive for somatic variants (750). Somatic variants were classified into clinically actionable (60%) and non-actionable (40%).The use of NGS in routine clinical laboratory practice allowed for the detection of tumor profiles that are essential for the selection of targeted therapies and identification of applicable clinical trials, contributing to the improvement of personalized patient care in oncology.

Project description:BackgroundNext-generation sequencing allows genome-wide analysis of changes in chromatin states and gene expression. Data analysis of these increasingly used methods either requires multiple analysis steps, or extensive computational time. We sought to develop a tool for rapid quantification of sequencing peaks from diverse experimental sources and an efficient method to produce coverage tracks for accurate visualization that can be intuitively displayed and interpreted by experimentalists with minimal bioinformatics background. We demonstrate its strength and usability by integrating data from several types of sequencing approaches.ResultsWe have developed BAMscale, a one-step tool that processes a wide set of sequencing datasets. To demonstrate the usefulness of BAMscale, we analyzed multiple sequencing datasets from chromatin immunoprecipitation sequencing data (ChIP-seq), chromatin state change data (assay for transposase-accessible chromatin using sequencing: ATAC-seq, DNA double-strand break mapping sequencing: END-seq), DNA replication data (Okazaki fragments sequencing: OK-seq, nascent-strand sequencing: NS-seq, single-cell replication timing sequencing: scRepli-seq) and RNA-seq data. The outputs consist of raw and normalized peak scores (multiple normalizations) in text format and scaled bigWig coverage tracks that are directly accessible to data visualization programs. BAMScale also includes a visualization module facilitating direct, on-demand quantitative peak comparisons that can be used by experimentalists. Our tool can effectively analyze large sequencing datasets (~?100 Gb size) in minutes, outperforming currently available tools.ConclusionsBAMscale accurately quantifies and normalizes identified peaks directly from BAM files, and creates coverage tracks for visualization in genome browsers. BAMScale can be implemented for a wide set of methods for calculating coverage tracks, including ChIP-seq and ATAC-seq, as well as methods that currently require specialized, separate tools for analyses, such as splice-aware RNA-seq, END-seq and OK-seq for which no dedicated software is available. BAMscale is freely available on github (https://github.com/ncbi/BAMscale).

Dataset Information

Effective machine-learning assembly for next-generation amplicon sequencing with very low coverage.

Publications

Effective machine-learning assembly for next-generation amplicon sequencing with very low coverage.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets