Tutorials › Statistical Genetics › Choosing a Coding From Genotype Patterns

Genetic Effects and Trait Models · Tutorial 29 of 989

Choosing a Coding From Genotype Patterns

Use the equality patterns among genotype means to choose a coding that represents the observed trait pattern without imposing unsupported constraints.

Beginner 10 min read

What You'll Learn

  • Express additive, dominant, recessive, and unrestricted models as constraints on three genotype means
  • Compare adjacent mean differences to check for an additive pattern
  • Identify the mean equalities implied by dominant and recessive coding
  • Choose a two-degree-of-freedom genotypic coding when simpler patterns do not fit
  • Distinguish a descriptive pattern in estimated means from evidence that a population model is true

Reading a Genotype-Mean Table

A genotype-mean table puts three trait means side by side: one for each genotype at a biallelic locus. Its pattern can suggest a useful coding for the genetic effect. The key is to compare the means, not simply to choose a familiar coding or to focus on which genotype has the largest mean.

Use \(aa\), \(Aa\), and \(AA\) for the three genotypes, and let \(m_0\), \(m_1\), and \(m_2\) denote their respective mean trait values. The dosage of allele \(A\) is 0 for \(aa\), 1 for \(Aa\), and 2 for \(AA\). Throughout, the trait is quantitative and the genotype-specific means are measured on the same scale and in the same units.

Definition: A genotype-mean pattern is additive when the heterozygote mean lies exactly halfway between the two homozygote means; dominant for allele \(A\) when the \(Aa\) and \(AA\) means are equal; and recessive for allele \(A\) when the \(aa\) and \(Aa\) means are equal. An unrestricted genotypic model allows all three means to differ without imposing these equalities.

These labels describe restrictions on genotype means. They do not, by themselves, establish a biological mechanism or show that a genotype causes a trait difference. In real data, the means are usually estimated from individuals; a close visual pattern is not proof that the corresponding population means are exactly equal.

The Mean Constraints Behind Each Coding

The usual coding schemes make their assumptions visible as equations. Under additive coding, the mean changes by the same amount for each additional copy of \(A\). Under dominant or recessive coding, two genotypes are assigned the same coded value, so the model requires their means to be equal. Unrestricted coding makes no equality constraint.

$$ \begin{aligned} \text{Additive:}\quad &m_0=\mu,\quad m_1=\mu+\beta,\quad m_2=\mu+2\beta,\\ \text{Dominant for }A:\quad &m_0=\mu,\quad m_1=m_2=\mu+\delta,\\ \text{Recessive for }A:\quad &m_0=m_1=\mu,\quad m_2=\mu+\delta,\\ \text{Unrestricted:}\quad &m_0=\mu,\quad m_1=\mu+\delta_1,\quad m_2=\mu+\delta_2. \end{aligned} $$

Here \(\mu\) is the mean for the reference genotype in the stated coding. In the additive model, \(\beta\) is the change in mean per additional \(A\) allele. In the dominant and recessive models, \(\delta\) is the difference between the two coded groups. In the unrestricted model, \(\delta_1\) and \(\delta_2\) are separate differences from the \(aa\) mean. Their units are the same as the trait.

The patterns can also be checked directly using differences. Define the first step as \(m_1-m_0\) and the second as \(m_2-m_1\). Additivity requires the steps to be equal. Dominance for \(A\) requires the second step to be zero, while recessiveness for \(A\) requires the first step to be zero. These are exact statements about means; with estimated means, the observed differences may only be approximately equal to zero or to each other.

Worked Example: An Additive Pattern

Suppose a hypothetical study reports mean systolic blood pressure values of 110, 114, and 118 mmHg for \(aa\), \(Aa\), and \(AA\), respectively. Check whether additive coding matches these means.

The two adjacent steps are \(114-110=4\) mmHg and \(118-114=4\) mmHg. They are equal, so the heterozygote mean is halfway between the two homozygote means:

$$ m_1-m_0=114-110=4,\qquad m_2-m_1=118-114=4. $$

The additive model fits this table with \(\mu=110\) mmHg and \(\beta=4\) mmHg per \(A\) allele. Substitution gives \(m_0=110\), \(m_1=110+4=114\), and \(m_2=110+2(4)=118\) mmHg. This pattern is consistent with additive coding; the table alone does not establish that the population means follow an exact additive model.

Worked Example: A Dominant Pattern

Consider hypothetical mean enzyme activities of 5, 12, and 12 units per liter for \(aa\), \(Aa\), and \(AA\). Determine whether a dominant coding for \(A\) represents these means.

The heterozygote and \(AA\) means are equal: \(m_1=m_2=12\) units per liter. The dominant model groups both genotypes carrying \(A\) together. With \(aa\) as the reference, \(\mu=5\) units per liter and \(\delta=12-5=7\) units per liter. The predicted means are \(5\), \(5+7=12\), and \(5+7=12\) units per liter. They match the table exactly.

The first step is \(12-5=7\) units per liter, while the second is \(12-12=0\). That zero second step is the defining equality for dominance of \(A\) in this mean pattern.

Worked Example: A Recessive Pattern

Suppose a hypothetical trait has means of 30, 30, and 36 millimeters for \(aa\), \(Aa\), and \(AA\). Check the recessive coding for \(A\).

Here \(m_0=m_1=30\) mm, so the two genotypes with zero or one copy of \(A\) share a mean. The recessive coding uses an indicator that is 1 only for \(AA\). Its reference mean is \(\mu=30\) mm, and the difference for \(AA\) is \(\delta=36-30=6\) mm. It predicts \(30\), \(30\), and \(30+6=36\) mm, matching all three entries.

The adjacent steps confirm the pattern: \(30-30=0\) mm and \(36-30=6\) mm. The zero first step, rather than the size or direction of the final mean alone, identifies the recessive mean pattern for \(A\).

When the Three Means Need Separate Effects

Some tables do not satisfy any of the three simpler patterns. For instance, if the heterozygote mean is neither halfway between the homozygote means nor equal to either one, an additive, dominant, or recessive coding would impose a constraint not present in the table. An unrestricted genotypic model then represents the three means directly, using \(aa\) as the reference and separate contrasts for \(Aa\) and \(AA\).

Worked Example: Choosing Unrestricted Genotypic Coding

Suppose a hypothetical trait has means of 20, 21, and 29 centimeters for \(aa\), \(Aa\), and \(AA\). Compare the patterns and select a coding that represents the table without a false equality.

The adjacent steps are \(21-20=1\) cm and \(29-21=8\) cm, so they are not equal; the additive pattern does not fit. Also, \(m_1=21\) cm is not \(m_2=29\) cm, ruling out the exact dominant pattern, and \(m_0=20\) cm is not \(m_1=21\) cm, ruling out the exact recessive pattern.

The unrestricted model represents the means with \(\mu=20\) cm, \(\delta_1=21-20=1\) cm, and \(\delta_2=29-20=9\) cm. Substitution gives \(m_0=20\), \(m_1=20+1=21\), and \(m_2=20+9=29\) cm. The two separate contrasts reproduce all three means without requiring either adjacent difference to be equal or either pair of means to match.

A Short Justification: Why These Checks Work

The patterns follow directly from the coded predictor values. In an additive model, genotype dosage \(x\) takes values \(0,1,2\), and the mean is \(\mu+\beta x\). Thus the three predicted means are \(\mu\), \(\mu+\beta\), and \(\mu+2\beta\). Subtracting consecutive means gives \(\beta\) both times, so equal steps are necessary. Conversely, if the two steps are equal to some value \(\beta\), setting \(\mu=m_0\) gives exactly those three means under the additive formula.

For dominant coding, the predictor is 0 for \(aa\) and 1 for both \(Aa\) and \(AA\). A linear mean model in this predictor must give the same predicted mean to those latter genotypes. For recessive coding, the predictor is 0 for \(aa\) and \(Aa\), and 1 for \(AA\), requiring the first two means to match. In unrestricted coding, use separate indicators for \(Aa\) and \(AA\), with \(aa\) as the reference. The model has two contrasts and can therefore reproduce any three specified means: choose \(\mu=m_0\), \(\delta_1=m_1-m_0\), and \(\delta_2=m_2-m_0\).

$$ \begin{aligned} \text{Additive check:}\quad &m_1-m_0=m_2-m_1,\\ \text{Dominant-}A\text{ check:}\quad &m_1=m_2,\\ \text{Recessive-}A\text{ check:}\quad &m_0=m_1. \end{aligned} $$

Choosing a Coding Carefully

A practical decision begins with the scientific question and the pattern of the means. If the aim is to describe a pattern that follows allele dosage, check whether the adjacent steps are similar. If a carrier effect is the intended representation, check whether the two \(A\)-carrying genotypes have similar means. If an effect is expected only for \(AA\), check whether \(aa\) and \(Aa\) have similar means. Use unrestricted genotypic coding when you need to allow the three genotypes to differ without those constraints.

For estimated means, “similar” is not the same as “exactly equal.” Sampling fluctuation can make estimates differ even when population means are equal, or make them appear equal when they are not. A table is useful for proposing or describing a coding, but deciding whether a restriction is statistically compatible with the data requires accounting for uncertainty. Also check that all three genotype groups are represented and that their means are defined; an absent genotype group cannot supply an observed mean for this comparison.

The allele named in the pattern matters. “Dominant for \(A\)” means that \(Aa\) and \(AA\) share a mean; “recessive for \(A\)” means that \(AA\) differs while \(aa\) and \(Aa\) share a mean. Reversing which allele is named changes which two genotypes are grouped. State the allele and the required equality rather than relying on the label alone.

Common Mistake

A common mistake is to call a pattern additive because the \(AA\) mean is larger than the \(Aa\) mean, or dominant because the two means look close after rounding. Neither conclusion follows from that comparison alone. Additivity requires equal adjacent steps, and dominance for \(A\) requires equality of \(m_1\) and \(m_2\). For example, means of 20, 21, and 29 cm have an increasing trend, but their steps are 1 and 8 cm, and neither dominant nor recessive equality holds. The correction is to write down the relevant differences or equality before choosing the coding.

Key takeaway: Match the coding to the mean constraints it imposes: equal adjacent steps for additive coding, equal \(Aa\) and \(AA\) means for dominance of \(A\), equal \(aa\) and \(Aa\) means for recessiveness of \(A\), and separate contrasts for unrestricted genotypic coding. Treat observed matches as descriptive until uncertainty has been considered.

Check Your Understanding

For each question, focus on the genotype-mean constraints rather than on the labels alone.

  1. For means \(8\), \(11\), and \(14\) units for \(aa\), \(Aa\), and \(AA\), calculate both adjacent steps and identify the matching pattern.
  2. What equality between genotype means is required by a dominant coding for \(A\)?
  3. For means \(40\), \(40\), and \(35\) points, which simpler coding pattern matches, and what is the difference from the reference mean for \(AA\)?
  4. Why can an unrestricted genotypic model reproduce any three genotype means when \(aa\) is the reference?
  5. Why should approximate equality among sample means not automatically be treated as proof of a population-level dominant or recessive effect?
Clinical Biostats