Tutorials › Statistical Genetics › Reference Panel Mismatch and Ancestry Bias

Ancestry and Population Structure · Tutorial 152 of 989

Reference Panel Mismatch and Ancestry Bias

Learn to diagnose ancestry estimates that depend on an unrepresentative reference panel, and to assess the effects of panel choice with explicit comparisons and sensitivity checks.

Intermediate 10 min read

What You'll Learn

  • Define reference-panel mismatch and distinguish it from uncertainty caused by small panels.
  • Derive how allele-frequency mismatch can bias a two-source ancestry estimate.
  • Identify signs that an ancestry source is omitted from the reference panel.
  • Calculate how mismatched allele frequencies can alter a local ancestry call.
  • Compare estimates across plausible panels and interpret their sensitivity.

Why Reference-Panel Match Matters

An ancestry estimate is conditional on the reference populations used to construct it. The panel is not a neutral backdrop: its allele frequencies help define what the model means by each ancestry component. If a reference group is a poor proxy for a source population, the model may assign too much or too little ancestry to that component. If a source is omitted entirely, the model may explain its genetic variation using the available components, even when none is a good fit.

This is a problem of model specification, not merely sampling variability. A larger panel can reduce uncertainty about the allele frequencies of the sampled reference groups, but it cannot make those groups representative of a different source population. The practical question is therefore not only whether an estimate is numerically precise, but also whether it is stable under reasonable changes to the panel.

Definition: Reference-panel mismatch occurs when the allele frequencies or haplotypes represented by a panel do not adequately represent the ancestry sources relevant to the target sample. Omission is the limiting case in which a relevant source is absent from the panel. An ancestry estimate is panel-sensitive when it changes materially under plausible alternative panels or panel definitions.

For a simple illustration, suppose a target population has ancestry proportion \(q\) from source A and \(1-q\) from source B. At marker \(i\), let \(f_{iA}\) and \(f_{iB}\) be the allele frequencies in the two true sources. Under a two-source mixture model, the expected target allele frequency is

$$ p_i=qf_{iA}+(1-q)f_{iB}. $$

This relation assumes that the target is represented by the stated two sources, that the allele frequencies refer to the same marker and allele, and that the mixture proportion applies across the markers being considered. In real estimation, many markers contribute, genotypes are sampled, and the model may incorporate additional structure. The equation isolates the basic issue: the model substitutes panel frequencies for the source frequencies, and a mismatch changes the fitted proportions.

How Mismatch Creates Bias

At one marker, suppose the observed target frequency is \(p_i\), while the panel provides frequencies \(f^{*}_{iA}\) and \(f^{*}_{iB}\). If the fitted two-source model matches that frequency exactly, its estimate for the A proportion is

$$ \hat q_i=\frac{p_i-f^{*}_{iB}}{f^{*}_{iA}-f^{*}_{iB}}, \qquad f^{*}_{iA}\ne f^{*}_{iB}. $$

The denominator condition matters: if the panel frequencies are equal at this marker, the marker contains no information for distinguishing the two sources. Even when they differ, a small difference makes the estimate highly sensitive to small frequency errors. Estimates based on many markers combine information across loci, but the same principle remains: the fitted ancestry proportions depend on the reference frequencies and on how informative the markers are.

Worked Example: A Shifted Reference Frequency Inflates an Ancestry Estimate

Suppose a hypothetical target has true ancestry proportions \(q=0.60\) from source A and \(0.40\) from source B. At one marker, the true source allele frequencies are \(f_A=0.20\) and \(f_B=0.80\). The expected target frequency is

$$ p=(0.60)(0.20)+(0.40)(0.80)=0.12+0.32=0.44. $$

Now suppose the panel represents A with frequency \(f_A^{*}=0.30\), while its B frequency remains \(f_B^{*}=0.80\). Applying the two-source frequency equation gives

$$ \hat q=\frac{0.44-0.80}{0.30-0.80} =\frac{-0.36}{-0.50}=0.72. $$

The panel-based estimate is 0.72 for source A, although the hypothetical generating proportion is 0.60. Here the panel frequency for A has moved toward the frequency for B, making the target frequency appear to require more A to match it. This calculation isolates systematic mismatch; it does not include finite-sample uncertainty.

There is no universal direction to this bias. It depends on which allele frequencies are shifted, how those shifts relate to the target, the markers included, and the fitting procedure. Some loci can push an estimate upward and others downward. A model that reports one genome-wide proportion can therefore conceal substantial disagreement among marker-level signals.

Omitted Sources and Poor Fit

If the target contains ancestry from a source absent from the panel, a model with only the available sources must allocate that variation somewhere. The resulting proportions may be numerically well-defined while still failing to describe the target adequately. In a simple frequency model, an omitted source can make the target frequency fall outside the range attainable by mixtures of the represented sources. With multiple markers, the target may be compatible with the panel at some loci but not others, so no single mixture proportion fits them all.

Worked Example: An Omitted Source Produces Inconsistent Fits

Consider a hypothetical target with true proportions 0.50 from A, 0.30 from B, and 0.20 from C. A two-source panel includes A and B but omits C. At marker 1, the source allele frequencies are \(f_A=0.10\), \(f_B=0.70\), and \(f_C=0.99\). The target frequency is

$$ p_1=(0.50)(0.10)+(0.30)(0.70)+(0.20)(0.99) =0.050+0.210+0.198=0.458. $$

If the two-source model tries to match marker 1, its implied A proportion is

$$ \hat q_1=\frac{0.458-0.70}{0.10-0.70} =\frac{-0.242}{-0.60}\approx 0.403. $$

At marker 2, let the frequencies be 0.20 in A, 0.40 in B, and 0.99 in C. The target frequency is

$$ p_2=(0.50)(0.20)+(0.30)(0.40)+(0.20)(0.99) =0.100+0.120+0.198=0.418. $$

A mixture of only A and B must have a frequency between 0.20 and 0.40 at this marker. Since \(0.418>0.40\), no two-source proportion can reproduce the target frequency there. The first marker yields a feasible but misleading two-source fit; the second reveals a lack of fit. In a many-marker analysis, analogous discrepancies may appear as residual patterns or unstable proportions, rather than as one obvious impossible frequency.

A poor fit does not by itself identify which population is missing or prove that a particular alternative panel is correct. It is evidence that the stated panel and model may not adequately represent the target. Other explanations, such as marker quality issues or frequency-estimation error, should also be considered. The useful response is to examine whether the mismatch signal persists across markers and under reasonable panel alternatives.

Local Ancestry Calls Can Also Shift

Global ancestry proportions summarize contributions across the genome. Local ancestry methods assign or estimate ancestry along genomic segments, using marker information and, in many methods, information about neighboring positions. They rely on reference frequencies or haplotypes to distinguish sources. If those reference patterns are mismatched, the evidence at a segment can be misweighted: a segment may be assigned to the wrong represented source, or the evidence may become less decisive.

The following single-marker calculation illustrates the emission evidence, not a complete local ancestry algorithm. Assume one observed allele is modeled as Bernoulli within each source, and give the two sources equal prior probability at the position. Under these assumptions, the posterior probability for A is proportional to the A allele likelihood. Actual local ancestry methods may combine many markers and account for transitions between neighboring positions.

Worked Example: A Mismatched Panel Flips a Local Call

Suppose a hypothetical haploid allele observation is \(x=1\). The true source frequencies are 0.20 in A and 0.80 in B. With equal prior probabilities, the likelihoods for this observation are 0.20 under A and 0.80 under B, so the posterior probability of A is

$$ P(A\mid x=1)=\frac{0.20}{0.20+0.80}=0.20. $$

Using a 0.50 posterior threshold, the call is B. Now suppose the panel uses a mismatched A frequency of 0.90, while the B frequency remains 0.80. The panel-based posterior becomes

$$ P^{*}(A\mid x=1)=\frac{0.90}{0.90+0.80} =\frac{0.90}{1.70}\approx 0.529. $$

At the same threshold, the call is now A. The example shows how a changed reference frequency can alter the direction of the local evidence. It does not imply that every local call is determined by one marker; neighboring markers, transition assumptions, and uncertainty in estimated frequencies can all affect a real segment-level result.

Assessing Panel Sensitivity

A sensitivity analysis asks whether the substantive conclusion persists when plausible choices about the reference panel change. The alternatives should be scientifically defensible, not selected only because they produce a preferred estimate. Depending on the data and question, useful comparisons include a more geographically or genetically relevant panel, a panel with a relevant group added, a panel with a questionable proxy removed, and panels formed by reasonable subsets of reference samples. A comparison should keep other analysis choices as consistent as possible so that differences can be attributed to panel choice.

For global ancestry, report the estimates from each panel and the change in each component in its natural units, such as percentage points. For local ancestry, compare segment or position-level posterior probabilities, call changes, and regions where the evidence changes. A change in a hard call may be more informative when accompanied by the corresponding probabilities: a shift from 0.51 to 0.49 changes a thresholded label but indicates weak evidence in both analyses.

1
Specify plausible panels.
State which reference groups are included, excluded, or substituted, and why each choice is defensible for the target population.
2
Hold the analysis fixed where possible.
Use the same target data, markers, and estimation settings across comparisons so panel changes are interpretable.
3
Compare estimates and fit.
Examine global proportions, local probabilities or calls, and any available residual or model-fit diagnostics. Look for marker- or segment-specific disagreement as well as changes in summaries.
4
Report dependence, not just a preferred result.
Describe conclusions that remain stable, conclusions that vary, and the range of estimates across the defensible panels considered.

A stable estimate across tested panels is reassuring about sensitivity to those particular choices, but it does not prove that the panel represents every relevant source. Conversely, a changing estimate is not automatically evidence that one panel is wrong; it shows that the result depends on panel specification and deserves qualification. Panel sensitivity is distinct from uncertainty caused by small reference samples: the latter concerns imprecision in estimated frequencies, while mismatch concerns whether the represented frequencies are appropriate in the first place. Both can occur together.

Common Mistake and Practical Interpretation

Wrong reasoning: “The ancestry estimate has a small standard error, so the reference panel must be appropriate.” A small standard error can describe precision conditional on the fitted model and its inputs. It does not establish that the model includes the relevant sources or that its reference frequencies match them.

Correction: Treat precision and representativeness as separate questions. Examine alternative panels and fit diagnostics, and state clearly that ancestry proportions are estimates relative to the specified reference framework. They are not direct measurements of a person's complete ancestry, and a change in a component is not evidence of a biological cause or of a change in identity.

These checks matter because ancestry estimates may be used to describe samples, compare groups, or inform downstream genetic analyses. When panel mismatch is ignored, apparent differences between people or groups can partly reflect the chosen reference rather than the target samples alone. Transparent panel comparisons make that dependence visible and help prevent overconfident interpretation.

Key takeaway: Reference-panel mismatch can systematically distort global proportions and local ancestry evidence; omission can also produce poor fit. Compare defensible alternative panels, examine fit across markers or segments, and report how sensitive the conclusions are to panel choice.

Check Your Understanding

Answer these questions using the equations and reasoning in this tutorial.

  1. In the two-source frequency equation, why can a marker with nearly equal panel allele frequencies provide an unstable ancestry estimate?
  2. A target frequency is 0.60, and the panel frequencies are 0.20 for A and 0.80 for B. What A proportion matches this marker under the two-source model?
  3. Why can adding many individuals to an unrepresentative reference group fail to remove mismatch bias?
  4. In the omitted-source example, what does it mean that the target frequency at marker 2 exceeds the range attainable by mixtures of A and B?
  5. Why should local ancestry probabilities be compared alongside thresholded calls in a panel sensitivity analysis?
Clinical Biostats