Discovery of Genomic Structural Variations with Next-Generation Sequencing Data

Transcription

1 Discovery of Genomic Structural Variations with Next-Generation Sequencing Data Advanced Topics in Computational Genomics Slides from Marcel H. Schulz, Tobias Rausch (EMBL), and Kai Ye (Leiden University)

2 Computational Methods

3 Detecting Genomic Rearrangements Reference Mate-pair or paired-end mapping abnormalities Split-Read alignments Read depth signals courtesy of Tobias Rausch (EMBL)

4 Detecting Genomic Rearrangements Reference Unmapped or single-anchored reads Mate-pair or paired-end mapping abnormalities Split-Read alignments Read depth signals Local assembly courtesy of Tobias Rausch (EMBL)

5 courtesy of Tobias Rausch (EMBL)

6 courtesy of Tobias Rausch (EMBL)

7 Insertions Deletions courtesy of Tobias Rausch (EMBL)

8 Lee et al. (2009) courtesy of Tobias Rausch (EMBL) Korbel et al. (2007)

9 courtesy of Tobias Rausch (EMBL)

10 courtesy of Tobias Rausch (EMBL)

11 courtesy of Tobias Rausch (EMBL)

12 courtesy of Tobias Rausch (EMBL)

13 1 Copy 1 Copy 0 Copy 2 Copy 2 Copy courtesy of Tobias Rausch (EMBL) Chiang et al. (2009)

14 Down-Syndrom Partial Trisomie 21 courtesy of Tobias Rausch (EMBL) Xie et al. (2009)

15 Human cancer cell lines compared to normal cell lines (SeqSeq algorithm, no fixed window size, multiple change points method ) Chiang et al. (2009)

16 With reads of length bps are we able to find the exact breakpoint of a structural variation?

17 With reads of length bps are we able to find the exact breakpoint of a structural variation? Yes using split-read mapping Donor Reference Example for read of length 40: Expected random matches for a 12bp read-prefix in the human genome?

18 With reads of length bps are we able to find the exact breakpoint of a structural variation? Yes using split-read mapping Donor Reference Example for read of length 40: Expected random matches for a 12bp read-prefix in the human genome?

19 With reads of length bps are we able to find the exact breakpoint of a structural variation? Yes using anchored split-read mapping Donor Reference mappable read mate provides anchor to narrow down search space Medvedev et al. (2009)

20 The Pindel algorithm (Deletions) How to do that? Ye et al. (2009)

21 The Pindel algorithm (Deletions) Use 3 end of left read as anchor point Use pattern growth to search for minimum and maximum unique substrings from the 3 end of the unmapped read (<=2x insert size) Ye et al. (2009)

22 #&)-./!'0&12-./!(3!%0&&$).!/)45&2 ATGCA ATCAAGTATGCTTAGC!"!#$%&$'($)!*!++ +, courtesy of Kai Ye (Leiden U.)

23 #&)-./!'0&12-./!(3!%0&&$).!/)45&2 ATGCA ATCAAGTATGCTTAGC!"!#$%&$'($)!*!++ +, courtesy of Kai Ye (Leiden U.)

24 #&)-./!'0&12-./!(3!%0&&$).!/)45&2 ATGCA ATCAAGTATGCTTAGC!"!#$%&$'($)!*!++ +, courtesy of Kai Ye (Leiden U.)

25 #&)-./!'0&12-./!(3!%0&&$).!/)45&2 ATGCA ATCAAGTATGCTTAGC!"!#$%&$'($)!*!++ +, courtesy of Kai Ye (Leiden U.)

26 #&),-.!'/&01,-.!(2!%/&&$)-!.)34&1 ATGCA ATCAAGTATGCTTAGC 5,-,'6'!6-,76$!86(8&),-.9!:;< 5/=,'6'!6-,76$!86(8&),-.9!:;<>!"!#$%&$'($)!*!++ *! courtesy of Kai Ye (Leiden U.)

27 The Pindel algorithm (Deletions) Use 3 end of left read as anchor point Use pattern growth to search for minimum and maximum unique substrings from the 3 end of the unmapped read (<=2x insert size) Use pattern growth to search for minimum and maximum unique substrings from the 5 end of the unmapped read (read length + Max_D) starting from mapped end in step 2 Ye et al. (2009)

28 The Pindel algorithm (Deletions) Use 3 end of left read as anchor point Use pattern growth to search for minimum and maximum unique substrings from the 3 end of the unmapped read (<=2x insert size) Use pattern growth to search for minimum and maximum unique substrings from the 5 end of the unmapped read (read length + Max_D) starting from mapped end in step 2 check if complete unmapped read can be combined from 3 and 5 end substrings matches Ye et al. (2009)

29 The Pindel algorithm (Insertions) Use 3 end of left read as anchor point Use pattern growth to search for minimum and maximum unique substrings from the 3 end of the unmapped read (<=2x insert size) Use pattern growth to search for minimum and maximum unique substrings from the 5 end of the unmapped read (read length -1) starting from mapped end in step 2 check if complete unmapped read can be combined from 3 and 5 end substrings matches Ye et al. (2009)

30 The Pindel algorithm (Insertions) Use 3 end of left read as anchor point Use pattern growth to search for minimum and maximum unique substrings from the 3 end of the unmapped read (<=2x insert size) Use pattern growth to search for minimum and maximum unique substrings from the 5 end of the unmapped read (read length -1) starting from mapped end in step 2 check if complete unmapped read can be combined from 3 and 5 end substrings matches In initial Pindel version exact matches to reference where required Ye et al. (2009)

31 The Pindel algorithm (Real Data) Ye et al. (2009)

32 The Pindel algorithm (Real Data) Ye et al. (2009)

33 The Pindel algorithm for complex variants a) large deletion b) tandem duplication c) inversion d-f) same as a-c with non-template sequence (yellow part) Ye et al. Pindel manual

34 Acknowledgements Tobias Rausch (EMBL) Kai Ye (Leiden University Medical Center) Anne-Katrin Emde (Freie Universität Berlin) References Kai Ye, Marcel H. Schulz, Quan Long, Rolf Apweiler, and Zemin Ning Pindel: a pattern growth approach to detect break points of large deletions and medium sized insertions from paired-end short reads. Bioinformatics (2009) 25(21): Pindel homepage: SplazerS homepage: