Tutorials › Statistical Genetics › Population Structure As Instrument Confounding

Mendelian Randomization Foundations · Tutorial 845 of 989

Population Structure As Instrument Confounding

See how population structure can confound a genetic instrument, and how ancestry-aware analysis and study design can address the problem.

Advanced 10 min read

What You'll Learn

  • Explain how ancestry can associate an instrument with potential outcomes through differences in allele frequencies and outcome risks.
  • Distinguish population-structure confounding from horizontal and vertical pleiotropy.
  • Calculate how ancestry-specific risks combine into a confounded marginal instrument–outcome association.
  • Interpret ancestry adjustment and standardization as ways to compare like ancestry distributions.
  • Identify design safeguards and their limitations, including ancestry restriction and within-family analysis.

Why Population Structure Matters for Mendelian Randomization

A genetic variant is not automatically independent of the factors that influence an outcome. In a population with ancestry-related structure, allele frequencies can differ across ancestry groups, while outcome risks can also differ because of environmental, social, geographic, or other background factors. The variant can then be associated with outcome risk through ancestry, even when the variant itself has no causal effect on the outcome.

This is a violation of the instrument independence assumption, not necessarily evidence that the variant acts through an unintended biological pathway. The distinction matters: population structure creates a noncausal association through differences between groups, whereas horizontal pleiotropy concerns a variant’s effect on the outcome through a pathway other than the exposure of interest.

Definition: Population structure is systematic variation in ancestry or related population background among individuals in a study. It is a source of instrument confounding when the instrument \(Z\) is associated with ancestry \(A\), and ancestry is associated with the outcome or its potential outcomes. In that case, the instrument may not be independent of the outcome determinants, as required for a valid instrumental-variable analysis.

For clarity, consider a binary instrument \(Z\), ancestry group \(A\), exposure \(X\), and outcome \(Y\). The same issue applies to a genetic dosage, such as the number of effect alleles, rather than a binary instrument. In a simplified causal framework, the instrument should be relevant to \(X\), independent of outcome-related potential outcomes, and have no effect on \(Y\) except through \(X\). Population structure threatens independence if \(Z\) and outcome determinants vary together across ancestry groups. It can also make the instrument–exposure association reflect ancestry differences, so relevance and interpretation should be assessed in the population and analysis being used.

How Ancestry Creates an Instrument–Outcome Association

Suppose allele frequencies differ across ancestry groups. Then \(P(A=a\mid Z=1)\) may differ from \(P(A=a\mid Z=0)\). If outcome risk also varies by ancestry, the two instrument groups have different mixtures of background outcome risk. Consequently, their marginal outcome risks can differ even if the instrument has no causal effect on the outcome. If the instrument does affect the exposure, that marginal association may combine the intended causal pathway with this ancestry-related imbalance.

A useful decomposition makes the source of the difference explicit. Let \(m_{za}=E[Y\mid Z=z,A=a]\), using a binary outcome so that this conditional mean is a risk. The marginal difference in outcome risks can be written by adding and subtracting the risk under \(Z=0\) averaged over the ancestry distribution among those with \(Z=1\):

$$ \begin{aligned} E[Y\mid Z=1]-E[Y\mid Z=0] &=\sum_a m_{1a}P(A=a\mid Z=1)-\sum_a m_{0a}P(A=a\mid Z=0)\\ &=\sum_a (m_{1a}-m_{0a})P(A=a\mid Z=1)\\ &\quad+\sum_a m_{0a}\{P(A=a\mid Z=1)-P(A=a\mid Z=0)\}. \end{aligned} $$

The first term compares instrument groups within ancestry strata, with a common set of weights. The second is a composition term: it arises because ancestry distributions differ between the instrument groups, and it is nonzero when ancestry-specific baseline risks differ. This is an associational decomposition; a causal interpretation of the first term requires the relevant assumptions to hold within strata.

Worked Example: A Marginal Association With Within-Ancestry Effects

Consider hypothetical data with two ancestry groups, \(A_1\) and \(A_2\). The joint proportions of ancestry and instrument status are \(P(A_1,Z=1)=0.30\), \(P(A_1,Z=0)=0.20\), \(P(A_2,Z=1)=0.05\), and \(P(A_2,Z=0)=0.45\). Suppose the outcome risks are 0.30 and 0.20 for \(Z=1\) and \(Z=0\), respectively, in \(A_1\); in \(A_2\), they are 0.15 and 0.05. These are hypothetical risks, not results from a real study.

The instrument prevalence is \(P(Z=1)=0.30+0.05=0.35\), and \(P(Z=0)=0.20+0.45=0.65\). Thus, among \(Z=1\), the ancestry proportions are \(0.30/0.35\) for \(A_1\) and \(0.05/0.35\) for \(A_2\). Among \(Z=0\), they are \(0.20/0.65\) and \(0.45/0.65\), respectively.

The marginal risk among those with \(Z=1\) is:

$$ \left(\frac{0.30}{0.35}\right)(0.30)+ \left(\frac{0.05}{0.35}\right)(0.15) =\frac{0.0975}{0.35} \approx 0.2786. $$

The marginal risk among those with \(Z=0\) is:

$$ \left(\frac{0.20}{0.65}\right)(0.20)+ \left(\frac{0.45}{0.65}\right)(0.05) =\frac{0.0625}{0.65} \approx 0.09615. $$

The marginal risk difference is \(0.2785714-0.0961538\approx 0.1824\), or about 18.24 percentage points. Yet the within-ancestry risk differences are \(0.30-0.20=0.10\) in \(A_1\) and \(0.15-0.05=0.10\) in \(A_2\). The marginal contrast is larger because the \(Z=1\) group contains more people from the higher-risk ancestry group. This example shows how a marginal instrument–outcome association can be distorted by population structure even when the within-ancestry contrasts agree.

A Pure Confounding Example

Population structure can produce an instrument–outcome association even when there is no within-ancestry association at all. The next example isolates that possibility: outcome risk depends on ancestry but not on instrument status within either ancestry group.

Worked Example: An Association With No Within-Ancestry Difference

Use the same hypothetical joint proportions: 0.30 in \(A_1\) with \(Z=1\), 0.20 in \(A_1\) with \(Z=0\), 0.05 in \(A_2\) with \(Z=1\), and 0.45 in \(A_2\) with \(Z=0\). Suppose the outcome risk is 0.20 in \(A_1\) for both instrument groups, and 0.05 in \(A_2\) for both instrument groups.

Among \(Z=1\), the marginal risk is:

$$ \left(\frac{0.30}{0.35}\right)(0.20)+ \left(\frac{0.05}{0.35}\right)(0.05) =\frac{0.0625}{0.35} \approx 0.1786. $$

Among \(Z=0\), it is:

$$ \left(\frac{0.20}{0.65}\right)(0.20)+ \left(\frac{0.45}{0.65}\right)(0.05) =\frac{0.0625}{0.65} \approx 0.09615. $$

The marginal risk difference is \(0.1785714-0.0961538\approx 0.08242\), or 8.242 percentage points. Within \(A_1\), the risk difference is \(0.20-0.20=0\); within \(A_2\), it is \(0.05-0.05=0\). The marginal association arises entirely because \(A_1\), the higher-risk group in this example, is more common among \(Z=1\) than among \(Z=0\). It is not evidence of an effect of \(Z\) on \(Y\).

What Ancestry Adjustment Does

Adjustment for ancestry aims to compare instrument groups at the same measured ancestry distribution, rather than comparing groups whose ancestry composition differs. In practice, ancestry may be represented using measured covariates such as self-reported background, genetic principal components, geographic variables, or combinations of these, depending on the study. These measures are not interchangeable, and no single adjustment automatically captures every relevant source of structure.

One way to express the idea is direct standardization. Choose common weights \(q(a)\) that sum to one across ancestry strata. Compute a standardized risk for each instrument value using those same weights, then contrast the standardized risks. For example, in the pure-confounding example, use equal weights of \(0.5\) for \(A_1\) and \(A_2\):

$$ \begin{aligned} R_{Z=1}^{\mathrm{std}}&=(0.5)(0.20)+(0.5)(0.05)=0.125,\\ R_{Z=0}^{\mathrm{std}}&=(0.5)(0.20)+(0.5)(0.05)=0.125,\\ R_{Z=1}^{\mathrm{std}}-R_{Z=0}^{\mathrm{std}}&=0. \end{aligned} $$

The standardized contrast is zero because, within each ancestry stratum, the risks are identical across instrument groups. More generally, regression adjustment or stratification can be used to condition on measured ancestry variables. The exact model and estimand depend on the outcome and analysis, but the central requirement is that comparisons are made among people with comparable measured ancestry rather than relying only on an unadjusted marginal association.

Worked Example: Checking Whether the Instrument Tracks Ancestry

In a hypothetical sample, suppose half the participants are in \(A_1\) and half in \(A_2\), with the joint proportions from the examples above. The instrument prevalence is \(0.30/0.50=0.60\) in \(A_1\) and \(0.05/0.50=0.10\) in \(A_2\). The instrument is therefore much more common in \(A_1\).

The overall instrument prevalence is \(0.30+0.05=0.35\). Among instrument carriers, the proportion in \(A_1\) is \(0.30/0.35\approx0.8571\). Among noncarriers, it is \(0.20/0.65\approx0.3077\). The difference in these ancestry proportions is approximately \(0.8571-0.3077=0.5494\), a substantial imbalance in this hypothetical sample.

This calculation is a descriptive check, not a test that establishes whether independence holds. It shows why ancestry adjustment is relevant here: the instrument groups have different ancestry composition, and the previous example showed that ancestry-specific outcome risk can convert this imbalance into a marginal instrument–outcome association.

Study-Design Safeguards and Their Limits

Adjustment is not the only response. Design choices can reduce population structure before estimation, and analysis should reflect the population from which participants were selected.

  • Restricting or matching on ancestry-related characteristics can make comparison groups more similar. Restriction may reduce generalizability and does not guarantee that meaningful structure within the selected group has been removed.
  • Recruiting across relevant groups and recording ancestry-related covariates allows the analysis to account for measured structure. Adequate overlap matters: if a variant is nearly absent in a stratum, the data may provide little information for a within-stratum comparison.
  • Using genetic principal components or related measures can capture some broad axes of genetic variation. They are proxies, not complete measurements of every social, geographic, or genetic difference relevant to the outcome.
  • Within-family designs compare genetic differences among relatives, often siblings. Because relatives share much of their ancestry and family background, these designs can reduce confounding from population structure and some family-level factors. They have their own assumptions and may be less statistically precise.

Ancestry is not a simple biological category that fully determines a person’s characteristics. Genetic variation is continuous and overlapping, while social and environmental conditions can also be patterned by geography, family history, and social processes. Adjustment should therefore be justified by the study context, carefully described, and not treated as a guarantee that all confounding has disappeared.

The causal interpretation also depends on the Mendelian randomization assumptions beyond population structure. Ancestry adjustment cannot repair a failure of the exclusion restriction, and a strong instrument does not by itself establish independence. Conversely, an association between ancestry and the instrument is not proof that the MR estimate is biased: bias depends on whether the imbalance also connects the instrument to outcome determinants and on how the analysis handles that structure.

Common Mistake

Common mistake: “The variant is inherited, so it is randomized; therefore its association with the outcome cannot be confounded.” The correction is that random segregation does not ensure that an instrument is independent of ancestry in a population sample. If allele frequencies differ across ancestry groups and outcome risks also differ, the marginal instrument–outcome association can reflect ancestry composition. Assess and address population structure as part of the instrument-independence argument.

A related error is to assume that adjustment for a few principal components proves that independence holds. Such adjustment can reduce measured structure, but adequacy depends on the variables measured, model specification, sample composition, and the causal structure of the study. The conclusion should be calibrated to the design and evidence rather than stated as certainty.

Key takeaway: Population structure threatens instrument independence when ancestry predicts both instrument frequency and outcome determinants. Compare instrument groups within or standardized to a common ancestry distribution, and use design safeguards where possible; interpret any adjustment in light of what it can and cannot measure.

Check Your Understanding

Answer each question using the definitions and reasoning in this tutorial.

  1. Explain how an allele-frequency difference between ancestry groups can create an instrument–outcome association when the instrument has no effect on the outcome within either group.
  2. In the first worked example, why is the marginal risk difference not equal to either within-ancestry risk difference?
  3. What does direct standardization change about the comparison of instrument groups, and what did it show in the pure-confounding example?
  4. Why does genetic principal-component adjustment reduce some population-structure concerns without proving that instrument independence holds?
  5. Name one way a within-family design can reduce population-structure confounding and one limitation it may have.
Clinical Biostats