Complete(Genome(analysis:( Structural(varia5on(detec5on( ( - PowerPoint PPT Presentation

Complete(Genome(analysis:( Structural(varia5on(detec5on( ( SVDetect(tutorial( Alban&LERMINE&–&Ins/tut&Curie,&INSERM&U900,&Mines&ParisTech&

Workshop(outlines( 1) SVDetect(tool(presenta5on( ( (Methodology(&(func5onali5es( 2)(SVDetect(tutorial(on(Galaxy( ( (Hands(on(workshop(

SVDetect( SVdetect(is(a(tool(dedicated(to(detec5on(and(predic5on(of(inter/intra(chromosomal( rearrangements(from(pairedFend/mateFpair(sequencing(data( ( This(tool(detect(structural(varia5on((SV's)(by(using(sliding(windows(and(clustering( strategies,(and(also(allows(you(to(visualize(them(at(genome(scale(

Different(classes(of(SV’s( Lot(of(Tumors(have(chromosomal(muta5ons(in(their(genome.(This(can(affect( the(number(of(chromosomes,(and(also(its(structure,(implica5ng(large(events( such(as(dele5ons,(inser5ons,(duplica5ons,(inverisons,transloca5ons,(etc.( Alkan&et&al&2012&

SV's(discovery(from(sequencing(data( 4(different(methods:( ( F(Based(on(coverage(depth(differences((ControlF(FREEC)( F(SpliVed(reads(method,(local(reads(realignment(around(breaking(point( (Pindel)( F(De(novo(assembly,(new(sequences(F(inser5ons(F( F(PairedFend/mateFpair(sequences(constraints((Breakdancer,(PEMer,( SVDetect)( ( Reviews(on(thema:( Computa/onal&methods&for&discovering&structural&varia/on&with&nextKgenera/on& sequencing ( Paul(Medvedev(et(al,(Nature(Methods(6,(S13(F(S20((2009) ( ( ( Genome&structural&varia/on&discovery&and&genotyping.&& Alkan(C(et(al,(Nat(Rev(Genet.((2012)((

SVDetect(>(strategy( PairedFend/mateFpair(based(methods,(such(as(SVDetect,(evaluate(reads(distance,( order,(and(orienta5on(a^er(alignment(on(a(reference(genome.(Abnormal(or( discordants(pairs(are(clustered(to(define(their(type(of(SV ( Reference ALIGNMENT5OF5PAIRED5ENDS PairedFends reads INTRA&CHROMOSOMAL-REARRANGEMENTS- Abnormally mapped paired in5the5same chromosome CONSTRAINTS Ordering Distance5<5cutoff Distance5>5cutoff Orientation 1 2 Reference Mapping Sample Distance LARGE5INSERTION5? LARGE5DELETION5? Normally mapped pairs INTER&CHROMOSOMAL-REARRANGEMENTS- Abnormally5mapped5paired5in5two5different5chromosomes Reference CHR5A Reference Mapping Ex:(SOAPdenovo( Mapping Sample CHR5B NO5SV Zeitouni&et&al&2010& TRANSLOCATION5?

SVDetect(>(implementa?on( SVDETECT' 4'different'funcBons & COMPARE( BAM' preprocessing' LINKS2CIRCOS( FILTERING( LINKING( LINKS2BED( DETECTION( Reference'genome' SVDETECT' ConfiguraBon'file' Genome(frac?onned((5kb(fragments)( (perl)' ' Abnormal'pairs' C'DetecBon'parameters' 'BAM'file'' mandatory'or'opBonal' Chromosomal(links(crea?on(between(each(fragments(pairs( C'Filtering'parameters' supported(by(paired(reads(>(clusters( C'VisualizaBon'opBons' C'Parallelized'computaBon & Values(for(median(reads(distance,(order(and( orienta?on( FILTERING( Cluster(filtering:(chromosomes,(number(of(pairs,( CIRCOS( distance/order/orienta?on(coherence( COMPARE( Circos(ou(Bed(output( BED( SV(poten?al(type(defini?on( BAL/UNBAL(transloca?on( SV's(comparison(for(2(samples:(common,(specific(

Tes?ng(data( (Biological(data(coming(from(neuroblastoma(tumor( (Two(samples:(Tumoral(vs(Control((Cons5tu5onal(DNA)( (( (Illumina(GAII(mateFpair(sequencing,(50x50(bp,(insert(size(of(3kb( (( (Alignments(files(in(binary(format:(BAM( (( (SVDetect(need(2(types(of(files(as(input:( (( ( (F(Bam(file(containing(all(apparied(and(aligned(reads( ( (>(«(sample_mates.bam(»(&(«(reference_mates.bam(»( ( (( (F(LEN(file(giving(reference(chromosome(size( ( (>(«(hs18_chr5_chr11.len(»( ( ( ( 1 ! chr5 ! ! 180857966 ! ! 2 ! chr11 ! ! 134452384 !

Let’s(start(hands(on(SVDetect(tutorial!(