Reprever: Resolving low-copy duplicated sequences using template driven assembly

Sangwoo Kim, Paul Medvedev, Tara A. Paton, Vineet Bafna

Research output: Contribution to journalArticlepeer-review

5 Scopus citations

Abstract

Genomic sequence duplication is an important mechanism for genome evolution, often resulting in large sequence variations with implications for disease progression. Although paired-end sequencing technologies are commonly used for structural variation discovery, the discovery of novel duplicated sequences remains an unmet challenge. We analyze duplicons starting from identified high-copy number variants. Given paired-end mapped reads, and a candidate high-copy region, our tool, Reprever, identifies (a) the insertion breakpoints where the extra duplicons inserted into the donor genome and (b) the actual sequence of the duplicon. Reprever resolves ambiguous mapping signatures from existing homologs, repetitive elements and sequencing errors to identify breakpoint. At each breakpoint, Reprever reconstructs the inserted sequence using profile hidden Markov model (PHMM)-based guided assembly. In a test on 1000 artificial genomes with simulated duplication, Reprever could identify novel duplicates up to 97% of genomes within 3 bp positional and 1% sequence errors. Validation on 680 fosmid sequences identified and reconstructed eight duplicated sequences with high accuracy. We applied Reprever to reanalyzing a re-sequenced data set from the African individual NA18507 to identify >800 novel duplicates, including insertions in genes and insertions with additional variation. polymerase chain reaction followed by capillary sequencing validated both the insertion locations of the strongest predictions and their predicted sequence.

Original languageEnglish (US)
Pages (from-to)e128
JournalNucleic acids research
Volume41
Issue number12
DOIs
StatePublished - Jul 2013

All Science Journal Classification (ASJC) codes

  • Genetics

Fingerprint Dive into the research topics of 'Reprever: Resolving low-copy duplicated sequences using template driven assembly'. Together they form a unique fingerprint.

Cite this