Tutorials › Statistical Genetics › Affected-Pedigree-Member Sharing Statistics

Nonparametric Linkage Methods · Tutorial 251 of 989

Affected-Pedigree-Member Sharing Statistics

Calculate and interpret weighted allele-sharing statistics while keeping marker frequencies and pedigree structure in the null model.

Intermediate 10 min read

What You'll Learn

  • Define a weighted allele-sharing score for a pair of affected relatives
  • Explain why rare alleles can contribute more to an affected-pedigree-member statistic
  • Calculate expected sharing for independent marker-allele copies
  • Explain why related relatives require a pedigree-specific null distribution
  • Standardize an observed sharing score and interpret it without overclaiming
  • Identify why allele sharing by state is not the same as sharing identical by descent

Why Compare Sharing With a Pedigree-Specific Expectation?

If affected relatives share marker alleles more often than expected, that can be evidence that a marker is linked to a disease locus: relatives who inherit the same chromosomal segment may also inherit the same marker alleles. An affected-pedigree-member (APM) statistic summarizes this sharing across affected relatives. Its central comparison is not simply “how much sharing was observed?” but “how much sharing was observed relative to what this marker and this pedigree would produce under the null model?”

Two features matter from the outset. First, two allele copies can match by state (IBS) without having been inherited from the same ancestral copy; marker genotypes alone do not make those concepts interchangeable. Second, relatives are not independent observations. Their inheritance patterns are connected by the pedigree, so an appropriate null expectation and variance must account for that structure.

Definition: An APM sharing statistic is a summary of marker-allele sharing among affected members of a pedigree, often with each shared allele weighted according to its population frequency. Evidence is assessed by comparing the observed statistic with its expectation under a specified no-linkage model that retains the marker allele frequencies and pedigree inheritance structure.

A Weighted Pairwise Sharing Score

One common way to construct an APM statistic begins with a score for each pair of affected relatives. Let \(p_a\) denote the population frequency of marker allele \(a\). For illustration, use the inverse-square-root weight \(w(p_a)=1/\sqrt{p_a}\): a match involving a less frequent allele receives a larger contribution. Other APM implementations can use different weighting functions, so the weight must be stated when a statistic is reported.

Suppose relatives \(i\) and \(j\) each have two observed marker-allele copies. Define their IBS score as the sum over the four cross-person copy comparisons, adding \(w(p_a)\) for each pair of copies that both carry allele \(a\). If an allele occurs twice in one person and once in the other, it contributes twice. This copy-pair convention makes the score explicit and avoids ambiguity about whether “sharing” means one match or two.

$$ S_{ij}=\sum_{r=1}^{2}\sum_{s=1}^{2} \mathbf{1}(G_{ir}=G_{js})\,w\!\left(p_{G_{ir}}\right), \qquad T=\sum_{i<j;\ i,j\ \mathrm{affected}} S_{ij}. $$

Here \(G_{ir}\) is copy \(r\) in person \(i\), and \(\mathbf{1}(\cdot)\) is 1 when its condition is true and 0 otherwise. \(T\) is the total pairwise sharing score for the affected members being analyzed. The same pair can contribute several matching copy pairs, and a person in a larger affected pedigree can participate in many pairs.

Worked Example: Scoring One Pair of Genotypes

Consider a hypothetical marker with alleles A and B, with frequencies \(p_A=0.75\) and \(p_B=0.25\). Two affected relatives have genotypes AA and AB. Use \(w(p)=1/\sqrt{p}\) and the four-copy-comparison score above.

The two copies in the first relative are A, and the second relative has one A and one B. There are two matching A copy pairs and no matching B copy pair. Thus

$$ S_{ij}=2\left(\frac{1}{\sqrt{0.75}}\right)+0\left(\frac{1}{\sqrt{0.25}}\right) =2(1.1547)=2.309. $$

The result is a weighted IBS score, not a probability and not a count of copies known to be identical by descent (IBD). The weighting makes the A matches contribute their specified frequency-based amount.

Expected Sharing Depends on Frequencies

A match is more likely for a common allele than for a rare one. The null expectation therefore has to account for allele frequencies; otherwise, a marker with many common alleles could appear to show excess sharing simply because matches are frequent at that marker.

For a simple baseline, suppose the two copies being compared are independent draws from the marker population. The probability that both copies are allele \(a\) is \(p_a^2\). Since a match on \(a\) contributes \(w(p_a)\), the expected weighted contribution from that comparison is \(p_a^2w(p_a)\). There are four cross-person copy comparisons in \(S_{ij}\), so the expected pair score for two independent genotypes is \(4\sum_a p_a^2w(p_a)\). This baseline is useful for understanding the frequency effect, but it is not the correct null expectation for related pedigree members.

Worked Example: A Frequency-Based Baseline

For a hypothetical two-allele marker, let \(p_A=0.75\), \(p_B=0.25\), and \(w(p)=1/\sqrt p\). Suppose, just for this baseline calculation, that the four cross-person copy comparisons each involve independent population draws. The expected pair score is

$$ E(S_{ij})=4\left[p_A^2\frac{1}{\sqrt{p_A}}+ p_B^2\frac{1}{\sqrt{p_B}}\right] =4\left[0.75^{3/2}+0.25^{3/2}\right]. $$

Since \(0.75^{3/2}\approx0.6495\) and \(0.25^{3/2}=0.1250\),

$$ E(S_{ij})=4(0.7745)=3.098. $$

This is the independent-copy baseline under the stated frequencies and weight. It should not be substituted for the null mean of a pair of siblings, whose marker alleles are correlated through their parents and transmissions.

Building the Null Model Around the Pedigree

For a pedigree, the null model must preserve the relationships among its members. A standard approach specifies the marker allele frequencies in the founders and the Mendelian transmission process through the pedigree. Under the no-linkage null, disease status does not select marker inheritance patterns through linkage to the disease locus; nevertheless, the pedigree still induces dependence among marker genotypes. The expected score and its variance must be computed under those pedigree inheritance rules, not by treating every relative pair as an independent pair of population draws.

In practice, the null distribution can be obtained by enumerating possible founder genotypes and transmissions when feasible, or by simulation from the specified pedigree model. For every null replicate, calculate the same \(T\) from the resulting genotypes. The null mean is \(E_0(T)\), and the null variance is \(\operatorname{Var}_0(T)=E_0[(T-E_0(T))^2]\). The subscript 0 indicates that these quantities are calculated under the no-linkage null.

Key result: A standardized APM score is \(Z=[T_{\mathrm{obs}}-E_0(T)]/\sqrt{\operatorname{Var}_0(T)}\), provided the null variance is positive. This reports the observed excess or deficit in null standard-deviation units. A normal approximation to \(Z\) is not automatically valid; the null distribution may be discrete or skewed, especially in a small pedigree.

Worked Example: Standardizing a Three-Sibling Score

Consider a hypothetical pedigree of three affected siblings typed at a two-allele marker with frequencies 0.80 and 0.20. Under a specified Mendelian inheritance model, suppose exact enumeration of the null inheritance patterns for the chosen weighted score gives \(E_0(T)=15\) and \(E_0(T^2)=262.5\). These are pedigree- and marker-model-specific moments, not universal constants.

First calculate the null variance:

$$ \operatorname{Var}_0(T)=E_0(T^2)-[E_0(T)]^2 =262.5-15^2=262.5-225=37.5. $$

Suppose the hypothetical observed score is \(T_{\mathrm{obs}}=18.75\). Its standardized excess is

$$ Z=\frac{18.75-15}{\sqrt{37.5}} =\frac{3.75}{6.1237} \approx0.6124. $$

Thus the observed score is about 0.6124 null standard deviations above its mean. This calculation alone does not establish a p-value: the exact null distribution or a justified approximation is needed for that inference. In particular, using a variance that ignores the shared parental transmissions would change the standardization and could misstate the evidence.

Why Pedigree Structure Changes the Answer

If two affected siblings share a marker allele, their match is not surprising in the same way as a match between unrelated people. Siblings have a chance of inheriting the same parental allele copy, and multiple sibling pairs in a larger family are dependent because they reuse the same individuals and parental transmissions. Consequently, the number of affected pairs is not the number of independent observations. Increasing the number of affected relatives can increase \(T\), but it can also change the null mean and variance.

The marker itself matters as well. With inverse-square-root weighting, a match on an allele with frequency 0.10 contributes \(1/\sqrt{0.10}\approx3.162\), while a match on an allele with frequency 0.90 contributes \(1/\sqrt{0.90}\approx1.054\). Rare-allele matches therefore count more under this weighting rule. But a rare allele also has a different probability of appearing in the pedigree, so its contribution must be evaluated against a null distribution using the same allele frequencies. A large weighted score is not, by itself, evidence of linkage.

The no-linkage model also determines what is conditioned on. The pedigree, observed affected-member set, marker frequencies, missing genotypes, and any ascertainment assumptions should be represented appropriately in the null calculation. Changing those choices can change the distribution of \(T\). Thus a reported APM result is interpretable only alongside the scoring rule and the null model used to calibrate it.

A Short Justification for the Null Mean

The frequency contribution to the independent-copy baseline follows directly by summing over mutually exclusive matching alleles. For one comparison, the event “both copies are \(a\)” has probability \(p_a^2\), and its score contribution is \(w(p_a)\). The expectation of that comparison is therefore the sum of the score times the probability over all alleles:

$$ E\!\left[\mathbf{1}(G_1=G_2)w(p_{G_1})\right] =\sum_a P(G_1=a,G_2=a)w(p_a) =\sum_a p_a^2w(p_a). $$

Adding the expectations for four comparisons gives \(4\sum_a p_a^2w(p_a)\) by linearity of expectation; independence among the four comparisons is not required for that mean calculation. For a pedigree, the same logic applies once the relevant joint probabilities are calculated under its inheritance model. Variance is more demanding: it includes covariances among score components, so it cannot generally be found by adding pairwise variances as if relatives were independent.

Common Mistakes and Practical Interpretation

  • Using the unrelated-person baseline for relatives. The formula \(4\sum_a p_a^2w(p_a)\) assumes independent allele-copy draws. Use a pedigree-specific null for siblings and other relatives.
  • Treating every pair as independent. Pairs in a family overlap in people and inherited alleles. Their dependence affects the variance of the total score.
  • Confusing IBS with IBD. Equal observed marker alleles are IBS. They can be shared by descent, but genotype matching alone does not prove that they are.
  • Ignoring the weighting rule. A score using \(1/\sqrt p\) is not interchangeable with an unweighted count or a statistic using a different weight.
  • Reading a standardized score as a p-value. \(Z\) gives a distance from the null mean in null standard-deviation units. A p-value requires a calibrated null distribution and a stated tail direction.
  • Interpreting excess sharing as causation. Evidence of linkage concerns co-inheritance in the studied pedigrees; it does not show that the marker allele causes disease.
Key takeaway: An APM statistic summarizes weighted marker-allele sharing among affected relatives. Its interpretation depends jointly on the allele-frequency weighting, the observed sharing, and a null distribution that respects the pedigree’s inheritance and dependence structure.

Check Your Understanding

Use the stated weighting rule and distinguish the scoring step from the null-model step.

  1. For \(w(p)=1/\sqrt p\), how much does one match contribute when its allele frequency is 0.16?
  2. Why is \(4\sum_a p_a^2w(p_a)\) not generally the correct expected score for a sibling pair?
  3. What does \(Z=1.2\) say about an observed score, and what additional information is needed to obtain a p-value?
  4. Explain how two marker alleles can be IBS without being IBD.
  5. Why can adding affected relatives increase both the observed total score and its null variance?
Clinical Biostats