Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Correcting palindromes in long reads after whole-genome amplification.

ABSTRACT: BACKGROUND:Next-generation sequencing requires sufficient DNA to be available. If limited, whole-genome amplification is applied to generate additional amounts of DNA. Such amplification often results in many chimeric DNA fragments, in particular artificial palindromic sequences, which limit the usefulness of long sequencing reads. RESULTS:Here, we present Pacasus, a tool for correcting such errors. Two datasets show that it markedly improves read mapping and de novo assembly, yielding results similar to these that would be obtained with non-amplified DNA. CONCLUSIONS:With Pacasus long-read technologies become available for sequencing targets with very small amounts of DNA, such as single cells or even single chromosomes.

SUBMITTER: Warris S

PROVIDER: S-EPMC6218980 | biostudies-literature | 2018 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Correcting palindromes in long reads after whole-genome amplification.

Warris Sven S Schijlen Elio E van de Geest Henri H Vegesna Rahulsimham R Hesselink Thamara T Te Lintel Hekkert Bas B Sanchez Perez Gabino G Medvedev Paul P Makova Kateryna D KD de Ridder Dick D

BMC genomics 20181106 1

<h4>Background</h4>Next-generation sequencing requires sufficient DNA to be available. If limited, whole-genome amplification is applied to generate additional amounts of DNA. Such amplification often results in many chimeric DNA fragments, in particular artificial palindromic sequences, which limit the usefulness of long sequencing reads.<h4>Results</h4>Here, we present Pacasus, a tool for correcting such errors. Two datasets show that it markedly improves read mapping and de novo assembly, yie ...[more]

PMID: 30400848

Similar Datasets

Whole-genome haplotyping using long reads and statistical methods.

Project description:The rapid growth of sequencing technologies has greatly contributed to our understanding of human genetics. Yet, despite this growth, mainstream technologies have not been fully able to resolve the diploid nature of the human genome. Here we describe statistically aided, long-read haplotyping (SLRH), a rapid, accurate method that uses a statistical algorithm to take advantage of the partially phased information contained in long genomic fragments analyzed by short-read sequencing. For a human sample, as little as 30 Gbp of additional sequencing data are needed to phase genotypes identified by 50× coverage whole-genome sequencing. Using SLRH, we phase 99% of single-nucleotide variants in three human genomes into long haplotype blocks 0.2-1 Mbp in length. We apply our method to determine allele-specific methylation patterns in a human genome and identify hundreds of differentially methylated regions that were previously unknown. SLRH should facilitate population-scale haplotyping of human genomes.

| S-EPMC4073643 | biostudies-literature

VeChat: correcting errors in long reads using variation graphs.

Project description:Error correction is the canonical first step in long-read sequencing data analysis. Current self-correction methods, however, are affected by consensus sequence induced biases that mask true variants in haplotypes of lower frequency showing in mixed samples. Unlike consensus sequence templates, graph-based reference systems are not affected by such biases, so do not mistakenly mask true variants as errors. We present VeChat, as an approach to implement this idea: VeChat is based on variation graphs, as a popular type of data structure for pangenome reference systems. Extensive benchmarking experiments demonstrate that long reads corrected by VeChat contain 4 to 15 (Pacific Biosciences) and 1 to 10 times (Oxford Nanopore Technologies) less errors than when being corrected by state of the art approaches. Further, using VeChat prior to long-read assembly significantly improves the haplotype awareness of the assemblies. VeChat is an easy-to-use open-source tool and publicly available at https://github.com/HaploKit/vechat .

| S-EPMC9636371 | biostudies-literature

Whole genome complete resequencing of Bacillus subtilis natto by combining long reads with high-quality short reads.

Project description:De novo microbial genome sequencing reached a turning point with third-generation sequencing (TGS) platforms, and several microbial genomes have been improved by TGS long reads. Bacillus subtilis natto is closely related to the laboratory standard strain B. subtilis Marburg 168, and it has a function in the production of the traditional Japanese fermented food "natto." The B. subtilis natto BEST195 genome was previously sequenced with short reads, but it included some incomplete regions. We resequenced the BEST195 genome using a PacBio RS sequencer, and we successfully obtained a complete genome sequence from one scaffold without any gaps, and we also applied Illumina MiSeq short reads to enhance quality. Compared with the previous BEST195 draft genome and Marburg 168 genome, we found that incomplete regions in the previous genome sequence were attributed to GC-bias and repetitive sequences, and we also identified some novel genes that are found only in the new genome.

| S-EPMC4199626 | biostudies-literature

Comparing vertebrate whole-genome shotgun reads to the human genome.

Project description:Multi-species sequence comparisons are a very efficient way to reveal conserved genes. Because sequence finishing is expensive and time consuming, many genome sequences are likely to stay incomplete. A challenge is to use these fragmented data for understanding the human genome. Methods for using cross-species whole-genome shotgun sequence (WGS) for genome annotation are described in this paper. About one-half million high-quality rat WGS reads (covering 7.5% of the rat genome) generated at the Baylor College of Medicine Human Genome Sequencing Center were compared with the human genome. Using computer-generated random reads as a negative control, a set of parameters was determined for reliable interpretation of BLAST search results. About 10% of the rat reads contain regions that are conserved in the human genomic sequence and about one-third of these include known gene-coding regions. Mapping the conserved regions to human chromosomes showed a 23-fold enrichment for coding regions compared with noncoding regions. This approach can also be applied to other mammalian genomes for gene finding. These data predicted approximately 42,500 genes in the human, slightly more than reported previously.

| S-EPMC311156 | biostudies-literature

Primase-based whole genome amplification.

Project description:In vitro DNA amplification methods, such as polymerase chain reaction (PCR), rely on synthetic oligonucleotide primers for initiation of the reaction. In vivo, primers are synthesized on-template by DNA primase. The bacteriophage T7 gene 4 protein (gp4) has both primase and helicase activities. In this study, we report the development of a primase-based Whole Genome Amplification (pWGA) method, which utilizes gp4 primase to synthesize primers, eliminating the requirement of adding synthetic primers. Typical yield of pWGA from 1 ng to 10 ng of human genomic DNA input is in the microgram range, reaching over a thousand-fold amplification after 1 h of incubation at 37 degrees C. The amplification bias on human genomic DNA is 6.3-fold among 20 loci on different chromosomes. In addition to amplifying total genomic DNA, pWGA can also be used for detection and quantification of contaminant DNA in a sample when combined with a fluorescent reporter dye. When circular DNA is used as template in pWGA, 10(8)-fold of amplification is observed from as low as 100 copies of input. The high efficiency of pWGA in amplifying circular DNA makes it a potential tool in diagnosis and genotyping of circular human DNA viruses such as human papillomavirus (HPV).

| S-EPMC2490742 | biostudies-literature

Pacasus: Long reads and whole-genome amplication

Project description:Whole-genome amplification (WGA) create specific chimeric fragments, which consist mainly of palindrome sequences. We developed the tool Pacasus to detect and correct these palindromic sequences in long reads, for example from PacBio and Nanopore.

| PRJEB21791 | ENA

Minimizer-space de Bruijn graphs: Whole-genome assembly of long reads in minutes on a personal computer.

Project description:DNA sequencing data continue to progress toward longer reads with increasingly lower sequencing error rates. Here, we define an algorithmic approach, mdBG, that makes use of minimizer-space de Bruijn graphs to enable long-read genome assembly. mdBG achieves orders-of-magnitude improvement in both speed and memory usage over existing methods without compromising accuracy. A human genome is assembled in under 10 min using 8 cores and 10 GB RAM, and 60 Gbp of metagenome reads are assembled in 4 min using 1 GB RAM. In addition, we constructed a minimizer-space de Bruijn graph-based representation of 661,405 bacterial genomes, comprising 16 million nodes and 45 million edges, and successfully search it for anti-microbial resistance (AMR) genes in 12 min. We expect our advances to be essential to sequence analysis, given the rise of long-read sequencing in genomics, metagenomics, and pangenomics. Code for constructing mdBGs is freely available for download at https://github.com/ekimb/rust-mdbg/.

| S-EPMC8562525 | biostudies-literature

Automated reconstruction of whole-genome phylogenies from short-sequence reads.

Project description:Studies of microbial evolutionary dynamics are being transformed by the availability of affordable high-throughput sequencing technologies, which allow whole-genome sequencing of hundreds of related taxa in a single study. Reconstructing a phylogenetic tree of these taxa is generally a crucial step in any evolutionary analysis. Instead of constructing genome assemblies for all taxa, annotating these assemblies, and aligning orthologous genes, many recent studies 1) directly map raw sequencing reads to a single reference sequence, 2) extract single nucleotide polymorphisms (SNPs), and 3) infer the phylogenetic tree using maximum likelihood methods from the aligned SNP positions. However, here we show that, when using such methods to reconstruct phylogenies from sets of simulated sequences, both the exclusion of nonpolymorphic positions and the alignment to a single reference genome, introduce systematic biases and errors in phylogeny reconstruction. To address these problems, we developed a new method that combines alignments from mappings to multiple reference sequences and show that this successfully removes biases from the reconstructed phylogenies. We implemented this method as a web server named REALPHY (Reference sequence Alignment-based Phylogeny builder), which fully automates phylogenetic reconstruction from raw sequencing reads.

| S-EPMC3995342 | biostudies-literature

GenoFrag: software to design primers optimized for whole genome scanning by long-range PCR amplification.

Project description:Genome sequence data can be used to analyze genome plasticity by whole genome PCR scanning. Small sized chromosomes can indeed be fully amplified by long-range PCR with a set of primers designed using a reference strain and applied to several other strains. Analysis of the resulting patterns can reveal the genome plasticity. To facilitate such analysis, we have developed GenoFrag, a software package for the design of primers optimized for whole genome scanning by long-range PCR. GenoFrag was developed for the analysis of Staphylococcus aureus genome plasticity by whole genome amplification in approximately 10 kb-long fragments. A set of primers was generated from the genome sequence of S.aureus N315, employed here as a reference strain. Two subsets of primers were successfully used to amplify two portions of the N315 chromosome. This experimental validation demonstrates that GenoFrag is a robust and reliable tool for primer design and that whole genome PCR scanning can be envisaged for the analysis of genome diversity in S.aureus, one of the major public health concerns worldwide.

| S-EPMC373259 | biostudies-literature

Assessment of palindromes as platforms for DNA amplification in breast cancer.

Project description:DNA amplification, particularly of chromosomes 8 and 11, occurs frequently in breast cancer and is a key factor in tumorigenesis, often associated with poor prognosis. The mechanisms involved in the amplification of these regions are not fully understood. Studies from model systems have demonstrated that palindrome formation can be an early step in DNA amplification, most notably seen in the breakage-fusion-bridge (BFB) cycle. Therefore, palindromes might be associated with gene amplicons in breast cancer. To address this possibility, we coupled high-resolution palindrome profiling by the Genome-wide Analysis of Palindrome Formation (GAPF) assay with genome-wide copy-number analyses on a set of breast cancer cell lines and primary tumors to spatially associate palindromes and copy-number gains. We identified GAPF-positive regions distributed nonrandomly throughout cell line and tumor genomes, often in clusters, and associated with copy-number gains. Commonly amplified regions in breast cancer, chromosomes 8q and 11q, had GAPF-positive regions flanking and throughout the copy-number gains. We also identified amplification-associated GAPF-positive regions at similar locations in subsets of breast cancers with similar characteristics (e.g., ERBB2 amplification). These shared positive regions offer the potential to evaluate the utility of palindromes as prognostic markers, particularly in premalignant breast lesions. Our results implicate palindrome formation in the amplification of regions with key roles in breast tumorigenesis, particularly in subsets of breast cancers.

| S-EPMC3266031 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data