Tutorials › Statistical Genetics › Transmission Tests for Quantitative Traits

Family-Based Association Tests · Tutorial 266 of 989

Transmission Tests for Quantitative Traits

Build and interpret a quantitative-trait transmission score by comparing each child’s inherited allele count with what the parents’ genotypes predict.

Intermediate 10 min read

What You'll Learn

  • Define the transmitted-allele count and its expectation from the parental genotypes
  • Construct a quantitative-trait score using centered offspring phenotypes
  • Derive the score’s zero expectation under the conditional null
  • Calculate score contributions and conditional variances for informative trios
  • Interpret the score’s direction and recognize when its normal approximation is unreliable

Why Use a Transmission Score for a Quantitative Trait?

A family-based association test can ask whether an allele is transmitted to offspring in a way that tracks a measured trait, such as blood pressure or height. For a quantitative trait, the evidence is not simply whether an offspring inherited a particular allele. The score combines how the child’s inheritance differs from what the parents could have transmitted with how far the child’s trait value lies from a reference level.

The comparison is made within families. Parental genotypes determine the possible transmissions, and Mendelian inheritance supplies their expected frequencies. This makes the score different from a comparison of allele frequencies between unrelated people. It does not, by itself, establish that the tested allele causes the trait: an association can reflect a causal effect or correlation with a nearby causal variant.

Definition: At a biallelic marker, let \(A\) be the allele being tested. For offspring \(i\), let \(G_i\) be the number of \(A\) alleles inherited, so \(G_i\) is 0, 1, or 2. Let \(m_i=E(G_i\mid P_i)\), where \(P_i\) denotes the two parental genotypes. For a fixed trait reference value \(c\), the quantitative transmission score is $$ U=\sum_{i=1}^{n}(G_i-m_i)(Y_i-c), $$ where \(Y_i\) is the offspring’s quantitative phenotype. A trio is informative for this score when its parental genotypes allow \(G_i\) to vary across possible transmissions.

The centering term \(G_i-m_i\) measures the offspring’s excess or deficit of transmitted \(A\) alleles relative to the parental-genotype expectation. The factor \(Y_i-c\) gives that transmission deviation more weight when the child’s phenotype is farther from the reference value. The value \(c\) must be chosen and reported as part of the analysis; it is expressed in the same units as the phenotype. A common choice is a prespecified trait mean or a mean estimated under the null model.

Expected Allele Counts and Informative Parents

At a biallelic marker, write parental genotypes as counts of \(A\): \(AA\) has count 2, \(Aa\) has count 1, and \(aa\) has count 0. Under Mendelian segregation, a heterozygous parent transmits \(A\) with probability \(1/2\); a homozygous parent transmits the same allele in every meiosis. Thus, if the mother and father have \(a_m\) and \(a_f\) copies of \(A\), respectively, the expected transmitted count is

$$ m_i=E(G_i\mid P_i)=\frac{a_m+a_f}{2}. $$

This formula follows because each parent contributes one allele. A homozygous parent’s contribution is fixed, while a heterozygous parent’s expected contribution is \(1/2\). For example, \(AA\times Aa\) gives \(m_i=2/2+1/2=1.5\) only if using parental allele counts as written would be incorrect: each parent’s transmitted-allele expectation is its allele count divided by two, so the correct value is \((a_m+a_f)/2\). For \(AA\times Aa\), this is \((2+1)/2=1.5\). The offspring can have one or two copies of \(A\), with mean 1.5.

If \(h_i\) is the number of heterozygous parents, each contributes a Bernoulli transmission with probability \(1/2\). Assuming the parents’ meioses are independent,

$$ \operatorname{Var}(G_i\mid P_i)=\frac{h_i}{4}. $$

A trio with two homozygous parents has no transmission variation: \(G_i=m_i\) for every possible child, so its score contribution is zero. A trio with one or two heterozygous parents can be informative. Informative does not mean that every such trio contributes a nonzero observed score; the realized transmitted count might equal its expectation.

Worked Example: One Trio With Two Heterozygous Parents

Suppose both parents are \(Aa\), the offspring has genotype \(AA\), and the trait is measured as 72 units. Use \(c=60\) units.

Each parent transmits \(A\) with probability \(1/2\), so the expected offspring count is \(m_i=0.5+0.5=1\). The observed count is \(G_i=2\), giving a transmission deviation of \(G_i-m_i=2-1=1\). The phenotype deviation is \(Y_i-c=72-60=12\) units.

The trio’s score contribution is $$ (G_i-m_i)(Y_i-c)=(2-1)(72-60)=1\times12=12. $$

The contribution is positive: this offspring received one more \(A\) allele than expected from the parents, and the trait value is 12 units above the reference. With two heterozygous parents, the conditional variance of \(G_i\) is \(2/4=0.5\). The contribution’s conditional variance, holding the phenotype fixed, is \(0.5(12)^2=72\) score-units squared.

Why the Conditional Null Expectation Is Zero

The relevant null hypothesis is that, conditional on the parental genotypes, transmission of the tested allele is not associated with the offspring’s phenotype. Under this conditional null, and assuming Mendelian segregation, the observed transmission count has expectation \(m_i\). For a trio whose phenotype is held fixed in the conditioning, \(Y_i-c\) is a constant. Therefore,

$$ E\{(G_i-m_i)(Y_i-c)\mid P_i,Y_i\} =(Y_i-c)E(G_i-m_i\mid P_i,Y_i)=0. $$

The last equality follows because the conditional null gives \(E(G_i\mid P_i,Y_i)=E(G_i\mid P_i)=m_i\). Summing over trios gives \(E(U\mid \text{parental genotypes and phenotypes})=0\). A nonzero observed score is possible under the null because transmissions vary randomly; the zero is the expected score over repetitions of those transmissions, not a claim that every sample score must equal zero.

If transmissions are independent across trios, the conditional variance follows from the same setup. Each phenotype deviation is fixed, so multiplying a transmission deviation by \(Y_i-c\) multiplies its variance by \((Y_i-c)^2\). The variances add across independent trios:

$$ \operatorname{Var}(U\mid P,Y) =\sum_{i=1}^{n}(Y_i-c)^2\operatorname{Var}(G_i\mid P_i) =\sum_{i=1}^{n}(Y_i-c)^2\frac{h_i}{4}. $$

This variance is conditional on the parental genotypes and the observed phenotypes. It is not the variance of the trait itself. A standardized score can be formed as \(Z=U/\sqrt{\operatorname{Var}(U\mid P,Y)}\), provided the denominator is positive. A standard normal approximation for \(Z\) is most credible when information is spread across many independent transmissions; a few trios do not generally justify relying on that approximation for a p-value.

Worked Example: Summing Scores Across Four Hypothetical Trios

Consider four hypothetical trios, with trait values in units and \(c=50\) units. The table gives the offspring genotype, so \(G_i\) is the observed number of \(A\) alleles.

ParentsOffspring\(Y_i\)\(m_i\)\((G_i-m_i)(Y_i-c)\)Variance contribution
\(AA\times Aa\)\(Aa\)621.5\((1-1.5)(62-50)=-6\)\((12)^2(1/4)=36\)
\(Aa\times aa\)\(aa\)380.5\((0-0.5)(38-50)=6\)\((-12)^2(1/4)=36\)
\(Aa\times Aa\)\(Aa\)701\((1-1)(70-50)=0\)\((20)^2(2/4)=200\)
\(AA\times AA\)\(AA\)552\((2-2)(55-50)=0\)\((5)^2(0)=0\)

For the first trio, the \(AA\) parent always transmits \(A\), while the \(Aa\) parent transmits it with probability \(1/2\). Thus \(m_1=1+0.5=1.5\), and the observed \(Aa\) child has \(G_1=1\). Its contribution is \((-0.5)(12)=-6\). For the second trio, \(m_2=0.5\), while the \(aa\) child has \(G_2=0\); its contribution is \((-0.5)(-12)=6\).

The third trio has expected and observed count both equal to 1, so its score contribution is zero, although its transmissions are variable and its variance contribution is positive. The fourth trio has no heterozygous parent, so its transmission is fixed and it contributes neither score nor variance.

Adding the contributions gives \(U=-6+6+0+0=0\). Under independent meioses, the conditional variance is \(36+36+200+0=272\) score-units squared, and the standardized score is \(0/\sqrt{272}=0\). This small hypothetical set is used to show the arithmetic, not to claim that a normal approximation or an inferential conclusion is warranted.

Direction, Scale, and Interpretation

The sign of \(U\) summarizes how transmission deviations align with trait deviations from \(c\). A positive contribution results when an offspring above \(c\) inherits more \(A\) alleles than expected, or when an offspring below \(c\) inherits fewer than expected. A negative contribution results from the opposite alignments. The total score combines these contributions across families; it is not an estimated change in the trait per allele and should not be reported in trait units alone. Its units are allele-count deviations multiplied by trait units.

Worked Example: Contributions That Reinforce Each Other

Take \(c=50\) units and two hypothetical \(Aa\times aa\) trios. In the first, the child is \(Aa\), so \(G_1=1\), and has \(Y_1=65\). In the second, the child is \(aa\), so \(G_2=0\), and has \(Y_2=35\).

For each trio, one parent is heterozygous and the other is \(aa\), so \(m_1=m_2=0.5\). The first contribution is $$ (1-0.5)(65-50)=0.5\times15=7.5. $$

The second contribution is $$ (0-0.5)(35-50)=(-0.5)(-15)=7.5. $$

The combined score is \(U=7.5+7.5=15\) allele-count-deviation times trait units. Both trios contribute positively: the higher-than-reference value coincides with excess \(A\) transmission in the first, and the lower-than-reference value coincides with deficient \(A\) transmission in the second. This describes alignment in these hypothetical families; it does not establish a causal effect.

Assumptions and a Common Pitfall

Interpreting the score requires that the parental genotypes and offspring phenotype be measured correctly, that the relevant transmissions follow Mendelian probabilities under the null, and that the conditional null of no transmission–phenotype association is appropriate. The variance formula above additionally assumes independent meioses across the trios being summed. Related offspring or other sources of dependence require a variance calculation that accounts for that dependence; simply adding individual variances would then be unjustified.

A positive score is evidence in a particular direction only relative to the chosen allele \(A\), the trait reference \(c\), and the stated score definition. Changing which allele is tested reverses the sign of each transmission deviation. Changing \(c\) can also change the score, so the reference must be chosen transparently rather than selected after inspecting which value gives the desired result.

  • Mistake: “The parents are heterozygous, so the child’s observed allele count is the expected count.” Correction: Heterozygous parents define a probability distribution over transmissions; the child’s realized count can differ from its expectation.
  • Mistake: “A zero contribution means the trio provided no transmission information.” Correction: The observed deviation \(G_i-m_i\) may be zero even when the parental genotypes allow more than one transmission outcome. In that case the trio can still contribute to the score variance.
  • Mistake: “A positive score proves that allele \(A\) raises the trait.” Correction: The score summarizes an association pattern under a specified coding and reference. It is not a causal effect estimate, and linkage to another variant can also produce a pattern.
  • Mistake: “The standardized score is automatically standard normal.” Correction: Standardization alone does not ensure a normal null distribution. The approximation needs enough independent transmission information and should not be trusted uncritically in a very small sample.
Key takeaway: For each trio, subtract the parental-genotype expectation from the child’s transmitted \(A\)-allele count, multiply by the child’s trait deviation from a specified reference, and sum. Under Mendelian transmission and the conditional null, the score has expectation zero; informative transmissions determine its conditional variance.

Check Your Understanding

Use the definitions and conditional reasoning above to answer the following questions.

  1. For parents \(Aa\times aa\), what are the expected transmitted \(A\)-allele count and its conditional variance?
  2. If a child’s observed count equals \(m_i\), must the trio be uninformative? Explain.
  3. For a trio with one heterozygous parent, \(Y_i-c=-8\) units, and \(G_i-m_i=0.5\), calculate the score contribution and state its units.
  4. Why is the conditional expectation of a trio’s score contribution zero under the null?
  5. What additional assumption allows the conditional variances of separate trios to be added using the stated formula?
Clinical Biostats