Association Evidence Is Not the Same as Variant Resolution
Combining genetic association results across ancestry groups can increase precision and make a locus easier to detect. But finding a strong association is not the same as identifying which variant at that locus is causal. A correlated, noncausal variant can show a strong association because it is in linkage disequilibrium (LD) with a causal variant. The pattern of that correlation can differ across ancestry groups.
The key distinction is what gets combined. A conventional meta-analysis often combines an estimate and standard error separately for each variant. This can be useful for testing association, but a set of pooled statistics does not, by itself, describe how each ancestry’s LD generated its pattern of associations. Joint trans-ancestry fine-mapping retains ancestry-specific association patterns and their LD information in the likelihood.
These approaches are not opposites in every implementation. Meta-analysis can be one step in a broader analysis, and a meta-analysis result may be followed by fine-mapping that uses suitable LD information. The important distinction is whether ancestry-specific patterns and LD are retained for the fine-mapping calculation, rather than reducing the evidence to pooled association statistics alone.
What a Fixed-Effect Meta-Analysis Combines
Suppose ancestry group \(a\) provides an estimate \(\hat{\beta}_a\) for the effect of the same harmonized effect allele on the same trait, with standard error \(SE_a\). Under a fixed-effect model, the true effect is assumed common across groups, and the estimates are treated as independent with approximately normal sampling distributions. The inverse-variance estimate and its standard error are:
The corresponding Wald statistic is \(Z_{\text{meta}}=\hat{\beta}_{\text{meta}}/SE_{\text{meta}}\). It summarizes evidence for association at that variant under the model. It does not say whether that variant is causal, or whether a correlated variant generated the association.
This calculation requires estimates to refer to the same allele and trait scale, and the groups to be independent for the stated variance formula. The fixed-effect interpretation additionally assumes a common underlying effect. If effects genuinely differ, a pooled estimate can still be computed, but it describes a weighted average rather than a common effect. Harmonization and heterogeneity diagnostics are therefore relevant to interpretation, not optional details.
Worked Example: Pooling Two Variant-Level Estimates
Consider hypothetical results for one harmonized variant, with effects measured in trait standard deviations per effect allele. An ancestry group A has \(\hat{\beta}_A=0.10\) and \(SE_A=0.025\); an independent group B has \(\hat{\beta}_B=0.10\) and \(SE_B=0.050\). Assume for this calculation that the fixed-effect model is appropriate.
The weights are \(w_A=1/0.025^2=1600\) and \(w_B=1/0.050^2=400\). Thus,
The pooled standard error is \(SE_{\text{meta}}=1/\sqrt{2000}=0.02236\) trait standard deviations per effect allele. The statistic is \(Z_{\text{meta}}=0.100/0.02236=4.472\). The pooled estimate is more precise than either group’s estimate, because both groups contribute information. This calculation establishes strong association evidence at this variant under the stated assumptions; it does not establish that this is the causal variant.
Why LD Changes the Fine-Mapping Question
For a region containing \(p\) variants, let \(z_a\) denote the vector of ancestry-specific association statistics and \(R_a\) the corresponding LD correlation matrix. In a common large-sample approximation with standardized genotypes and a quantitative trait, the summary-statistic distribution can be written as:
Here \(n_a\) is the sample size for group \(a\), and \(\gamma_a\) is the vector of standardized genetic effects. This approximation assumes the summary statistics and LD reference describe comparable samples and allele coding; actual implementations may use other scales or likelihoods. Its essential message is that LD affects both the expected association pattern and the covariance among statistics. If variant \(j\) is the only causal variant, the mean pattern is proportional to column \(j\) of \(R_a\), not simply a spike at \(j\).
With independent ancestry groups, the joint likelihood combines their contributions while preserving their different \(R_a\) matrices. In a simple shared-effect illustration, \(\gamma_a=\gamma\) for each group. The likelihood is proportional to the product of the ancestry-specific distributions, so the observed pattern in each group can be evaluated against the LD pattern expected under each causal candidate. Other models can allow group-specific effects; the central point here is that the LD structure can remain ancestry-specific.
Worked Example: A Shared Variant Produces Different Association Patterns
Suppose two hypothetical ancestry groups each have \(n_a=100\) and the same standardized effect \(\gamma=0.40\) at causal variant 1. There are two variants. Their LD matrices are \(R_A=\begin{pmatrix}1&0.90\\0.90&1\end{pmatrix}\) and \(R_B=\begin{pmatrix}1&0.10\\0.10&1\end{pmatrix}\). Both matrices are valid correlation matrices. Under the approximation above, \(\sqrt{n_a}\gamma=10(0.40)=4\). If variant 1 is causal, the expected association-statistic vectors are:
These are expected patterns, not claims that real observed statistics must equal their means. In group A, the noncausal variant also has a strong expected statistic because it is highly correlated with the causal variant. In group B, its expected statistic is much smaller. If the two groups instead had identical LD, the association pattern would be more alike across them.
For comparison, if variant 2 were causal with the same effect, the expected vectors would be \((3.6,4)\) in group A and \((0.4,4)\) in group B: the relevant LD columns have switched. A joint analysis can compare each group’s observed vector to these candidate-specific patterns. The information is not just that one or both variants have small p-values; it is how their relative associations align with the LD in each group.
What Pooling Retains and What It Can Lose
A pooled statistic can retain substantial information about whether a variant is associated. It can also preserve some ability to distinguish variants, especially when the pooled association patterns differ clearly and the fine-mapping procedure uses an appropriate LD model. But pooling the statistics can discard which ancestry-specific pattern contributed to the result.
For illustration, with two equally precise independent groups, a common-effect meta-analysis of \(Z\)-statistics at one variant is \(Z_{\text{meta}}=(Z_A+Z_B)/\sqrt{2}\). The scaling gives a statistic with unit variance under the null. The operation is a sum: it does not record the two input values separately. For instance, ancestry-specific pairs \((4,0)\) and \((2,2)\) produce the same pooled value, \(4/\sqrt{2}\), even though the distribution of evidence across groups differs. If their LD structures differ, those patterns may have different implications for candidate causal variants.
This is a direct information-loss argument. If a pooled statistic is \(T=Z_A+Z_B\), then for any value \(d\), the pair \((Z_A+d,Z_B-d)\) has the same \(T\). The map from two ancestry-specific values to one pooled value is therefore not one-to-one. A joint likelihood can use the separate values and their ancestry-specific covariance structures; the pooled value alone cannot reconstruct them.
Worked Example: A Strong Pooled Signal With Little Resolution
Consider hypothetical data for two variants in two equally sized groups. In group A their LD correlation is \(0.98\), and in group B it is \(0.96\). Suppose the observed ancestry-specific statistics are \(z_A=(5.0,4.9)\) and \(z_B=(5.0,4.8)\). These values are compatible with a strong association at either candidate when the variants are so highly correlated.
For each variant, the equal-information pooled statistic is the sum across groups divided by \(\sqrt{2}\). For variant 1 it is \((5.0+5.0)/\sqrt{2}=7.071\). For variant 2 it is \((4.9+4.8)/\sqrt{2}=6.859\). Both are strong association signals. Yet their pooled statistics are close, and the ancestry-specific LD values indicate that either variant can be a strong proxy for the other.
This example does not prove that a fine-mapping analysis must fail: actual resolution depends on the observed data, sample sizes, LD estimates, and model. It shows why small p-values alone are insufficient. To assess which variant is better supported as causal, an analysis needs candidate-specific likelihoods or posterior evidence that account for LD, not merely a declaration that the most significant variant is causal.
A Short Derivation of the Joint Versus Pooled Distinction
Under the normal approximation, the ancestry-specific log likelihood for a proposed effect vector \(\gamma_a\), apart from constants, is:
For independent groups, the joint log likelihood is \(\ell_{\text{joint}}=\sum_a\ell_a\). Thus each group contributes a residual measured against its own expected association pattern and covariance. If the candidate causal variant changes, the relevant column of each \(R_a\) changes, potentially by a different amount in each group.
By contrast, a pooled statistic combines observations before this ancestry-specific comparison. Because summing is many-to-one, the pooled value does not generally determine all the \(z_a\) vectors. When LD differs, those missing vectors can be informative about which candidate better explains the evidence. This justification does not imply that every joint analysis is more precise or that every meta-analysis loses all useful information. It explains why the two procedures answer distinct questions and why a pooled association statistic alone is not a substitute for a joint fine-mapping likelihood.
When the Distinction Matters
The distinction is most important when variants are correlated, LD differs appreciably across ancestry groups, and the analysis aims to prioritize causal variants rather than simply test association. In that setting, a large group may dominate an inverse-variance pooled result, even if a smaller group’s distinct LD pattern is valuable for separating candidates. Conversely, if LD is similar across groups or the smaller group has little information, retaining ancestry-specific data may not materially change resolution.
Before interpreting results, check that variant alleles and trait scales are aligned, that the LD matrices match the relevant ancestry groups and variant set, and that the assumed effect model is suitable. Incomplete variant coverage can also change which candidates can be compared; a variant absent or poorly imputed in one group cannot contribute the same evidence there. A technically correct pooled association analysis may therefore still be inadequate for the causal-resolution question.
Check Your Understanding
Use the definitions and examples above to answer these questions.
- What fixed-effect assumption underlies inverse-variance pooling, and what does the pooled estimate represent if effects differ across groups?
- Why can a noncausal variant have a strong association statistic in one ancestry group?
- In the two-group example with LD correlations 0.90 and 0.10, what feature of the ancestry-specific patterns helps distinguish the two candidate variants?
- Why does a pooled statistic generally fail to recover the separate ancestry-specific statistics?
- What additional information must a fine-mapping analysis use before a strong association can support causal-variant resolution?