Tutorials › Statistical Genetics › Worked Pedigree From Diagram to Kinship and Covariance

Pedigrees and Genetic Relatedness · Tutorial 79 of 989

Worked Pedigree From Diagram to Kinship and Covariance

Trace a multigeneration pedigree systematically, then use its kinship and inbreeding coefficients to calculate additive genetic covariances.

Beginner 9 min read

What You'll Learn

  • State founder and inheritance assumptions before calculating relatedness
  • Trace shared-ancestor paths and calculate kinship coefficients
  • Find an offspring’s inbreeding coefficient from its parents’ kinship
  • Distinguish kinship coefficients from additive relationship coefficients
  • Convert a pedigree relationship matrix into additive genetic covariances
  • Interpret genetic variances and covariances in squared trait units

From Pedigree to Covariance

A pedigree is more than a family diagram: under stated assumptions, it gives the expected sharing of alleles identical by descent between individuals. Those expectations can be organized as kinship coefficients, converted to additive relationship coefficients, and then scaled by an additive genetic variance to obtain genetic covariances. A reliable calculation therefore begins with the pedigree’s assumptions and proceeds in a fixed order.

The worked pedigree below includes full siblings whose descendants mate, creating an inbred individual, and then a parent–offspring mating that creates a further inbred descendant. It is hypothetical. The calculations concern expected genetic relationships under the pedigree model, not observed genotypes at a particular variant.

Definition: The kinship coefficient \(\phi_{ij}\) is the probability that one allele sampled at random from individual \(i\) and one sampled at random from individual \(j\) are identical by descent. For an individual’s two alleles, the inbreeding coefficient \(F_i\) is their probability of being identical by descent. The additive relationship coefficient is \(A_{ij}=2\phi_{ij}\) for distinct individuals, and \(A_{ii}=1+F_i\).

The pedigree calculations assume autosomal inheritance, correct parentage, and Mendelian transmission: each offspring receives one allele at a locus from each parent, with either parental allele having probability one-half. Unless specified otherwise, founders are assumed unrelated to one another and noninbred. Here, “unrelated” means no identity by descent within the pedigree being modeled; it is a modeling assumption, not a claim that the founders could not share any distant ancestors in reality. The additive covariance calculations additionally assume a common additive genetic variance, denoted \(\sigma_A^2\), and include no shared environmental covariance.

The Hypothetical Pedigree

There are five founders: \(A\), \(B\), \(E\), \(G\), and \(K\). They are assumed mutually unrelated and noninbred. The other individuals and their parents are as follows.

IndividualParentsRelationship relevant to the calculation
CA and BFull sibling of D
DA and BFull sibling of C
FC and EDescendant of C
HD and GDescendant of D
IF and HDescendant of first cousins F and H
LI and KChild of I and an unrelated founder
MI and LChild of a parent–offspring mating

In generation order, \(A\) and \(B\) have children \(C\) and \(D\). Individual \(C\) has a child \(F\) with founder \(E\), while \(D\) has a child \(H\) with founder \(G\). The first cousins \(F\) and \(H\) have child \(I\). Individual \(I\) has child \(L\) with founder \(K\), and then \(I\) and \(L\) have child \(M\). This ordering makes the assumptions visible and ensures each nonfounder’s parents are identified before the individual’s coefficients are calculated.

Recursion: For distinct individuals \(i\) and \(j\), if \(i\) has parents \(s_i\) and \(d_i\), then \(\phi_{ij}=(\phi_{s_i j}+\phi_{d_i j})/2\). For an individual \(i\), \(\phi_{ii}=(1+F_i)/2\), and if \(i\) is a nonfounder with parents \(s_i,d_i\), then \(F_i=\phi_{s_i d_i}\).

The off-diagonal recursion follows because a randomly sampled allele from \(i\) came from either parent with probability one-half. Conditional on its parental source, the chance of identity by descent with an allele sampled from \(j\) is the corresponding parental kinship. For the diagonal, two alleles sampled from one individual are the same homolog with probability one-half; otherwise they are the individual’s two distinct alleles, which are identical by descent with probability \(F_i\). Thus \(\phi_{ii}=(1+F_i)/2\). These rules also make clear why an inbreeding coefficient is not itself a kinship coefficient between two different people.

Worked Example: Kinship of the First Cousins

Worked Example: Kinship of the First Cousins

Calculate \(\phi_{FH}\), the kinship of \(F\) and \(H\), using their shared ancestors. First verify the pedigree features that determine the calculation.

Individuals \(C\) and \(D\) are full siblings because they have the same two parents, \(A\) and \(B\). Their children \(F\) and \(H\) are therefore first cousins. Both \(F\) and \(H\) are noninbred: \(F\)’s parents \(C\) and \(E\) are unrelated under the founder assumptions, and \(H\)’s parents \(D\) and \(G\) are likewise unrelated. Each of \(A\) and \(B\) is a common ancestor of \(F\) and \(H\).

For either common ancestor, the path from \(F\) to that ancestor has two parent–child steps, and so does the path from \(H\). Because the common ancestors are noninbred and there are no additional shared ancestral paths, each contributes \((1/2)^{2+2+1}=1/32\). Add the contributions from \(A\) and \(B\):

$$ \phi_{FH}=\frac{1}{32}+\frac{1}{32}=\frac{1}{16}=0.0625. $$

The result is a kinship coefficient, not an additive relationship coefficient. For these distinct individuals, the additive relationship is \(A_{FH}=2\phi_{FH}=1/8=0.125\). This is the pedigree-based expected relationship under the stated assumptions.

Inbreeding and Kinship Along the Pedigree

The inbreeding coefficient of a nonfounder is the kinship of that individual’s parents. Since \(I\)’s parents are \(F\) and \(H\), the cousin calculation gives \(F_I=\phi_{FH}=1/16\). By contrast, \(F\) and \(H\) themselves have \(F_F=F_H=0\) under the assumed unrelated-parent pairings. Their child \(I\) is inbred because the parents are related; neither parent needs to be inbred for their child to be inbred.

Worked Example: Calculating Inbreeding and Parent–Offspring Kinship

Calculate \(F_I\), then calculate the kinship between \(I\) and \(L\). Individual \(L\) has parents \(I\) and founder \(K\).

The parents of \(I\) are \(F\) and \(H\), whose kinship was calculated as \(1/16\). Therefore,

$$ F_I=\phi_{FH}=\frac{1}{16}=0.0625. $$

The self-kinship of \(I\) is consequently \(\phi_{II}=(1+F_I)/2=(1+1/16)/2=17/32\). Because \(K\) is unrelated to \(I\), \(\phi_{KI}=0\). Apply the recursion to \(L\), whose parents are \(I\) and \(K\), paired with \(I\):

$$ \phi_{LI}=\frac{\phi_{II}+\phi_{KI}}{2} =\frac{17/32+0}{2} =\frac{17}{64} =0.265625. $$

This is the kinship coefficient between parent \(I\) and child \(L\). It is not \(17/32\): that value is the additive relationship \(A_{IL}=2\phi_{IL}\). The distinction matters again for \(M\), whose parents are \(I\) and \(L\). Its inbreeding coefficient is \(F_M=\phi_{IL}=17/64=0.265625\).

The same recursion can fill selected entries of an additive relationship matrix. For example, the relationship between \(I\) and \(F\) can be found from \(I\)’s parents \(F\) and \(H\). Since \(F\) is noninbred, \(\phi_{FF}=1/2\), and the first worked example gives \(\phi_{HF}=1/16\). Therefore,

$$ \phi_{IF}=\frac{\phi_{FF}+\phi_{HF}}{2} =\frac{1/2+1/16}{2} =\frac{9}{32}, \qquad A_{IF}=2\phi_{IF}=\frac{9}{16}. $$

For \(M\), its parents are \(I\) and \(L\). We have \(\phi_{II}=17/32\) and \(\phi_{LI}=17/64\), so \(\phi_{MI}=(17/32+17/64)/2=51/128\), and \(A_{MI}=51/64\). The relationship is larger than the ordinary parent–offspring value for an unrelated mating because \(M\)’s other parent, \(L\), is also related to \(I\).

Worked Example: From Relationships to Genetic Covariances

Worked Example: From Relationships to Genetic Covariances

Suppose the additive genetic variance in a hypothetical trait is \(\sigma_A^2=4\) trait-units squared. Calculate the additive genetic covariance between \(I\) and \(L\), the covariance between \(I\) and \(M\), and the additive genetic variance of \(M\).

Under the additive model, the covariance between two individuals’ additive genetic values is \(\operatorname{Cov}(g_i,g_j)=A_{ij}\sigma_A^2\). For distinct individuals \(I,L\), \(A_{IL}=2\phi_{IL}=2(17/64)=17/32\). Hence,

$$ \operatorname{Cov}(g_I,g_L) =A_{IL}\sigma_A^2 =\frac{17}{32}(4) =\frac{17}{8} =2.125\ \text{trait-units}^2. $$

For \(I\) and \(M\), the relationship calculated above is \(A_{IM}=51/64\). Thus,

$$ \operatorname{Cov}(g_I,g_M) =A_{IM}\sigma_A^2 =\frac{51}{64}(4) =\frac{51}{16} =3.1875\ \text{trait-units}^2. $$

For an individual’s diagonal entry, \(A_{MM}=1+F_M=1+17/64=81/64\). The additive genetic variance of \(M\) is therefore

$$ \operatorname{Var}(g_M) =A_{MM}\sigma_A^2 =\frac{81}{64}(4) =\frac{81}{16} =5.0625\ \text{trait-units}^2. $$

All three results are in squared trait units. They describe additive genetic covariance and variance under this model; they do not include environmental covariance or establish that genetic relatedness causes a difference in the measured trait.

A Compact Matrix Check

A small matrix is a useful way to check that the pedigree calculations use one convention consistently. For individuals \(I,L,M\), the diagonal entries are \(1+F_i\), and the off-diagonal entries are twice the corresponding kinship coefficients. To obtain the remaining entry, use \(M\)’s parents \(I,L\): \(\phi_{ML}=(\phi_{IL}+\phi_{LL})/2=(17/64+1/2)/2=49/128\), since \(L\) is noninbred. Thus \(A_{LM}=49/64\).

$$ \mathbf{A}_{I,L,M}= \begin{pmatrix} 17/16 & 17/32 & 51/64\\ 17/32 & 1 & 49/64\\ 51/64 & 49/64 & 81/64 \end{pmatrix}. $$

The matrix is symmetric, as a covariance or relationship matrix must be. Its diagonal is not all one: the inbred individuals \(I\) and \(M\) have diagonal entries greater than one. Multiplying the whole matrix by \(\sigma_A^2=4\) gives the additive genetic covariance matrix in trait-units squared. This matrix check is not a substitute for tracing parentage carefully, but it can reveal a misplaced entry or a kinship-versus-relationship factor-of-two error.

Common Mistakes and Interpretation

  • Leaving the founder assumption implicit. A founder is not automatically unrelated to every other founder. State the assumption; if founders are related, their kinship entries must reflect that information.
  • Using kinship where a relationship coefficient is required. For distinct individuals, \(A_{ij}=2\phi_{ij}\). Multiplying covariance by \(\phi_{ij}\) alone understates it by a factor of two under the stated additive model.
  • Setting every diagonal relationship to one. The diagonal is \(1+F_i\), so an inbred individual has a larger additive genetic variance than a noninbred individual when the same \(\sigma_A^2\) is used.
  • Confusing an individual’s inbreeding with the parents’ relationship scale. The offspring’s \(F\) equals the parents’ kinship \(\phi\), not twice that value. Twice the kinship is the additive relationship between distinct parents.
  • Reading a genetic covariance as a total trait covariance. The calculation \(\operatorname{Cov}(g_i,g_j)=A_{ij}\sigma_A^2\) is for additive genetic values. A total phenotypic covariance may also involve environmental or other modeled components.
Key takeaway: Specify founder assumptions, work through the parentage in generation order, calculate offspring inbreeding from parental kinship, and convert kinship to additive relationships before scaling by additive genetic variance. Each step has a different quantity and interpretation.

Check Your Understanding

Use the pedigree and assumptions in this tutorial to answer the questions.

  1. Why are \(F\) and \(H\) first cousins, and which individuals are their shared ancestors?
  2. What is \(F_I\), and which coefficient in the pedigree does it equal?
  3. Explain why \(\phi_{IL}=17/64\) is a kinship coefficient, while \(A_{IL}=17/32\) is an additive relationship coefficient.
  4. If \(\sigma_A^2=4\) trait-units squared, what is the additive genetic covariance between \(I\) and \(L\)? State the units.
  5. Why is \(A_{MM}\) greater than one, and what does \(A_{MM}\sigma_A^2\) represent?
Clinical Biostats