Tutorials › Statistical Genetics › Horizontal Pleiotropy And The Exclusion Restriction

Mendelian Randomization Foundations · Tutorial 843 of 989

Horizontal Pleiotropy And The Exclusion Restriction

See why a valid genetic instrument must affect the outcome only through the exposure, and quantify the Wald-ratio bias when that condition fails.

Advanced 10 min read

What You'll Learn

  • State the relevance, independence, and exclusion requirements for a genetic instrument
  • Distinguish an effect on the outcome outside the exposure pathway from an effect carried through the exposure
  • Derive the Wald ratio when a variant has a direct outcome effect
  • Calculate the size and direction of direct-effect bias in hypothetical examples
  • Explain why residual association with the variant must be ruled out in the association derivation
  • Identify what a single variant’s Wald ratio cannot establish about instrument validity

Why the Exclusion Restriction Matters

Mendelian randomization (MR) uses genetic variation as an instrumental variable to estimate the causal effect of an exposure \(X\) on an outcome \(Y\). The central appeal is that a genetic variant \(G\) can provide variation in \(X\) that is less vulnerable to some sources of confounding than an ordinary comparison of people with different exposure levels. But this interpretation depends on what the variant does to the outcome. If \(G\) can affect \(Y\) by a route that does not pass through \(X\), its outcome association is not solely evidence about the causal effect of \(X\).

That failure is called horizontal pleiotropy when the variant affects the outcome through a pathway outside the exposure pathway. It violates the exclusion restriction. The consequence is concrete: the variant’s outcome association includes both the effect transmitted through \(X\) and an extra contribution, so the Wald ratio need not equal the causal effect.

Definition: A genetic variant is a valid instrument for the effect of \(X\) on \(Y\) when it is relevant to \(X\), independent of causes that would undermine the instrument–outcome comparison, and satisfies the exclusion restriction: it affects \(Y\) only through \(X\). A horizontal pleiotropic effect is an effect of the variant on \(Y\) that does not operate through \(X\), and therefore violates the exclusion restriction.

A Model for the Variant Associations

Use a linear, additive model with a single genetic variant. Let \(G\) denote the variant, often coded as the number of copies of a specified allele. Let \(\gamma\) be its effect on the exposure, \(\beta\) the causal effect of the exposure on the outcome, and \(\alpha\) the variant’s direct effect on the outcome that is not mediated by \(X\). The model is

$$ X=\gamma G+U_X,\qquad Y=\beta X+\alpha G+U_Y. $$

Here \(U_X\) and \(U_Y\) collect other determinants of the exposure and outcome. They may be associated with one another; for example, causes of \(X\) and \(Y\) can confound their observational association. For the marginal-association derivation below, however, assume explicitly that

$$ \operatorname{Cov}(G,U_X)=0,\qquad \operatorname{Cov}(G,U_Y)=0. $$

These conditions ensure that the residual terms do not themselves contribute to the marginal associations of \(G\) with \(X\) or \(Y\). They are stronger and more precise than merely saying that \(G\) is independent of exposure–outcome confounders: independence from those confounders alone does not guarantee that \(G\) is unrelated to every component of \(U_X\) and \(U_Y\). In the usual instrumental-variable interpretation, the relevant independence condition must make the instrument exogenous with respect to the causes represented in the residuals. The model also assumes the stated linear effects and no interaction terms.

Define \(\gamma_G\) as the population slope from regressing \(X\) on \(G\), and \(\Gamma_G\) as the population slope from regressing \(Y\) on \(G\). For a nonconstant \(G\), each slope is its covariance with the outcome divided by \(\operatorname{Var}(G)\). Thus, under the residual-orthogonality conditions,

$$ \gamma_G =\frac{\operatorname{Cov}(G,X)}{\operatorname{Var}(G)} =\gamma, \qquad \Gamma_G =\frac{\operatorname{Cov}(G,Y)}{\operatorname{Var}(G)} =\beta\gamma+\alpha. $$

The relevance condition requires \(\gamma\ne 0\): the variant must be associated with the exposure. The exclusion restriction corresponds to \(\alpha=0\) in this model. If \(\alpha\ne 0\), the variant has an outcome effect beyond the route through \(X\). The independence and residual-orthogonality conditions are distinct from exclusion: a variant can have no residual association and still violate exclusion through a nonzero \(\alpha\).

How Horizontal Pleiotropy Changes the Wald Ratio

For one relevant variant, the Wald ratio divides its outcome association by its exposure association. In the population model above, the ratio is

$$ \frac{\Gamma_G}{\gamma_G} =\frac{\beta\gamma+\alpha}{\gamma} =\beta+\frac{\alpha}{\gamma}. $$

The term \(\alpha/\gamma\) is the departure from the causal effect \(\beta\). Its sign depends on both the direct outcome effect and the exposure association. For a fixed \(\alpha\), a smaller absolute value of \(\gamma\) makes the departure larger. A stronger association with the exposure does not prove that \(\alpha=0\); it only reduces this particular ratio’s sensitivity to a fixed direct effect.

Key result: Under the linear model and the stated residual-orthogonality conditions, the Wald ratio equals the causal effect \(\beta\) when the variant is relevant and satisfies exclusion (\(\gamma\ne 0\) and \(\alpha=0\)). With a horizontal direct effect, it equals \(\beta+\alpha/\gamma\), not \(\beta\). This is a statement about population associations; sampling error can add further uncertainty to an estimated ratio.

Worked Example: A Direct Effect Shifts the Ratio

Suppose hypothetical population associations are generated by a variant whose exposure effect is \(\gamma=0.40\) exposure units per allele. Let the causal effect be \(\beta=2.0\) outcome units per exposure unit, and let the variant also have a direct outcome effect \(\alpha=1.5\) outcome units per allele. Assume the model’s residuals are uncorrelated with \(G\), as required for the association equations above.

The outcome association combines the mediated contribution and the direct contribution:

$$ \Gamma_G=\beta\gamma+\alpha =(2.0)(0.40)+1.5 =0.80+1.5 =2.30 \quad\text{outcome units per allele}. $$

The Wald ratio is therefore

$$ \frac{\Gamma_G}{\gamma_G} =\frac{2.30}{0.40} =5.75 \quad\text{outcome units per exposure unit}. $$

The causal effect specified in the model is \(2.0\) outcome units per exposure unit. The ratio’s difference from that value is \(5.75-2.0=3.75\), which also follows from \(\alpha/\gamma=1.5/0.40=3.75\). The ratio is biased upward relative to \(\beta\) in this example because both \(\alpha\) and \(\gamma\) are positive.

Worked Example: A Relevant Variant With No Direct Effect

Consider a second hypothetical model in which \(\gamma=0.25\) exposure units per allele, \(\beta=-0.80\) outcome units per exposure unit, and \(\alpha=0\) outcome units per allele. Assume again that \(G\) is uncorrelated with both model residuals. The variant is relevant because its exposure association is nonzero, and it satisfies exclusion in this model because its direct effect is zero.

The variant–outcome association is

$$ \Gamma_G=\beta\gamma+\alpha =(-0.80)(0.25)+0 =-0.20 \quad\text{outcome units per allele}. $$

Dividing by the exposure association gives

$$ \frac{\Gamma_G}{\gamma_G} =\frac{-0.20}{0.25} =-0.80 \quad\text{outcome units per exposure unit}. $$

In this model, the population Wald ratio equals the specified causal effect. This agreement follows from the assumptions, including exclusion; it is not guaranteed just because a calculated ratio is numerically plausible. With real data, estimates of the two associations have sampling error, so an estimated ratio will not generally equal the underlying \(\beta\) exactly.

Why Residual Orthogonality Must Be Stated

The expression \(\Gamma_G=\beta\gamma+\alpha\) does not follow from the structural equations alone. It also requires that neither residual contribute to the marginal association with \(G\). Substituting the equation for \(X\) into the equation for \(Y\) makes the dependence clear:

$$ Y=(\beta\gamma+\alpha)G+\beta U_X+U_Y. $$

Taking covariance with \(G\) and dividing by its variance gives the general association expressions

$$ \gamma_G=\gamma+\frac{\operatorname{Cov}(G,U_X)}{\operatorname{Var}(G)}, \qquad \Gamma_G=\beta\gamma+\alpha+ \frac{\beta\operatorname{Cov}(G,U_X)+\operatorname{Cov}(G,U_Y)} {\operatorname{Var}(G)}. $$

Consequently, if either residual is associated with the variant, the observed associations include extra terms. In that case, one cannot simply attribute the whole outcome association to \(\beta\gamma+\alpha\), or interpret the ratio’s departure from \(\beta\) as \(\alpha/\gamma\). This is why the assumptions supporting a formula matter as much as the formula itself.

Worked Example: Residual Association Adds Another Term

For a hypothetical illustration, take \(\operatorname{Var}(G)=1\), structural exposure coefficient \(\gamma=0.30\) exposure units per allele, \(\beta=2.0\) outcome units per exposure unit, and direct effect \(\alpha=0.50\) outcome units per allele. Suppose \(\operatorname{Cov}(G,U_X)=0.20\) and \(\operatorname{Cov}(G,U_Y)=0.10\), using the corresponding covariance units. These values deliberately violate the residual-orthogonality condition.

The observed exposure association is not \(\gamma\); it is

$$ \gamma_G=0.30+\frac{0.20}{1}=0.50 \quad\text{exposure units per allele}. $$

The outcome association is

$$ \Gamma_G=(2.0)(0.30)+0.50+ \frac{(2.0)(0.20)+0.10}{1} =0.60+0.50+0.50 =1.60 \quad\text{outcome units per allele}. $$

The observed Wald ratio is \(1.60/0.50=3.20\) outcome units per exposure unit. It is not obtained by treating the observed exposure association as the structural \(\gamma=0.30\) and applying \(\beta+\alpha/\gamma\), which would give \(2.0+0.50/0.30\approx3.67\). The residual-association terms change both marginal associations. This example demonstrates why the residual conditions must be checked before using the simpler direct-effect bias formula.

Interpretation and a Common Pitfall

The exclusion restriction is a claim about the causal routes from the variant to the outcome, not a property that can be established by the Wald ratio’s arithmetic. In the model, a pathway from \(G\) to \(Y\) that bypasses \(X\) is represented by \(\alpha G\). If that effect is present, the outcome association contains information about both the exposure-mediated effect and the direct effect. A single variant’s two associations provide a ratio, but do not by themselves separate \(\beta\) from \(\alpha\).

A common mistake is to say, “The variant is strongly associated with the exposure, so the Wald ratio is a causal estimate.” Relevance is necessary, but it does not imply exclusion. Another mistake is to attribute every difference between an estimated ratio and a proposed causal effect to horizontal pleiotropy. Sampling error, residual association, or a mismatch between the model and the data could also matter. The algebra identifies the consequence of \(\alpha\ne0\) only under the stated model and residual conditions.

Takeaway: A valid instrument must be relevant, satisfy the independence condition, and affect the outcome only through the exposure. Under the linear model with residuals uncorrelated with the variant, horizontal pleiotropy changes the Wald ratio by \(\alpha/\gamma\). The ratio alone cannot verify that the direct effect is zero.

Check Your Understanding

Answer each question using the model and assumptions stated above.

  1. In \(Y=\beta X+\alpha G+U_Y\), what does \(\alpha\ne0\) mean for the exclusion restriction?
  2. If \(\gamma=0.5\), \(\beta=3\), and \(\alpha=-1\), calculate \(\Gamma_G\) and the Wald ratio. Give the ratio’s units relative to the units of \(X\) and \(Y\).
  3. Why is relevance alone insufficient to make a variant a valid instrument?
  4. Which residual covariance conditions are needed to derive \(\gamma_G=\gamma\) and \(\Gamma_G=\beta\gamma+\alpha\) from the displayed model?
  5. Suppose \(\alpha=0\) and the residual conditions hold. What does the population Wald ratio equal, provided \(\gamma\ne0\)?
Clinical Biostats