corso di Genomica a a 2010 2011 lezione

  • Slides: 16
Download presentation
corso di Genomica a. a. 2010 -2011 lezione 35 -36 • laurea magistrale Biotecnologia

corso di Genomica a. a. 2010 -2011 lezione 35 -36 • laurea magistrale Biotecnologia Industriale martedì 25 Gennaio 2011 aula 6 A orario : Martedì ore 14. 00 - 16. 00 Giovedì ore 13. 00 - 15. 00 Esami: 24 Febbraio, 3 Marzo, 23 Marzo Lezioni fino al 10 Febbraio D. Frezza

Chromosome contact probability Calcolo della probabilità di contatti intercromosomici In(s) di coppie di loci

Chromosome contact probability Calcolo della probabilità di contatti intercromosomici In(s) di coppie di loci separati dalla distanza s sul cromosoma n Il valore decresce in maniera monotonica su ogni cromosoma con comportamento dei polimeri: - la distanza 3 D aumenta con la distanza sul cromosoma - valori in accordo con analisi 3 C e FISH (fluorescence in situ hybridization) - a distanze maggiori di 200 Mb la probabilità media di un’interazione è sempre maggiore della probabilità di interazione tra due cromosomi diversi: questa evidenza implica l’esistenza di territori cromosomici

Inter chromosome contact probability In figura 2 B si evidenzia che tra cromosomi più

Inter chromosome contact probability In figura 2 B si evidenzia che tra cromosomi più piccoli a partire dal 16 al 22 aumenta la frequenza di interazioni - in accordo con risultati FISH questi cromosomi spesso colocalizzano al centro del nucleo. -è particolare il chrms 18 che è piccolo ma povero di geni non interagisce spesso con gli altri chrms piccoli in accordo con la FISH che ne osserva la localizzazione frequente in periferia del nucleo

Zoom on individual chromosomes Matrice della frequenza dei contatti genome wide i cui valori

Zoom on individual chromosomes Matrice della frequenza dei contatti genome wide i cui valori sono il numero di prodotti di ligasi della library tra il locus “i” e “j” (SOM) - per vedere se esistono contatti più frequenti - normalizzazione della matrice M* dividendo i valori (entry) per il valore medio della probabilità di contatto per i loci ad una certa distanza genomica SOM - la matrice normalizzata mostra molti grandi blocchi di interazioni arricchite o impoverite che generano una forma a “plaid” a scacchi - se due loci (in questo caso regioni di 1 Mb) sono vicine nello spazio, si pensa che condividano i bordi ed abbiano condivisione di profili di interazione

matrice di correlazione C Definizione della matrice di correlazione C di ci-j Pearson tra

matrice di correlazione C Definizione della matrice di correlazione C di ci-j Pearson tra la riga ith e la colonna jth Questa trasformazione rende molto più a fuoco la matrice a “pattern” a scacchi (fig 3 C) - il 71% dei valori rappresenta una significatività con p<0. 05

Classificazione dei cromosomi Arbitrariamente parte A e B se i loci sono contatti arricchiti

Classificazione dei cromosomi Arbitrariamente parte A e B se i loci sono contatti arricchiti o impoveriti - analisi della componente principale (PC) fig S 1 - questa analisi coincide col pattern a “plaid” della matrice C - marcatura simile su cromosomi diversi hanno comportamento simile di interazioni - marcature dissimili su chrms diversi hanno comportamento dissimile - conclusione: l’intero genome può essere suddiviso in due tipi di compartimenti spaziali all’interno di ognuno dei quali si hanno interazioni più frequenti che non tra compartimenti diversi = territori cromosomici

Controllo sperimentale 3 D-FISH su chrms 14 Questa prova sperimentale di controllo per verificare

Controllo sperimentale 3 D-FISH su chrms 14 Questa prova sperimentale di controllo per verificare la predittività della matrice e del metodo HI-C su chrms 14 - presi i domini L 1, L 3 = compart A, L 2, L 4 = compart B - posizione sul chrms = L 1 - L 2 - L 3 - L 4 - prove 3 D FISH mostrano : L 3 -L 1 > L 2 - L 2 -L 4 > L 3 - su chrms 22 risultati analoghi Alta correlazione tra risultati 3 D-FISH ed HI-C Può servire per misurare le distanze

Comportamenti dei compartimenti A e B Loci del cmprtm B hanno ad una certa

Comportamenti dei compartimenti A e B Loci del cmprtm B hanno ad una certa distanza interazioni più frequenti del cmprtm A (fig S 4) - inferenza: il compartimento B è impacchettato più densamento del cmprtm A - evidenze FISH: loci cmprtm B hanno localizzazioni spaziali più vicine

Genoma compatto e rilassato Plot su assi log-log - comportamento medio della probabilità di

Genoma compatto e rilassato Plot su assi log-log - comportamento medio della probabilità di contatti intracromosomici come funzione della distanza sulla distribuzione dell’intero genoma I(s) - questo valore è più significativo nell’intervallo tra 500 kb e 7 Mb nella scalla di probabilità di contatto s-1 - l’intervallo è quello della grandezza dei domini genomici di cromatina aperta e chiusa

Comportamento tipo polimero (polymer like) Regioni cromosomiche a modello di equilibrio globulare Originalmente descritti

Comportamento tipo polimero (polymer like) Regioni cromosomiche a modello di equilibrio globulare Originalmente descritti come polimeri a modello annodato in un solvente scarso all’equilibrio (poca solubilità) Distinzione tra “modello globulare” e modello globulare all’equilibrio Modello alternativo di Grosberg : i polimeri come DNA interfasico si può organizzare come conformazione a lunga durata non all’equilibrio descritto come “globulo frattale” - stato condensato è formato da un polimero non intricato (districabile) quando si condensa in una serie di piccoli globuli in una configurazione “palline su una fibra”

modello a palline condensate Le palline servono come monomeri in cicli successivi di condensazione

modello a palline condensate Le palline servono come monomeri in cicli successivi di condensazione spontanea finchè rimane solo un singolo globulo di globuli La struttura risultante somiglia ad una curva di Peano Un frattale a traettoria continua che riempie uno spazio tridimensionale senza incrociarsi (intrigarsiannodarsi)

Modello a frattali per chrms Il modello a globuli frattali per spiegare strutture di

Modello a frattali per chrms Il modello a globuli frattali per spiegare strutture di regioni di cromatina è attraente: - manca di nodi e facilita il ripiegamento e spiegamento durante le attività genetiche di attivazione genica, repressione e ciclo cellulare - nel globulo frattale le regioni contigue tendono a formare settori la cui grandezza corrisponde alle grandezze originali delle regioni (fig 4 C) - in contrasto un globulo all’equilibrio è altamente annodato e manca di questi settori; invece le posizioni lineari e spaziali non sono correlate anche dopo poche megabasi - il globulo frattale non è mai stato osservato prima

conclusioni su scala di alcune Mb Analisi con simulazione Monte Carlo di 500 globuli

conclusioni su scala di alcune Mb Analisi con simulazione Monte Carlo di 500 globuli frattali all’equilibrio Possibilità di aumentare il numero di “reads” n 2 del sequenziamento per definire e mappare interazioni specifiche a lunga distanza tra enhancers, silenziatori ed isolatori

La struttura risultante somiglia ad una curva di Peano http: //library. wolfram. com/images/infocenter/Peano. Curve

La struttura risultante somiglia ad una curva di Peano http: //library. wolfram. com/images/infocenter/Peano. Curve 794. gif

from wide genome to genome Dalla descrizione globale delle interazioni alle strutture che interagiscono

from wide genome to genome Dalla descrizione globale delle interazioni alle strutture che interagiscono Quali sono e come sono organizzate sul genoma L’esempio del cluster della catena pesante delle Ig La scoperta degli enhancers, l’attivazione a livello epigenetico, l’organizzazione e le differenze tra uomo-topo Come è evoluto il cluster da roditori a primati

cluster Ig topo-uomo (duplicazione) TOR VERGATA centromero 12 F 1 BAC 199 M 11(AF

cluster Ig topo-uomo (duplicazione) TOR VERGATA centromero 12 F 1 BAC 199 M 11(AF 450245) J D V 14 q 32. 33 Ig. H 3’RR-1 g 4 g 2 g hs 4 hs 1, 2 hs 3 mta 1 3’a 2 g elk 2. 1 m J DV Em (cloni instabili) a 2 e mta 1 crip 2 crip 1 hole Human Igh cluster Chromosome 14 g 3 d a 1 e g 1 g 3 d m J DV Em hs 4 hs 1, 2 hs 3 J D V g 1 trascrizione hs 4 hs 3 B hs 1, 2 hs 3 A mta 1 g 2 a g 2 b a e mta 1 crip 2 crip 1 hole Mouse Igh cluster Chromosome 12 telomero 3’a 1 hs 4 hs 1. 2 hs 3 20 bp a 1 3’a 1 86, 035, 000 Ig. H 3’RR-2 86, 040, 000 86, 055, 000 elk 2. 2 hole 85, 894, 500 86, 045, 000 85, 904, 500 85, 914, 500 85, 924, 500 86, 060, 000 86, 065, 000 86, 070, 000 hs 4 hs 1. 2 85, 934, 500 85, 944, 500 86, 075, 000 hs 3 20 bp a 2 85, 954, 500 3’a 2