splitthreader:-a-graphical-algorithm-for- analysis...
TRANSCRIPT
![Page 1: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/1.jpg)
SplitThreader: A graphical algorithm for analysis of highly rearranged and amplified cancer genomes
Maria Nattestad
Mike Schatz lab
Cold Spring Harbor Laboratory and Johns Hopkins University
![Page 2: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/2.jpg)
Complex genome to study: SK-‐BR-‐3
(Davidson et al, 2000)
OHen used for pre-‐clinical research on Her2-‐targeJng therapeuJcs such as HercepJn (Trastuzumab) and resistance to these therapies.
Most commonly used Her2-‐amplified breast cancer cell line
80 chromosomes instead of 46
![Page 3: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/3.jpg)
Align to reference with BWA-‐MEM
Copy number analysis
Structural variant calling using Sniffles
SplitThreader graphical analysis
PacBio SMRT sequencing
Assembly
Assembly-‐based variant-‐calling
using AssemblyJcs
EvoluJon Fusion Visualizer
POSTER: Sniffles by Fritz Sedlazeck at poster
1310
![Page 4: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/4.jpg)
661 structural variants called by Sniffles
Her2 oncogene
Why SplitThreader EvoluJon?
![Page 5: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/5.jpg)
Her2
Her2
Chr 17: 83 Mb
8 Mb
![Page 6: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/6.jpg)
Her2
Her2
![Page 7: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/7.jpg)
Her2 GSDMB
TATDN1
8 Mb
RARA
PKIA
Chromosome 17
POSTER: ValidaJon of these SV calls using BioNano – Sara Goodwin
at poster 607
![Page 8: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/8.jpg)
Can we infer the order of events?
instead of organisms: DNA sequences
instead of traits: variants
1 2 3
Copy number:
Variants:
Let’s apply parsimony to a toy example
More parsimonious More likely true
Less parsimonious Less likely true
![Page 9: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/9.jpg)
Good: Most parsimonious Less parsimonious
variant 1 variant 3
variant 2
variant 1 variant 3
1 2 3
1
variant 1
variant 2
3
1 2 3
1 3
1 3
1 2 3
Variant 1 evolving twice is highly unlikely
1 2 3
![Page 10: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/10.jpg)
SplitThreader: Building the graph
CHR 1
CHR 2
ATCGCCTA
GTCCATAG
80
100
20
ATCG CCGA
ATAG GTCC
CHR 1
CHR 2100
20
80
start port end port
![Page 11: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/11.jpg)
SplitThreader EvoluJon
2 2
Source/
Sink
6 13
7
8 0 0 6
0
0 7
0 0
![Page 12: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/12.jpg)
Her2 GSDMB
TATDN1
RARA
PKIA
Chr 17 Chr 8
1. Healthy chromosome 17 2. TranslocaJon into
chromosome 8 3. TranslocaJon within
chromosome 8 4. Complex variant and
inverted duplicaJon within chromosome 8
5. TranslocaJon within chromosome 8
SplitThreader EvoluJon infers the complex history of this region for the first Jme
![Page 13: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/13.jpg)
• Found candidate fusions using IsoSeq (PacBio’s RNA-‐seq) – 3 million reads, mean 3.5 kb read length
• Use SplitThreader to look for DNA evidence – 16 gene fusions with a direct link in the DNA – 6 gene fusions each of which is created by threading through 2 variants
SplitThreader Fusion
![Page 14: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/14.jpg)
SplitThreader Fusion: CYTH1-‐EIF3H example
CYTH1
EIF3H
Chr 17
Chr 8
CYTH1 EIF3H
30 IsoSeq reads
3.7 Mb
?
![Page 15: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/15.jpg)
SplitThreader Fusion: CYTH1-‐EIF3H
CYTH1
EIF3H
Criteria Red path Green path
Sense-‐sense stranded No Yes
Fewer variants 2 2
More split read DNA evidence 40, 11 27, 36
Shorter transcript (threshold < 1Mb)
1.1 Mb 164 kb
![Page 16: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/16.jpg)
SplitThreader Fusion: CYTH1-‐EIF3H results
CYTH1
EIF3H
Chr 17
Chr 8
CYTH1 EIF3H
30 IsoSeq reads
8 kb
3.7 Mb
27 36 ?
![Page 17: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/17.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 18: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/18.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 19: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/19.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 20: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/20.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 21: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/21.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 22: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/22.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 23: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/23.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 24: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/24.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 25: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/25.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 26: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/26.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 27: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/27.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 28: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/28.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 29: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/29.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 30: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/30.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 31: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/31.jpg)
hip://qb.cshl.edu/spliihreader
SplitThreader Visualizer
![Page 32: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/32.jpg)
Explain oncogene amplificaJons
Trace gene fusions
qb.cshl.edu/splitthreader
SplitThreader EvoluJon
SplitThreader Fusion
SplitThreader Visualizer
SplitThreader
Explore the genome
![Page 33: SplitThreader:-A-graphical-algorithm-for- analysis …schatzlab.cshl.edu/presentations/2016/2016.02.12.AGBT...2016/02/12 · SplitThreader:-A-graphical-algorithm-for-analysis-of-highly-rearranged-and-amplified-cancer-genomes--Maria](https://reader033.vdocument.in/reader033/viewer/2022060414/5f12805495659c54530de2aa/html5/thumbnails/33.jpg)
Acknowledgments
Sara Goodwin Fritz Sedlazeck Philipp Rescheneder Marley Alford Timour Baslan Tyler Garvin Han Fang James Gurtowski Elizabeth Huion Melissa Kramer Eric Antoniou James Hicks W. Richard McCombie Michael Schatz
Karen Ng Timothy Beck Yogi Sundaravadanam John McPherson
Elizabeth Tseng Jason Chin