LOD score
LOD score LOD score
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
- Page 2 and 3: La probabilità di trovare una fami
- Page 4 and 5: AaBb (P) AaBb AaBb (P) aabb (P) Aab
- Page 6: Il logaritmo del rapporto tra quest
- Page 10 and 11: II III I 1 2 1 2 3 4 5 6 7 8 9 EcoR
- Page 12 and 13: Il LOD SCORE passo per passo • Sc
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