Tutorials › High-Dimensional Genomic Data Analysis › Fixed-Effect Trans-Ancestry Meta-Analysis

Ancestry-Aware Association and Cross-Population Portability · Tutorial 300 of 795

Fixed-Effect Trans-Ancestry Meta-Analysis

You will derive and apply the fixed-effect inverse-variance estimator, interpret its uncertainty, and assess whether a shared trans-ancestry effect is a defensible summary.

Advanced 10 min read

What You'll Learn

  • State the shared-effect assumption required for fixed-effect trans-ancestry meta-analysis.
  • Derive the inverse-variance pooled estimate and its standard error.
  • Calculate a pooled association estimate, confidence interval, and test statistic.
  • Explain why allele-frequency differences can change precision without changing the assumed effect.
  • Recognize when a pooled fixed-effect estimate may conceal ancestry-specific effects.
  • Compare pooled analysis with reporting ancestry-stratified estimates.

Why Combine Ancestry-Stratified Estimates?

Association studies may estimate an effect separately in each ancestry group, then combine those estimates to obtain a more precise summary. This can be useful when the groups share the same underlying association effect on the chosen scale but have different sample sizes or genotype information. The fixed-effect method gives more weight to estimates with smaller standard errors.

The key qualification is substantive, not merely computational: a fixed-effect analysis assumes one common effect across the ancestry groups being combined. Differing allele frequencies can affect how precisely each group estimates that effect. They do not, by themselves, establish that the effects are equal or unequal. The earlier material on testing effect heterogeneity across ancestries addresses how to assess evidence of differences; here the focus is on combining estimates under a shared-effect assumption.

Definition: Suppose \(G\) ancestry groups provide estimates \(\hat{\beta}_g\) of an association effect, with standard errors \(SE_g\), for \(g=1,\ldots,G\). A fixed-effect trans-ancestry meta-analysis assumes that each group estimates the same effect \(\beta\), on the same outcome and effect scale. Under independent, approximately unbiased estimates with variances \(SE_g^2\), the inverse-variance pooled estimate is \[ \hat{\beta}_{FE}=\frac{\sum_{g=1}^{G}w_g\hat{\beta}_g}{\sum_{g=1}^{G}w_g}, \qquad w_g=\frac{1}{SE_g^2}. \] Its standard error is \(SE(\hat{\beta}_{FE})=1/\sqrt{\sum_g w_g}\), when the standard errors correctly represent the sampling uncertainty.

Model, Assumptions, and Derivation

For each ancestry group, suppose the association analysis estimates the change in a specified outcome associated with a one-unit increase in the same genetic predictor. The predictor coding, outcome definition, effect scale, and relevant model specification must be comparable across groups. For example, if the effect is a per-allele log-odds ratio, all estimates must refer to that same log-odds scale and the same effect allele. The estimates should also come from analyses that adequately control the relevant sources of confounding within each group.

A standard working model for the summary estimates is \(\hat{\beta}_g \sim N(\beta,SE_g^2)\), with independent estimates across groups. The normal approximation is generally justified by large-sample association analyses. The common \(\beta\) is the shared-effect assumption. Independence is reasonable for non-overlapping ancestry samples; if samples overlap, the estimates may be correlated and the ordinary inverse-variance formula is not appropriate without accounting for that covariance. In practice, \(SE_g\) is estimated, and the formula uses those estimated standard errors as plug-in measures of uncertainty.

The pooled estimator can be derived as the value of \(\beta\) that minimizes the weighted squared discrepancies between the group estimates and a common effect. With \(w_g=1/SE_g^2\), define \[ L(\beta)=\sum_{g=1}^{G}w_g(\hat{\beta}_g-\beta)^2. \] Differentiating and setting the derivative to zero gives \[ \frac{dL}{d\beta}=-2\sum_{g=1}^{G}w_g(\hat{\beta}_g-\beta)=0, \qquad \hat{\beta}_{FE}=\frac{\sum_g w_g\hat{\beta}_g}{\sum_g w_g}. \] Under the model, its variance follows from independence: \[ \operatorname{Var}(\hat{\beta}_{FE}) =\frac{\sum_g w_g^2 SE_g^2}{(\sum_g w_g)^2} =\frac{\sum_g w_g}{(\sum_g w_g)^2} =\frac{1}{\sum_g w_g}. \] This gives the stated standard error. A large-sample Wald interval is \(\hat{\beta}_{FE}\pm 1.96SE(\hat{\beta}_{FE})\), and a two-sided test of \(H_0:\beta=0\) uses \(z=\hat{\beta}_{FE}/SE(\hat{\beta}_{FE})\), compared with the standard normal distribution.

Worked Example: Pooling Two Estimates With Different Precision

Suppose two hypothetical ancestry-stratified analyses estimate the same per-allele log-odds effect for a binary outcome. Group A reports \(\hat{\beta}_1=0.20\) with \(SE_1=0.10\); Group B reports \(\hat{\beta}_2=0.10\) with \(SE_2=0.20\). Assume the estimates are independent and that a common log-odds effect is scientifically appropriate.

The inverse-variance weights are \(w_1=1/(0.10)^2=100\) and \(w_2=1/(0.20)^2=25\). Thus, \[ \hat{\beta}_{FE} =\frac{100(0.20)+25(0.10)}{100+25} =\frac{22.5}{125} =0.1800. \] The more precise estimate receives four times the weight, because its variance is one quarter as large.

The pooled standard error is \[ SE(\hat{\beta}_{FE})=\frac{1}{\sqrt{125}}=0.08944. \] The 95% Wald interval on the log-odds scale is \[ 0.1800\pm1.96(0.08944) =0.1800\pm0.1753 =(0.0047,\,0.3553). \] The test statistic is \(z=0.1800/0.08944=2.012\), giving a two-sided normal-approximation \(p\)-value of about 0.0442. Under the shared-effect model, the estimated per-allele odds ratio is \(\exp(0.1800)=1.197\): the estimated odds are about 19.7% higher per allele. The interval is an interval for the common log-odds effect under the model; it is not a probability statement about where the fixed parameter lies.

Allele Frequencies and Information

Ancestry groups can have different allele frequencies, and that can alter the precision of an association estimate. To see the mechanism in a simple setting, consider a hypothetical linear regression of a continuous outcome on genotype dosage \(G\), with an intercept, independent observations, and errors with variance \(\sigma^2\). If genotype dosage has variance \(\operatorname{Var}(G)\), the standard error of the estimated genotype coefficient is approximately \(\sigma/\sqrt{n\operatorname{Var}(G)}\). Under Hardy–Weinberg equilibrium for a biallelic genotype with effect-allele frequency \(f\), \(\operatorname{Var}(G)=2f(1-f)\). This illustration assumes no additional covariates; in an adjusted analysis, the relevant information depends on genotype variation remaining after adjustment.

Worked Example: Different Frequencies, Different Precision

Consider two hypothetical groups of \(n=1000\) individuals each. Assume a continuous outcome, residual standard deviation \(\sigma=1\), an unadjusted linear model with an intercept, independent observations, and Hardy–Weinberg equilibrium. The effect-allele frequency is \(f_A=0.30\) in Group A and \(f_B=0.05\) in Group B. For this illustration, suppose both groups have the same per-allele effect, \(\beta=0.40\) outcome units.

The genotype variances are \(2(0.30)(0.70)=0.42\) and \(2(0.05)(0.95)=0.095\). The approximate standard errors are therefore \[ SE_A=\frac{1}{\sqrt{1000(0.42)}}=0.04880, \qquad SE_B=\frac{1}{\sqrt{1000(0.095)}}=0.1026. \] The lower-frequency group has less genotype variation and consequently a less precise estimate in this simplified model.

Suppose the observed estimates are \(\hat{\beta}_A=0.40\) and \(\hat{\beta}_B=0.30\), in outcome units per allele. Their approximate weights are \(w_A=1/0.04880^2=420\) and \(w_B=1/0.1026^2=95.0\). The pooled estimate is \[ \hat{\beta}_{FE} =\frac{420(0.40)+95.0(0.30)}{420+95.0} =\frac{196.5}{515} =0.3816 \] outcome units per allele, with standard error \(1/\sqrt{515}=0.04407\) outcome units per allele. The higher-frequency group contributes more information here; its allele frequency does not make its effect intrinsically more important. The numerical standard errors rely on the simplified assumptions above, so an actual meta-analysis should use the standard errors reported by the group-specific association analyses.

Interpreting the Pooled Result and Checking Its Fit

The pooled estimate is a precision-weighted common-effect estimate, not an unweighted average of ancestry-specific effects and not an estimate of the effect in every individual. If the shared-effect assumption is credible, the result summarizes the common association on the specified scale. If true effects differ, the formula still produces a number, but that number is a precision-weighted summary of the supplied estimates; it need not correspond to a single effect that applies to each ancestry group.

Before interpreting the result, check that the group estimates target comparable quantities, that their directions use a consistent allele definition, and that their standard errors are valid. Review the estimates and intervals by group rather than relying only on the pooled result. Large discrepancies in direction or magnitude are a warning that the common-effect model may be inadequate. The earlier material on effect heterogeneity provides a formal assessment; a pooled \(p\)-value alone does not establish homogeneity.

Worked Example: A Precise Pooled Result Can Hide Disagreement

Suppose two hypothetical, independent ancestry groups report estimates on the same per-allele outcome scale: \(\hat{\beta}_1=0.30\), \(SE_1=0.05\), and \(\hat{\beta}_2=-0.10\), \(SE_2=0.10\). First, calculate the fixed-effect result mechanically: \[ w_1=1/(0.05)^2=400,\qquad w_2=1/(0.10)^2=100, \] \[ \hat{\beta}_{FE} =\frac{400(0.30)+100(-0.10)}{500} =\frac{110}{500} =0.2200, \qquad SE(\hat{\beta}_{FE})=\frac{1}{\sqrt{500}}=0.04472. \]

The pooled estimate is positive and precise under the assumed common-effect model. Yet the group-specific 95% Wald intervals are \(0.30\pm1.96(0.05)=(0.202,\,0.398)\) and \(-0.10\pm1.96(0.10)=(-0.296,\,0.096)\); they do not overlap. This is a conspicuous warning that the shared-effect assumption may not describe these data well. The calculation does not resolve the discrepancy: the pooled number should not be presented as evidence that both groups have an effect of 0.2200. Report the group estimates and investigate heterogeneity before treating the common-effect summary as the scientific target.

Choosing a Summary and Recognizing Limitations

Fixed-effect pooling is most useful when a common association effect is a defensible target and the group estimates are independent, comparable, and adequately calibrated. Its principal advantage is efficiency: it combines information while giving more influence to estimates with lower sampling variance. Its principal limitation follows from the same weighting rule. Groups with more information can dominate the pooled estimate, so it need not describe the effect in a smaller group or represent an equal-weight average across ancestries.

Reporting ancestry-stratified estimates alongside a pooled estimate makes the weighting and any disagreement visible. If the scientific goal is to describe group-specific effects, the stratified estimates may be more informative than a single summary. If there is credible between-group variation, a model that allows such variation is an alternative, but it answers a different question and requires additional assumptions. A fixed-effect analysis should not be used as a substitute for that model simply because it gives a convenient single estimate.

The method also depends on comparable effect scales. A shared numerical coefficient has a coherent interpretation only when the outcome, predictor, and coefficient scale are aligned. Differences in allele frequency primarily influence information in the idealized model above; other population differences can also alter the estimated association target. Thus, a similar-looking coefficient does not by itself prove that the same underlying biological or causal quantity has been estimated. Keep the interpretation at the level supported by the association model.

Common Pitfalls

  • Assuming that a pooled estimate proves a common effect. The formula returns a weighted average even when group estimates disagree. Check the ancestry-stratified results and assess whether the shared-effect assumption is plausible.
  • Using sample size alone as the weight. The standard fixed-effect weight is \(1/SE_g^2\). Sample size is only one determinant of precision; genotype variation and other features of the analysis matter too.
  • Interpreting different precision as different effect. A group with a larger standard error may simply contain less information. Precision and effect magnitude are distinct quantities.
  • Ignoring dependence between estimates. The standard error formula assumes independence. Overlapping samples can induce covariance, which must be accounted for rather than treated as independent information.
  • Reporting only the pooled result. A single precise number can conceal meaningful ancestry-specific disagreement. Include the stratum-level estimates when interpreting trans-ancestry evidence.
Key takeaway: Fixed-effect trans-ancestry meta-analysis combines comparable, independent estimates by inverse variance, and its standard error follows from that weighting. The pooled result is interpretable as one shared effect only when the shared-effect assumption is reasonable; differing allele frequencies can change precision without establishing effect homogeneity.

Check Your Understanding

Work each problem, then compare your reasoning with the solutions.

  1. Two independent groups estimate a common effect: \(\hat{\beta}_1=0.12\), \(SE_1=0.06\), and \(\hat{\beta}_2=0.20\), \(SE_2=0.12\). Calculate the fixed-effect estimate and standard error.
  2. In one sentence, state the shared-effect assumption and name two other conditions needed for the ordinary inverse-variance formula.
  3. Explain why a lower allele frequency can increase the standard error without implying a smaller per-allele effect.
  4. A pooled estimate is precise, but the group-specific estimates have opposite signs. What should be checked or reported before interpreting the pooled estimate?

Solutions. (1) The weights are \(1/0.06^2=277.8\) and \(1/0.12^2=69.44\). Thus \(\hat{\beta}_{FE}=[277.8(0.12)+69.44(0.20)]/347.2=0.1360\), and \(SE=1/\sqrt{347.2}=0.05367\). (2) The assumption is that each group estimates the same effect on the same scale; the estimates should also be independent and approximately unbiased with correctly estimated standard errors. (3) Lower allele frequency can reduce genotype variation and therefore information, increasing uncertainty; that concerns precision, not the assumed coefficient itself. (4) Examine and report the ancestry-specific estimates and intervals, check that they target comparable effects and are consistently coded, and assess heterogeneity. Do not treat the pooled result alone as evidence that the groups share one effect.

Clinical Biostats