Filogenesi molecolare

1. Filogenesi molecolare

2. Geni ortologhi e geni paraloghi Geni ortologhi: geni simili riscontrabili in organismi correlati tra loro. Il fenomeno della speciazione porta alla divergenza dei geni e quindi delle proteine che essi codificano. es. l’ α-globina di uomo e di topo hanno iniziato a divergere circa 80 milioni di anni fa, quando avvenne la divisione che dette vita ai primati e ai roditori. I due geni sono da considerarsi ortologhi. Geni paraloghi: geni originati dalla duplicazione di un unico gene nello stesso organismo. es. α-globina e β-globina umana hanno iniziato a divergere in seguito alla duplicazione di un gene globinico ancestrale. I due geni sono da considerarsi paraloghi.

3. Gene ancestrale duplicazione genica Gene A Gene B speciazione ortologhi Gene A2 Gene A1 paraloghi Gene B1 Gene B2 ortologhi Specie 1 Specie 2

4. Distanze genetiche Per la distinzione filogenetica di due sequenze, è necessario conoscere quanto esse divergono. Serve quindi un parametro ripetibile, la distanza genetica. Ottengo così il Numero di Sostituzioni per Sito che può essere una sottostima della distanza vera perché le mutazioni possono anche essere REVERSIONI ed avvenire sullo stesso sito (sito con multiple hits). Numero di sostituzioni osservate nell’allineamento Distanza = --------------------------------------------------------------- Lunghezza complessiva dell’allineamento per le proteine si usaspesso d = - ln ( 1 - p - 0,2 p2) (kimura) d: distanza p: frazione di amino acidi differenti

5. Proteine o acidi nucleici? Sequenze proteiche: - necessitano di matrici si sostituzione 20x20, molto complesse da trattare.- sono espressione di sole regioni codificanti.- aminoacidi identici possono essere espressione di più codoni Sequenze nucleotidiche:- sono descrivibili con matrici 4x4.- possono essere estratte da sequenze genomiche non codificanti- non hanno degenerazione né ridondanza. => Per la filogenesi molecolare è preferibile utilizzare le sequenze nucleotidiche

6. Assunzioni a priori Per formulare un modello evolutivo è necessario considerare alcuni aspetti generali che possono essere considerati assunzioni a priori del modello: 1. tutti i siti evolvono in modo indipendente2. tutti i siti mutano con la stessa probabilità3. tutte le sostituzioni sono ugualmente probabili4. la velocità di sostituzione è costante nel tempo5. la composizione delle basi è costantein generale si può affermare che: maggiore è il numero di assunzioni a priori - maggiore è la semplicità del modello - minore è l’attendibilità dei risultati

7. Alcuni modelli proposti JC69 - Jukes & Cantor (1969). Assunzioni: 1,2,3,4,5. KIM - Kimura (1980). Assunzioni: 1,2,4,5. Tamura (1992). Assunzioni: 1,2,4. F81 - Felsenstein (1981). Assunzioni: 1,2,3,4. Hasegawa (1985). Assunzioni: 1,2,4. GTR – Lanave (1984). Assunzioni: 1,2. Il modello GTR è il più complesso possibile, compatibilmente con una trattazione matematica adeguata. Tutti i modelli richiedono però che la composizione in basi sia stazionaria, altrimenti non possono essere applicati. Bisogna quindi verificarlo prima.

8. L’orologio molecolare L’evoluzione è un processo inevitabilmente divergente e il numero di mutazioni che si accumulano nel tempo è direttamente proporzionale al tempo intercorso dalla divergenza delle sequenze in analisi. (1965, Zuckerkandl e Pauling). Se questo è vero, data una distanza genetica calcolata osservando le divergenze, è possibile ottenere il tempo trascorso dal momento in cui due sequenze hanno cominciato a divergere. Inoltre, se la velocità di accumulo delle mutazioni è costante, è possibile la datazione degli organismi in base a un solo dato verificato di distanza temporale. K K es. Distanza uomo topo: 80 milioni di anniV = ----- => t = ------ Distanza tra le alpha-globine: 0.093 2t 2V V = 0.56 x 10-8 sostituzioni/sito anno => anche se l’orologio molecolare è vero, non è universale, perché siti diversi hanno diversi tassi di mutazione.

9. Alberi filogenetici Sono grafi costituiti da NODI, che rappresentano le unità tassonomiche e da RAMI che uniscono i nodi, rappresentando le distanze tra i due. Si definisce TOPOLOGIA la struttura generale di un albero. Se ai rami non si dà valenza di distanza evolutiva, ho un CLADOGRAMMA, altrimenti ho un FILOGRAMMA. Alberi CON RADICE: accetta come vera l’ipotesi dell’orologio molecolare e i nodi stanno in un preciso ordine temporale. Alberi SENZA RADICE: non prevede significato evolutivo in termini temporali. F A nodo nodo B E radice ramo C C D B D ramo E F A

10. Daterminazione della topologia E’ il primo passo dell’analisi filogenetica.Il numero complessivo di alberi che si possono costruire con N sequenze (denominate OTU, cioè operational taxonomic units) è (2N - 3)!Nr = ---------------- 2N-3 (N - 3)! Rooted: Unrooted: (2N - 5)!Nu = ---------------- 2N-3 (N - 3)! es. 10 OTU = ((2*10)-3)! / 27(7)! = circa 35 milioni

11. Metodi per la crezione degli alberi Metodologie: Algoritmi di clusterizzazione (UPMGA, Neighbour Joining): si basano sull’osservazione delle distanze genetiche calcolare in precedenza. Algoritmi di ottimizzazione (Minima evoluzione): ottimizzazione degli alberi in base a criteri obiettivi di qualità. Origine dei dati: Sequenze omologhe multiallineate: tempi di calcolo molto maggiori Distanze genetiche pre-calcolate: tempi di calcolo minori,

12. UPGMA: (Unweighted Pair Group Method with Aritmetic mean) E’ il metodo più semplice. Assume l’orologio molecolare. Calcolata una matrice di distanze, vengono prese le OTU più simili, raggruppare e ricalcolata la matrice di distanze. La prima coppia viene clusterizzata con quella che ha distanza minore con esse. Il punto di mezzo nell’albero è il nodo tra i clusters, la lunghezza del ramo è la distanza che separa i vari elementi del cluster. Chimp 0.007 0.009 0.007 Human 0.018 0.018 Gorilla Owl monkey 0.027 0.030 0.027 Spider monkey

13. Neighbour-Joining Si parte da una topologia a stella. Tra OTU vicine introduce un ramo corrispondente alla distanza tra le due, e ripete l’inserimento di rami per tutte le coppie di OTU. Tra tutti gli alberi possibili viene via via scelto quello che ha il totale di lunghezza dei rami monore. B C n A 1 D C B C 2 B D F D A E F E A E F

14. Massima parsimonia Non si lavora più con le distanze ma con le sequenze: si cerca l’albero che richiede il minor numero di sostituzioni che spieghino le differenze osservate tra le sequenze in analisi. Non si lavora su tutti i siti del multiallineamneto, ma solo su alcuni SIGNIFICATIVI, cioè che presentano almeno 2 residui uguali. 1 2 3A ATGCATAGCAGCATGCATCGB ATGCATAGCAGCTTGTTTCGC ATCCATACCAGCATGTGTCGD ATCCATACCAGCTTGTTTCG G C G G A C A B G C C C C B D D G C C C G G 1 A B D C C C A G A T A A A T A C A B A C A B G C C G T B T T C A T D D B T G D D T T T A T T A T G G 2 3 A B A B D D C T C C A C

15. Massima verosimiglianza Per ciascun sito del multiallineamento si calcola la probabilità che esso sia generato da un preciso albero. Estendendo la ricerca di probabilità a tutti i possibili alberi, ottengo l’albero a maggior probabilità di rappresentare il multiallineamento. E’ il metodo più corretto di valutare le significatività degli alberi, ma ha un grosso limite: L’ enorme mole di calcolo per ottenere il risultato => È impossibile ottene risultati per multiallineamenti di più di 20-30 sequenze Sono stati prodotti vari metodi di approssimazione, tra cui uno detto quartet puzzle (QP) in cui si fanno operazioni su gruppi di 4 sequenze, e i risultati vengono intersecati tra loro.

16. Bootstrapping E’ un metodo relativamente semplice per determinare la significatività di un risultato derivante da una analisi complessa: Si tratta di un RICAMPIONAMENTO casuale di sequenze a cui vengono applicate le stesse procedure applicate alle sequenze vere. In questo modo ottengo due valori: - uno effettivo, quello delle sequenze in analisi - uno che rappresenta il valore che otterrei con un campione casuale. Dal loro confronto è possibile ottenere una stima della veridicità del risultato. Maggiore è il numero di operazione di bootstrapping, maggiore è l’attendibilità dei valori. Maggiore è la percentuale di risultati di bootstrapping che si accordano col risultato vero, maggiore è la precisione della stima del valore vero.

17. Principali programmi per l’analisi filogenetica PHYLIP: è un pacchetto di programmi curati da Joseph Felsenstein che non hanno interfaccia grafica. Sono programmi rapidi, efficienti e versatili, ma un po’ scomodi da utilizzare a causa della diversità dei formati di file utilizzati e dell’approccio di separazione delle varie fasi. PAUP: come Phylip, ma con un numero maggiore di test statistici per la validazione degli alberi. Ma si paga. Markov: effettua il calcolo delle distanze genetiche utilizzando il modello GTR e effettua il test di stazionarietà della composizione del dataset. E’ on-line, ma per utilizzarlo bisogna registrarsi.

18. Principali programmi per l’analisi filogenetica MrBayes: utilizza per il calcolo delle distanze un metodo introdotto recentemente di inferenza bayesiana, che sta dando ottimi risultati in termini di qualità degli alberi generati. PhyloWin: dotato di interfaccia grafica, è piuttosto semplice e consente analisi con i metodi presenti anche in Phylip e PAUP. MEGA: attualmente è il programma più completo e facile da usare possibile. Ha un’interfaccia grafica molto intuitiva, una serie di strumenti grafici per la comprensione dei riusultati e comprende numerosi approcci per il calcolo delle distanze, per la costruzione degli alberi e per per la valutazione dell’attendibilità con numerosi test statistici. Inoltre il manuale di Mega è molto completo, quasi un libro sulla filogenesi molecolare.