La varianza de dominancia en un solo locus

Estaba leyendo el libro "Genética y análisis de rasgos cuantitativos", de Lynch y Walsh. I cómo la covarianza entre dos individuos con EII Θ se divide solo en la varianza aditiva y el componente de varianza de dominancia, incluso en el simple 1 caso de locus.

Aquí mi comprensión del modelado (para el caso simple de un alelo):

Dado un valor genotípico GRAMO i , j de media 0 , i , j { 0 , 1 } encontramos numeros α 0 y α 1 minimizando los mínimos cuadrados de la siguiente forma mi ( GRAMO i , j α i α j ) 2 , donde la expectativa está sobre la población.

A continuación definimos los términos de error en cada caso como d i , j = GRAMO i , j α i α j . De las propiedades, vistas como funciones de la población α i es independiente de d i , j , y ambos tienen media 0 .

La afirmación que se hace en el libro es que dados dos individuos, con EII Θ y probabilidad de que el genotipo sea igual Δ , la covarianza de los genotipos GRAMO i , j y GRAMO k , yo es dado por,

cov ( GRAMO i , j , GRAMO k , yo ) = 2 Θ σ A 2 + Δ σ D 2 ,

dónde σ A 2 = Var ( α i ) , y σ D 2 = Var ( d i , j ) .

Expandiendo la LHS de la expresión, mostrando que mi [ ( α i + α j ) ( α k + α yo ) ] = Θ σ A 2 es bastante fácil También parece seguirse de que los términos mi [ ( α i + α j ) d k , yo ] = 0 de la independencia de los errores de la α .

Al analizar mi [ d i , j d k , yo ] , vemos que si ambos genotipos son iguales, lo que ocurre con probabilidad Δ , entonces esto se reduce a σ D 2 . Esto nos da un término Θ σ D 2 . Además, si ambos i , j y k , yo no son EII entonces la covarianza es 0 . Sin embargo, cuando uno de los dos alelos es IBD, entonces no me queda claro que esta covarianza seguirá siendo 0 .

El libro parece afirmar que, a menos que ambos alelos sean EII, d i , j y d k , yo son independientes No veo por qué este es el caso. ¿Me estoy perdiendo algo aquí? Agradecería cualquier ayuda con esto.

@rg255: Es el capítulo 7 en la página 143-144. Ellos manejan el caso de múltiples alelos, pero aquí lo estaba simplificando al caso de un solo alelo.

Respuestas (1)

Parece que pude resolver la pregunta. Mi entendimiento es el siguiente.

Dejar ( X 1 , Y 1 ) , ( X 2 , Y 2 ) (corresponde a i , j y k , yo en la pregunta) sean los genotipos de dos individuos en algún lugar. Suponemos que ninguno de los individuos es consanguíneo. En este caso definir,

Δ 7 = PR ( Ambos alelos  { X 1 , Y 1 }  y  { X 2 , Y 2 }  son EII ) , Δ 8 = PR ( Exactamente un par de alelos de los cuatro pares es IBD ) , Δ 9 = PR (  No hay EII entre los dos individuos. ) .
Δ 7 , Δ 8 y Δ 9 se calculan a partir de los pedigríes. Por ejemplo, para hermanos no gemelos, Δ 7 = 1 4 , Δ 8 = 1 2 , Δ 9 = 1 4 .

Observamos que el coeficiente IBD (también llamado coeficiente de parentesco) Θ es

Θ = 1 2 Δ 7 + 1 4 Δ 8 .
También el coeficiente de fraternidad, Δ = Δ 7 . Así tenemos que la covarianza entre el valor genotípico de dos individuos está dada por
mi [ GRAMO ( X 1 , Y 1 ) GRAMO ( X 2 , Y 2 ) ] = Δ 7 mi [ GRAMO ( X 1 , Y 1 ) GRAMO ( X 2 , Y 2 ) | ambas EII ] + Δ 8 mi [ GRAMO ( X 1 , Y 1 ) GRAMO ( X 2 , Y 2 ) | una EII ] + Δ 9 mi [ GRAMO ( X 1 , Y 1 ) GRAMO ( X 2 , Y 2 ) | sin EII ] .

Tenga en cuenta además que,

mi [ GRAMO ( X 1 , Y 1 ) GRAMO ( X 2 , Y 2 ) | ambas EII ] = mi [ GRAMO 2 ( X 1 , Y 1 ) ] , = σ A 2 + σ D 2 .

mi [ GRAMO ( X 1 , Y 1 ) GRAMO ( X 2 , Y 2 ) | sin EII ] = mi [ GRAMO ( X 1 , Y 1 ) ] mi [ GRAMO ( X 2 , Y 2 ) ] , = 0.

mi [ GRAMO ( X 1 , Y 1 ) GRAMO ( X 2 , Y 2 ) | una EII ] = mi [ ( α ( X ) + α ( Y 1 ) + d ( X , Y 1 ) ) ( α ( X ) + α ( Y 2 ) + d ( X , Y 2 ) ) ] , = mi [ α 2 ( X ) ] + 2 mi ( α ( X ) d ( X , Y 1 ) ) + mi [ d ( X , Y 1 ) d ( X , Y 2 ) ] , = 1 2 σ A 2 + 2 mi [ α ( X ) mi [ d ( X , Y 1 ) | X ] ] + mi [ mi [ d ( X , Y 1 ) | X ] mi [ d ( X , Y 2 ) | X ] ] , = 1 2 σ A 2 .

Así tenemos eso,

mi [ GRAMO ( X 1 , Y 1 ) GRAMO ( X 2 , Y 2 ) ] = ( 1 2 Δ 7 + 1 4 Δ 8 ) 2 σ A 2 + Δ 7 σ D 2 , = 2 Θ σ A 2 + Δ σ D 2 ,
que es como se afirma.