Detecting copy number variation with mated short reads.
The development of high-throughput sequencing (HTS) technologies has opened the door to novel methods for detecting copy number variants (CNVs) in the human genome. While in the past CNVs have been detected based on array CGH data, recent studies have shown that depth-of-coverage information from HTS technologies can also be used for the reliable identification of large copy-variable regions. Such methods, however, are hindered by sequencing biases that lead certain regions of the genome to be over- or undersampled, lowering their resolution and ability to accurately identify the exact breakpoints of the variants. In this work, we develop a method for CNV detection that supplements the depth-of-coverage with paired-end mapping information, where mate pairs mapping discordantly to the reference serve to indicate the presence of variation. Our algorithm, called CNVer, combines this information within a unified computational framework called the donor graph, allowing us to better mitigate the sequencing biases that cause uneven local coverage and accurately predict CNVs. We use CNVer to detect 4879 CNVs in the recently described genome of a Yoruban individual. Most of the calls (77%) coincide with previously known variants within the Database of Genomic Variants, while 81% of deletion copy number variants previously known for this individual coincide with one of our loss calls. Furthermore, we demonstrate that CNVer can reconstruct the absolute copy counts of segments of the donor genome and evaluate the feasibility of using CNVer with low coverage datasets.
Medvedev
Paul
P
Department of Computer Science, University of Toronto, Toronto, Ontario M5R 3G4, Canada.
Fiume
Marc
M
Dzamba
Misko
M
Smith
Tim
T
Brudno
Michael
M
D000465:Algorithms; D020029:Base Pairing / Q000502:physiology*; D001483:Base Sequence / Q000502:physiology; D019457:Chromosome Breakage; D002874:Chromosome Mapping / Q000379:methods*; D056915:DNA Copy Number Variations / Q000235:genetics; D004252:DNA Mutational Analysis / Q000379:methods*; D032621:DNA Shuffling; D059014:High-Throughput Nucleotide Sequencing / Q000379:methods*; D006801:Humans; D015203:Reproducibility of Results
Nat Methods. 2009 Jul;6(7):473-4
19483690
Nat Methods. 2009 Jan;6(1):99-103
19043412
Genome Biol. 2003;4(8):R50
12914658
Genome Res. 2009 Sep;19(9):1586-92
19657104
Bioinformatics. 2003 Oct;19 Suppl 2:ii162-71
14534186
J Comput Biol. 2009 Aug;16(8):1101-16
19645596
Nature. 2008 Nov 6;456(7218):53-9
18987734
Nat Genet. 2004 Sep;36(9):949-51
15286789
Bioinformatics. 2010 Aug 15;26(16):1938-44
20562449
Nature. 2006 Nov 23;444(7118):444-54
17122850
Genome Res. 2009 Sep;19(9):1527-41
19546169
Nat Genet. 2008 Oct;40(10):1166-74
18776908
Nat Genet. 2008 Jun;40(6):722-9
18438408
Nat Methods. 2009 Sep;6(9):677-81
19668202
Nat Genet. 2008 Oct;40(10):1199-203
18776910
Nat Genet. 2005 Jul;37(7):727-32
15895083
Genome Biol. 2009;10(3):R25
19261174
Genome Res. 2009 Jul;19(7):1270-8
19447966
Bioinformatics. 2010 Feb 15;26(4):565-7
20022973
Nucleic Acids Res. 2008 Sep;36(16):e105
18660515
Science. 2004 Jul 23;305(5683):525-8
15273396
Nat Methods. 2009 Nov;6(11 Suppl):S13-20
19844226
Proc Natl Acad Sci U S A. 2003 Jun 24;100(13):7696-701
12788976
Bioinformatics. 2009 Nov 1;25(21):2872-7
19528083
Nature. 2008 May 1;453(7191):56-64
18451855
Bioinformatics. 2008 Jul 1;24(13):i59-67
18586745
Genome Res. 2004 Nov;14(11):2336-46
15520295
Genome Res. 2008 Nov;18(11):1814-28
18849524
Nat Genet. 2007 Jul;39(7 Suppl):S16-21
17597776
Genome Biol. 2009;10(3):R32
19327155
Nat Genet. 2007 Jul;39(7 Suppl):S37-42
17597780
Science. 2007 Oct 19;318(5849):420-6
17901297
Nat Genet. 2009 Oct;41(10):1061-7
19718026
Genome Res. 2004 Sep;14(9):1786-96
15342561