Es bastante difícil saber a priori qué
valores de LR se pueden obtener con cualquier par de hipótesis, con cualquier
conjunto de marcadores y con cualquier población. Pero a veces necesitamos
conocer esa información, ya sea antes de hacer los análisis en el Lab o después
de haberlos hecho.
Con Familias 3 podemos saberlo!! Y no sabéis lo útil que me
está resultando esta herramienta. Os explico aquí cómo funciona, y en un
comentario posterior veremos cómo hacerlo en Familias.
Imaginaros que tenemos un caso en el que sólo el hermano de
un niño está disponible y queremos saber si efectivamente son hermanos o si no
están relacionados. Los hermanos comparten:
– 1 alelo idéntico por descendencia en el 50% de los casos
– los 2 alelos idénticos por descendencia en el 25% de los
casos
– ningún alelo idéntico por descendencia en el 25% restante
Pues bien, Familias puede simular pares de
perfiles genéticos (teniendo en cuenta las frecuencias alélicas de nuestra
población, como ya vimos), de la siguiente forma:
– Pares de perfiles que cumplan los requisitos para ser
hermanos (siguiendo el dibujo anterior)
– Pares de perfiles que NO cumplan los requisitos para ser
hermanos
Tras las simulaciones, Familias calcula el LR para cada par
de perfiles (suponiendo las hipótesis definidas) y nos ofrece:
– El rango de valores de LR que se obtienen cuando los
perfiles pertenecen realmente a dos hermanos
– El rango de valores de LR que se obtienen cuando los
perfiles pertenecen a dos personas no relacionadas familiarmente
Y ¿qué utilidad tiene esto? Pues podéis llevar a cabo las
simulaciones:
a) Antes de hacer el genotipado en el Lab,
para saber qué podemos esperar de un caso, y si no nos convence pues solicitar
más familiares o ampliar nuestra batería de marcadores
b) Después de hacer el genotipado y de calcular
nuestro LR, para saber si el valor de LR obtenido en el caso que estamos
investigando es un valor esperado, es decir, que está dentro de los rangos de
LRs (distribuciones de LRs) resultantes de las simulaciones. Y si nuestro valor
no está en la distribución, pues algo tenemos mal… (las hipótesis no son las
correctas, las tasas de mutación no son las adecuadas, hemos usado frecuencias
de otra población, etc)
Yo uso mucho esta herramienta después de hacer el análisis en
el Lab, pues me gusta saber si el LR que he obtenido es un LR
«normal» con esas hipótesis, para ese número de marcadores, en esa
población…
Las gracias se las tenemos que dar a Daniel y a Thore…
Brillante!!!
Efectivamente John. Con nuestros electros no sabemos si los alelos son idénticos por descendencia o por estado. Pero sí es verdad que los hermanos comparten 1 alelo idéntico por descendencia en el 50% de los casos, ninguno en el 25%, etc…, aunque nosotros no podamos diferenciarlos. Esto es la teoría, lo que ocurre si dos personas son hermanos.
Si analizamos con secuenciación masiva podremos distinguir algo más, pero si aún así nos resultan dos alelos con la misma secuencia, no sabremos si son idénticos por descendencia o por estado.
Mil gracias por tu comentario!!
Abrazos
Hola, una pregunta. ¿Cómo defines un alelo idéntico por descendencia? ¿Cómo sabes que el alelo es idéntico por descendencia?, para eso no necesitaríamos el pedigrí para hallar un ancestro en común? y al no tenerlo simplemente si se tiene un único hermano no se puede hablar de alelo idéntico, o sí?
Gracias
Mil gracias por tu comentario Marta.
Pues no, nunca he probado eso, pero me temo que no afina tanto… Si tu LR no está dentro del rango esperado puede ser por muchos motivos, no sólo porque las frecuencias que estás usando no son las adecuadas. También pudiera ser que las hipótesis no son las correctas, es decir, que no tienen la relación de parentesco que reclaman, por ejemplo.
Abrazos!!
Sí, sí. Muy interesante. Y se me ocurre lo siguiente:
Quiere esto decir, que si desconociera el origen poblacional de los individuos analizados ¿podría utilizar esta herramienta para deducir la población? Osea, si mis resultados están dentro del rango de LR esperado al utilizar unas determinadas frecuencias pero «se salen» al utilizar las de otra población? Quizá en poblaciones muy diferentes??? Has probado alguna vez esto??