LOD score

LOD score LOD score

24.05.2013 Views

Ad ogni evento di fecondazione, la probabilità che lo zigote sia maschio oppure femmina è ½. p = probabilità che sia maschio, q = probabilità che sia femmina. La probabilità di ottenere questo risultato per puro caso ESATTAMENTE in quest’ordine (prima figlia femmina, 5 maschi di fila e ultima figlia femmina) è p 5 q 2 cioè (1/2) 5 x (1/2) 2 cioè 0,07813

Ad ogni evento di fecondazione, la probabilità che lo zigote sia maschio oppure<br />

femmina è ½. p = probabilità che sia maschio, q = probabilità che sia femmina.<br />

La probabilità di ottenere questo risultato per puro caso ESATTAMENTE in<br />

quest’ordine (prima figlia femmina, 5 maschi di fila e ultima figlia femmina) è<br />

p 5 q 2<br />

cioè (1/2) 5 x (1/2) 2 cioè 0,07813


La probabilità di trovare una famiglia di 7 figli con 5 maschi e 2 femmine<br />

IN QUALSIASI ORDINE è<br />

N!<br />

K!*(N-K)!<br />

* p 5 *q 2<br />

Dove N è il numero totale di figli (7) e K il numero di maschi (5)<br />

Quindi<br />

5040<br />

120*2<br />

21 p 5 * q 2 = 0,164<br />

= 21


Supponiamo di avere un soggetto eterozigote per due loci indipendenti<br />

A a<br />

B b<br />

Secondo la legge dell’indipendenza, ci aspettiamo che questo soggetto produca 4 tipi<br />

di gameti nelle stesse proporzioni (probabilità = 0,25 ciascuno). Se incrociamo questo<br />

eterozigote con un omozigote doppio recessivo (reincrocio o test-cross), ci aspettiamo<br />

nei figli proporzioni che rispettino le proporzioni tra i gameti.<br />

A<br />

B<br />

A<br />

b<br />

25% 25% 25% 25%<br />

Gameti prodotti dall’eterozigote (loci indipendenti)<br />

a<br />

B<br />

a<br />

b


AaBb<br />

(P)<br />

AaBb<br />

AaBb<br />

(P)<br />

aabb<br />

(P)<br />

Aabb<br />

(R)<br />

aabb<br />

AaBb<br />

(P)<br />

aabb<br />

(P)<br />

aaBb<br />

(R)<br />

In quest’albero genealogico troviamo valori che si discostano dal risultato atteso sulla<br />

base dell’ipotesi dell’indipendenza. Qual’è la probabilità di ottenere questo risultato<br />

semplicemente per caso? Per convenzione, definiamo AB e ab le combinazioni parentali,<br />

mentre Ab e aB le combinazioni ricombinanti.<br />

Analizziamo i figli: per 3 volte A è stato trasferito insieme a B e per 2 volte a insieme<br />

a b (5 parentali)<br />

per 2 volte A con b e B con a (2 ricombinanti).<br />

Secondo l’ipotesi dell’indipendenza, ci aspettiamo che i 4 fenotipi siano presenti<br />

in eguali proporzioni e quindi che la probabilità di ottenere per caso il risultato<br />

osservato sia 0,25 7 , cioè 6,103 * 10 -5


Vediamo ora di calcolare la probabilità di ottenere i risultati osservati sulla base di una<br />

effettiva associazione desunta dall’albero genealogico. Osserviamo che vi sono 5<br />

parentali e 2 ricombinanti. 2 ricombinanti su 7 = 0,286. Supponiamo che i loci in<br />

questione siano quindi associati e che distino tra loro 28,6 cM. Dal momento che non<br />

sappiamo se nell’eterozigote gli alleli sono in concordanza o in opposizione, calcoliamo<br />

le probabilità come segue:<br />

la probabilità q di avere un certo ricombinante è 0,286/2 = 0,143.<br />

La probabilità p di avere un certo parentale è quindi<br />

(1 – 0,286)/2 = 0,714/2 = 0,357<br />

La probabilità di ottenere dati come quelli osservati sulla base del linkage è<br />

½ x (0,357) 5 x (0,143) 2 + ½ x (0,357) 2 x (0,143) 5 perché non si conosce la fase<br />

(cioè non è noto se l’eterozigote lo è in cis o in trans)<br />

risultato = 0,001362 cioè 1,36 x 10 -3


Il logaritmo del rapporto tra questo risultato (ipotesi di linkage) e quello ottenuto in<br />

precedenza (ipotesi di indipendenza) viene definito <strong>LOD</strong> <strong>score</strong><br />

Log 10<br />

1,36 x 10 -3<br />

6,103 x 10 -5<br />

= 1,348<br />

Il <strong>LOD</strong> <strong>score</strong> è statisticamente significativo quando<br />

si avvicina almeno a 3 ( probabilità di solo 1 su 1.000<br />

che il risultato sia stato ottenuto per caso)


+ -<br />

SE/se se/se<br />

+ - - - - + -<br />

SE/se se/se se/se se/se se/se SE/se se/se<br />

+ - - + + + -<br />

SE/se se/se se/se SE/se SE/se SE/se SE/se<br />

Linkage del gene della distrofia miotonica con il locus secretore


II<br />

III<br />

I<br />

1 2<br />

1 2 3 4 5 6 7 8 9<br />

EcoRI EcoRI EcoRI<br />

3 kb 1.7 kb<br />

sonda<br />

1 2 3 4 5 6 7 8 9<br />

*<br />

* *<br />

Analisi del linkage della NF1 con una sonda del cromosoma 17<br />

4.7 kb<br />

3 kb


Il <strong>LOD</strong> Score è un parametro statistico usato per valutare la significatività dei risultati di una<br />

analisi di linkage. In pratica, è la probabilità di ottenere un determinato risultato per 2 loci che<br />

siano effettivamente in linkage rispetto alla probabilità che il risultato sia dovuto al caso.<br />

probabilità che i loci siano in linkage a θ = 0,117 (1/2 fase 1 + 1/2 fase 2)<br />

<strong>LOD</strong> = log 10 -----------------------------------------------------------------------------------------probabilità<br />

che i loci non siano in linkage<br />

1/2 (0,441 15 x 0,058 2 ) + 1/2 (0,058 15 x 0,441 2 )<br />

<strong>LOD</strong> = log 10 --------------------------------------------------------- = 2,127<br />

0,25 17<br />

Valori di <strong>LOD</strong> <strong>score</strong> elevati indicano una elevata probabilità di linkage tra i loci presi in<br />

considerazione.


Il <strong>LOD</strong> SCORE passo per passo<br />

• Scegliere un valore di θ (frequenza di ricombinazione) basandosi sui dati del pedigree.<br />

(ad esempio se ci sono 2 ric. apparentì su 17, la frequenza è 2/17, cioè 0,117)<br />

• Determinare il numero di possibili parentali e di possibili ricombinanti<br />

(15P e 2R)<br />

• Determinare la probabilità di generare un parentale o un ricombinante<br />

Probabilità di ottenere 1 ricombinante: θ/2 = 0,117/2 = 0,058<br />

Probabilità di ottenere 1 parentale: (1- θ)/2 = (1- 0,117)/2 = 0,441<br />

• Determinare la probabilità totale di ottenere i dati osservati sulla base di θ.<br />

-La probabilità totale è il prodotto delle singole probabilità<br />

-Se l’esempio è costituito da 15P e 2R, allora<br />

p = ½ (0,441) 15 x (0,058) 2 + ½ (0,058) 15 x (0,441) 2 = 7,8 x 10 -9


Determinare la probabilità di ottenere i dati osservati in assenza di linkage.<br />

In questo caso, p = (0,25) 17 = 5,8 x 10 -11<br />

in quanto i 4 genotipi sono ugualmente probabili<br />

Determinare il <strong>LOD</strong> Score<br />

<strong>LOD</strong> = Z = log 10 (7,8 x 10 -9 / 5,8 x 10 -11 ) = 2,127

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!