Why Ancestry Can Confound a Genetic Association
A genetic association is a relationship between a genetic measurement and a phenotype in observed data. It does not, by itself, show that changing the measured variant would change the phenotype. One reason is that genetic variants are not distributed identically in every population. If ancestry is also related to phenotype—for example, through ancestry-correlated environments or other genetic differences—then ancestry can connect genotype and phenotype even when the variant under study has no causal effect.
This mechanism is called population stratification when differences in population structure create confounding in a genetic association analysis. The central issue is not that ancestry is a biological explanation for every trait difference. Rather, ancestry can index population history and correlate with other influences on the phenotype. The resulting association may reflect those shared differences instead of an effect of the specific variant being tested.
The Causal Diagram
A directed acyclic graph (DAG) represents proposed causal relationships with arrows. In the diagram below, ancestry influences the distribution of the tested genotype and is also related to the phenotype. The arrow from \(G\) to \(Y\) is included to represent a possible causal effect; under the null scenario considered here, that arrow is absent.
The route \(G \leftarrow A \rightarrow Y\) is a backdoor path: it links genotype to phenotype through a common cause rather than through a causal effect of the tested variant. If both arrows from \(A\) are present, that path is open unless the analysis accounts for the relevant ancestry variation. Thus, an observed \(G\)-\(Y\) association can arise even when \(G\) does not cause \(Y\).
A simple structural description makes the point precise. Suppose \(G\) is generated partly by ancestry, and \(Y\) is generated by ancestry and other factors, but not by \(G\):
where \(U_G\) and \(U_Y\) represent other influences. There is no \(G\) argument in the equation for \(Y\), so this model contains no causal effect of the tested genotype on the phenotype. Nevertheless, different values of \(A\) can produce different genotype distributions and different phenotype distributions. Comparing genotype groups that contain different mixtures of ancestry groups can therefore compare people with different background risks.
This explanation depends on the causal structure, not on a claim that ancestry itself must directly cause the phenotype. For example, an ancestry grouping may correlate with environmental exposures, access to care, or other genetic variants that influence \(Y\). A causal diagram is a way to state and examine such assumptions; it is not a discovery that the arrows are true simply because the diagram has been drawn.
Worked Example: Reading the Confounding Path
Suppose a hypothetical study examines a variant and a blood biomarker. The study’s causal assumptions are that ancestry is related to the tested genotype, ancestry-correlated exposures influence the biomarker, and the tested variant has no effect on the biomarker. Identify the source of a possible association.
Step 1: Name the variables. Let \(A\) be the ancestry grouping, \(G\) the tested genotype, and \(Y\) the biomarker level. The assumptions imply \(A \rightarrow G\) and \(A \rightarrow Y\), with no \(G \rightarrow Y\) arrow.
Step 2: Trace the path. The path \(G \leftarrow A \rightarrow Y\) connects genotype and biomarker through ancestry. It is not a causal path from \(G\) to \(Y\), but it can generate an observed association.
Step 3: State the conclusion. If the genotype groups differ in ancestry composition and the biomarker differs across those ancestry groups, an unadjusted genotype–biomarker comparison can reflect the open backdoor path. The association alone does not establish a causal effect of the variant.
How Allele-Frequency Differences Enter
For a diploid autosomal variant, define \(G\) as the number of copies of a specified allele carried by an individual, so \(G\) takes values 0, 1, or 2. Let \(p_a\) be the frequency of that allele in ancestry group \(a\). If the two allele copies are sampled independently with probability \(p_a\), then the expected dosage in that group is \(2p_a\). This simple relationship shows why different allele frequencies create different genotype distributions.
The independence assumption in that calculation corresponds to random pairing of allele copies within the group, as in Hardy–Weinberg proportions. It is used here only to illustrate the expected dosage; real genotype distributions can depart from it because of factors such as inbreeding or nonrandom mating. The core confounding mechanism does not require this assumption: it requires that genotype distributions differ with ancestry and that phenotype distributions also differ with ancestry.
An allele-frequency difference alone is not confounding. If \(G\) varies with ancestry but \(Y\) does not, the ancestry-to-phenotype link needed for this path is absent. Likewise, if phenotype differs with ancestry but the tested genotype distribution does not, this particular path does not connect \(G\) and \(Y\). Confounding is plausible when the relevant links coexist, and the size of the resulting association depends on how strongly they operate and on how ancestry groups are represented in the sample.
Worked Example: From Allele Frequencies to Genotype Dosage
Consider a hypothetical diploid variant in two ancestry groups. The specified allele has frequency \(p_1=0.70\) in group 1 and \(p_2=0.30\) in group 2. Under independent allele copies, calculate the expected dosage in each group and describe what this implies.
For group 1, the expected dosage is
For group 2, it is
The expected dosage difference is \(1.40-0.60=0.80\) copies per person. Therefore, if the sample contains different proportions of the two groups, the genotype distribution can vary across ancestry composition. This calculation does not show that the variant affects any phenotype. To create the confounding path, ancestry must also be related to the phenotype in the study population.
A Hypothetical Pooled Association With No Within-Group Effect
The following example makes both links visible in observed counts. It is a small hypothetical dataset, not a real study. Let \(G\) be the number of copies of an allele, and let \(Y=1\) indicate a binary phenotype. Within each ancestry group, suppose the phenotype rate is the same for all three genotypes. The rates differ between ancestry groups, and genotype counts also differ between groups.
| Ancestry group | Genotype dosage \(G\) | Individuals | Individuals with \(Y=1\) | Phenotype rate |
|---|---|---|---|---|
| Group 1 | 0 | 5 | 1 | 20% |
| Group 1 | 1 | 30 | 6 | 20% |
| Group 1 | 2 | 45 | 9 | 20% |
| Group 2 | 0 | 45 | 27 | 60% |
| Group 2 | 1 | 30 | 18 | 60% |
| Group 2 | 2 | 5 | 3 | 60% |
Worked Example: Combining the Ancestry Groups
Using the hypothetical counts above, calculate the phenotype rate for each genotype after pooling the two ancestry groups. Then compare those rates with the within-group rates.
For \(G=0\): there are \(5+45=50\) individuals, of whom \(1+27=28\) have the phenotype. The pooled rate is \(28/50=0.56\), or 56%.
For \(G=1\): there are \(30+30=60\) individuals, of whom \(6+18=24\) have the phenotype. The pooled rate is \(24/60=0.40\), or 40%.
For \(G=2\): there are \(45+5=50\) individuals, of whom \(9+3=12\) have the phenotype. The pooled rate is \(12/50=0.24\), or 24%.
The pooled phenotype rate decreases from 56% to 24% across genotype dosages, which looks like a strong association. Yet within group 1 the rate is 20% for every dosage, and within group 2 it is 60% for every dosage. In these constructed data, genotype has no within-group association with the phenotype. The pooled pattern arises because the genotype categories have different ancestry compositions, while phenotype rates differ between the ancestry groups. It is not evidence that this allele protects against the phenotype.
What the Example Does—and Does Not—Establish
The example demonstrates a mechanism, not a universal outcome. Population structure can generate a spurious association, but the direction and size of the association are not fixed. They depend on the allele-frequency pattern, phenotype differences across ancestry, and the mixture of groups among the people being compared. Other sources of variation can also be present, including a real causal effect of the tested variant.
An ancestry grouping may also be too coarse to capture the relevant structure. People within one broad group can differ in ancestry, and a sample can include continuous or mixed ancestry. Conversely, group labels can correlate with social or environmental conditions for reasons that vary by setting. A label should not be treated as a complete measurement of ancestry or as a biological cause of a phenotype.
The appropriate conclusion from an observed association is therefore limited: genotype and phenotype differ together in the analyzed sample. To interpret that relationship causally, one must consider whether ancestry or other common causes could produce it, whether the study’s ancestry information captures the relevant structure, and whether the modeling assumptions are credible. This tutorial identifies the source of the concern; choosing and evaluating specific adjustment strategies requires additional tools.
Common Mistakes
- Equating association with causation. “The allele is associated with the phenotype, so the allele causes it” ignores the open path \(G \leftarrow A \rightarrow Y\). The observed association may be due to ancestry-related differences instead.
- Assuming every allele-frequency difference creates confounding. A frequency difference establishes a possible ancestry-to-genotype link, not an ancestry-to-phenotype link. Both are needed for this confounding path.
- Treating ancestry labels as complete or interchangeable with race. Ancestry is multidimensional, and broad labels may not represent the genetic or environmental structure relevant to a particular analysis.
- Assuming a null within-group association guarantees a null pooled association. The worked counts show why that reasoning can fail: genotype groups can contain different proportions of ancestry groups whose phenotype rates differ.
- Interpreting a confounded association as evidence about all members of a population. The pattern depends on the sample’s composition and may not transport to a population with different ancestry proportions or allele frequencies.
Check Your Understanding
Use the causal path and allele-frequency reasoning to answer the following questions.
- In the diagram \(G \leftarrow A \rightarrow Y\), which variable is the common cause, and why can this path produce an association without a causal effect of \(G\) on \(Y\)?
- If the specified allele frequency is 0.40 in an ancestry group, what is the expected dosage under independent allele copies? State the result in copies per person.
- Why does an allele-frequency difference between ancestry groups, by itself, not prove that a genotype–phenotype association is confounded?
- In the hypothetical pooled example, what were the phenotype rates for genotype dosages 0, 1, and 2 within group 1? What does their equality indicate about the within-group association?
- Name one reason that an ancestry grouping might not fully capture the structure relevant to a genetic association analysis.