Are you sure you want to leave this community? Leaving the community will revoke any permissions you have been granted in this community.
Structural variation is an important class of genetic variation in mammals. High-throughput sequencing (HTS) technologies promise to revolutionize copy-number variation (CNV) detection but present substantial analytic challenges. Converging evidence suggests that multiple types of CNV-informative data (e.g. read-depth, read-pair, split-read) need be considered, and that sophisticated methods are needed for more accurate CNV detection. We observed that various sources of experimental biases in HTS confound read-depth estimation, and note that bias correction has not been adequately addressed by existing methods. We present a novel read-depth-based method, GENSENG, which uses a hidden Markov model and negative binomial regression framework to identify regions of discrete copy-number changes while simultaneously accounting for the effects of multiple confounders. Based on extensive calibration using multiple HTS data sets, we conclude that our method outperforms existing read-depth-based CNV detection algorithms. The concept of simultaneous bias correction and CNV detection can serve as a basis for combining read-depth with other types of information such as read-pair or split-read in a single analysis. A user-friendly and computationally efficient implementation of our method is freely available.
Pubmed ID: 23275535
Publication data is provided by the National Library of Medicine ® and PubMed ®. Data is retrieved from PubMed ® on a weekly schedule. For terms and conditions see the National Library of Medicine Terms and Conditions.
Web search tool to find regions of similarity between biological sequences. Program compares nucleotide or protein sequences to sequence databases and calculates statistical significance. Used for identifying homologous sequences.
View all literature mentionsJava toolset for working with next generation sequencing data in the BAM format.
View all literature mentionsWeb tool to convert genome coordinates and genome annotation files between assemblies. Used to translate genomic coordinates from one assembly version into another and retrieves putative orthologous regions in other species using UCSC chained and netted alignments.
View all literature mentionsMus musculus with name C3H/HeJ from IMSR.
View all literature mentionsMus musculus with name CBA/J from IMSR.
View all literature mentionsCell line GM12878 is a Transformed cell line with a species of origin Homo sapiens (Human)
View all literature mentionsMus musculus with name DBA/2J from IMSR.
View all literature mentionsMus musculus with name AKR/J from IMSR.
View all literature mentionsMus musculus with name BALB/cJ from IMSR.
View all literature mentions