Tutorials › Statistical Genetics › Sample Overlap In Cross-Trait Summary Analyses

Genetic Correlation Across Traits · Tutorial 568 of 989

Sample Overlap In Cross-Trait Summary Analyses

Calculate the sampling covariance induced by overlapping participants and interpret it alongside the genetic signal and intercept diagnostics in cross-trait LDSC.

Advanced 10 min read

What You'll Learn

  • Derive the approximate covariance between two association statistics from their shared participants
  • Calculate the overlap contribution from sample sizes and the traits’ correlation
  • Interpret the sign and scale of the overlap contribution in context
  • Separate an approximately LD-independent intercept contribution from LD-dependent genetic covariance
  • Recognize why an LDSC intercept is a diagnostic, not a definitive count of overlapping participants

Why Sample Overlap Matters

Cross-trait summary analyses compare association statistics for two traits across many genetic variants. If the studies include some of the same participants, the two statistics can be correlated even when the variant has no causal effect on either trait. The shared participants contribute to both statistics, and correlated trait values in those people create a shared sampling fluctuation.

This overlap-induced covariance is distinct from covariance due to shared genetic effects. In cross-trait LD score regression (LDSC), the central diagnostic is how the expected product of the two association statistics varies with LD score: shared polygenic effects contribute an LD-dependent component, whereas overlap often contributes to the approximately LD-independent intercept. The distinction is useful, but not absolute: the intercept can also reflect other sources of confounding, and overlap information alone does not establish genetic sharing.

Definition: Let \(N_1\) and \(N_2\) be the sample sizes for traits 1 and 2, and let \(N_s\) be the number of participants present in both samples. The overlap-induced sampling covariance is the contribution to the covariance of the two association statistics that arises because the same participants contribute to both. Under a simple standardized-trait model, it is approximately \(N_s\rho_o/\sqrt{N_1N_2}\), where \(\rho_o\) is the correlation between the traits among overlapping participants.

The approximation assumes that the traits and genotypes have been handled on compatible scales, that association statistics are appropriately standardized, and that the relevant overlap participants have a common trait correlation. Real analyses may require adjustments for covariates, case-control sampling, or other design features. The formula is a useful baseline, not a substitute for checking how the studies and statistics were constructed.

Deriving the Overlap Contribution

Consider one variant and, for clarity, suppose each trait is measured on a standardized scale. Let \(g_i\) denote the centered, variance-one genotype for participant \(i\), and let \(u_{1i}\) and \(u_{2i}\) denote the standardized residual values for the two traits under a null model for this variant. In a simple score-test approximation, the association statistics have the form

$$ z_1 \approx \frac{1}{\sqrt{N_1}}\sum_{i\in S_1}g_i u_{1i}, \qquad z_2 \approx \frac{1}{\sqrt{N_2}}\sum_{i\in S_2}g_i u_{2i}. $$

Here \(S_1\) and \(S_2\) are the two study samples. Suppose participants are independent of one another, genotype is independent of the residual traits under the tested-variant null, and the residual correlation for an overlapping participant is \(\rho_o\). Terms from different participants have zero covariance under these assumptions. The only shared terms are from the \(N_s\) overlapping participants. For each such participant, the covariance of the two score contributions is

$$ \operatorname{Cov}(g_i u_{1i},g_i u_{2i}) = E(g_i^2 u_{1i}u_{2i}) = E(g_i^2)E(u_{1i}u_{2i}) = \rho_o, $$

because \(E(g_i^2)=1\) and the residual correlation is \(\rho_o\). Summing the \(N_s\) shared contributions and applying the two standardizing factors gives

$$ \operatorname{Cov}(z_1,z_2) \approx \frac{N_s\rho_o}{\sqrt{N_1N_2}}. $$

The sign follows the trait correlation among overlapping participants. Positive correlation produces positive overlap covariance; negative correlation produces negative covariance. The term is dimensionless, like the product of two \(z\)-statistics, and it is not itself a genetic correlation.

Worked Example: Calculating the Overlap Contribution

Suppose two hypothetical studies have \(N_1=80{,}000\) and \(N_2=50{,}000\) participants, with \(N_s=20{,}000\) participants in common. Assume the traits’ correlation among overlapping participants is \(\rho_o=0.30\). Under the simple standardized-trait approximation, calculate the overlap contribution to the covariance of their association statistics.

Use the overlap formula:

$$ \frac{N_s\rho_o}{\sqrt{N_1N_2}} = \frac{20{,}000(0.30)}{\sqrt{80{,}000(50{,}000)}}. $$

The numerator is \(6{,}000\), and the denominator is \(\sqrt{4{,}000{,}000{,}000}\approx 63{,}246\). Therefore,

$$ \frac{6{,}000}{63{,}246}\approx 0.09487. $$

The overlap contributes approximately \(0.0949\) to the covariance of the two \(z\)-statistics under this approximation. This is an expected covariance contribution, not a claim that every pair of observed \(z\)-statistics will have product \(0.0949\), and it does not measure the traits’ genetic correlation.

How Overlap Appears in Cross-Trait LDSC

For variant \(j\), let \(z_{1j}\) and \(z_{2j}\) be its association statistics in the two studies, and let \(\ell_j\) be its LD score, defined as the sum of squared correlations between that variant and the variants included in the analysis. Under the standard cross-trait LDSC approximation, the expected product of the association statistics can be written as

$$ E[z_{1j}z_{2j}] \approx \frac{\sqrt{N_1N_2}}{M}C_g\,\ell_j + \frac{N_s}{\sqrt{N_1N_2}}\rho_o. $$

Here \(M\) is the number of variants in the genetic-architecture model, and \(C_g\) is the total additive genetic covariance of the two standardized traits on the scale used by the model. The first term is the genetic component: its expected size changes with LD score. The second is the overlap contribution: in the simple model, it does not depend on the tested variant’s LD score. In practice, the fitted intercept can also absorb other approximately LD-independent contributions, such as certain forms of confounding.

This equation explains why overlap can matter even if the scientific target is genetic covariance. If it is ignored, an upward or downward shift in the association-statistic products can be mistaken for a genetic pattern, particularly when the analysis or its assumptions do not adequately separate the intercept from the slope. LDSC uses variation in LD score to estimate the LD-dependent slope and the intercept; it does not identify the cause of every intercept contribution.

Worked Example: Negative Trait Correlation Among Overlapping Participants

Consider hypothetical studies with \(N_1=100{,}000\), \(N_2=64{,}000\), and \(N_s=16{,}000\). Suppose the traits have correlation \(\rho_o=-0.25\) among participants in both studies. Calculate the overlap contribution and interpret its sign.

Substitute the values into the formula:

$$ \frac{N_s\rho_o}{\sqrt{N_1N_2}} = \frac{16{,}000(-0.25)}{\sqrt{100{,}000(64{,}000)}} = \frac{-4{,}000}{80{,}000} = -0.05000. $$

The expected overlap contribution is \(-0.0500\). Under the stated assumptions, overlap pushes the association-statistic covariance downward because the two traits are negatively correlated in the shared participants. This does not imply negative genetic covariance: the genetic component is a separate term and could be positive, negative, or near zero.

Using the Intercept as a Diagnostic

If the number of overlapping participants and their trait correlation are known, the formula provides a predicted overlap contribution to compare with the fitted cross-trait LDSC intercept. Agreement is consistent with overlap explaining much of the intercept under the assumed model. A discrepancy is a reason to investigate the assumptions and other sources of intercept inflation; it does not by itself prove that the overlap count or the LDSC analysis is wrong.

Conversely, an intercept can sometimes be used to make a rough overlap calculation if the overlap trait correlation is known and other intercept sources are negligible. Rearranging the formula gives

$$ N_s \approx \frac{I\sqrt{N_1N_2}}{\rho_o}, $$

where \(I\) is the intercept contribution being attributed to overlap. This inversion is unstable when \(\rho_o\) is close to zero, and it is not justified when the intercept also contains meaningful confounding or other contributions. An intercept is therefore a diagnostic, not a definitive census of shared participants.

Worked Example: Comparing Predicted Overlap With an LDSC Pattern

Suppose hypothetical studies have \(N_1=90{,}000\), \(N_2=40{,}000\), and \(N_s=18{,}000\), and the correlation among overlapping participants is estimated as \(\rho_o=0.25\). Calculate the predicted overlap contribution. Then suppose the mean products of association statistics at two LD-score values are \(0.095\) at \(\ell=5\) and \(0.135\) at \(\ell=15\). Use these values to illustrate the intercept-versus-slope distinction.

First calculate the predicted overlap contribution:

$$ \frac{18{,}000(0.25)}{\sqrt{90{,}000(40{,}000)}} = \frac{4{,}500}{60{,}000} = 0.07500. $$

The two illustrative mean products imply a change of \(0.135-0.095=0.040\) across an LD-score change of \(15-5=10\). The corresponding slope is

$$ \frac{0.040}{10}=0.00400 \quad\text{per unit of LD score}. $$

Extrapolating this straight-line pattern to LD score zero gives an intercept of \(0.095-(0.00400)(5)=0.0750\), matching the overlap contribution predicted from the supplied sample sizes and trait correlation. In an actual analysis, LDSC fits a regression using the variant-level data and its weighting and modeling assumptions; two illustrative points do not establish that overlap is the only intercept source. The increasing products with LD score are consistent with a positive LD-dependent component, which is the part relevant to genetic covariance under the model.

Assumptions, Interpretation, and Common Pitfalls

The simple calculation is most direct when the association statistics are comparably standardized, the sample sizes are correctly represented, and \(\rho_o\) describes the relevant trait correlation among the overlapping participants. If the shared participants differ from the full studies, a single correlation may be an imperfect summary. Covariate adjustment, relatedness, ascertainment, or different phenotype scales can also make the simple approximation less exact. In applied work, use overlap counts and trait-correlation estimates that match the analysis design, and interpret the LDSC intercept in light of other plausible sources of covariance.

A key distinction is between the correlation of observed trait values and genetic covariance. Correlated observed traits among shared participants can induce correlated sampling errors. Genetic covariance, by contrast, concerns shared additive genetic effects and is represented by the LD-dependent component in the cross-trait LDSC model. The two components can have different magnitudes or even opposite signs.

  • Mistake: “A positive intercept proves shared genetic effects.” The intercept can reflect sample overlap and other approximately LD-independent effects. The model associates shared polygenic effects with the LD-dependent slope, not with an intercept alone.
  • Mistake: “No overlap means no covariance between the statistics.” No overlap removes this particular shared-participant contribution under the simple model. It does not remove covariance caused by shared genetic effects or other sources that affect both studies.
  • Mistake: “The overlap formula gives the genetic correlation.” It gives an approximate sampling-covariance contribution to the association statistics. Genetic correlation is a distinct parameter based on genetic variances and genetic covariance.
  • Mistake: “The fitted intercept exactly equals the overlap term.” That equality requires the overlap model to be adequate and other intercept contributions to be negligible. Treat the comparison as a diagnostic rather than an identity guaranteed by the data.
Key takeaway: Under a simple standardized-trait model, overlapping participants contribute approximately \(N_s\rho_o/\sqrt{N_1N_2}\) to the covariance of two association statistics. In cross-trait LDSC, this is an approximately LD-independent contribution, while shared polygenic effects contribute an LD-dependent pattern. Interpret the intercept with overlap information and other possible confounding in view; neither overlap nor the intercept alone establishes the genetic relationship between traits.

Check Your Understanding

Use the definitions and assumptions in this tutorial to answer the following questions.

  1. Two studies have sample sizes \(N_1\) and \(N_2\), with \(N_s\) participants in common. What additional quantity is needed to calculate the simple overlap contribution to association-statistic covariance?
  2. If the traits are negatively correlated among overlapping participants, what is the expected sign of the overlap contribution? Does that sign determine the sign of genetic covariance?
  3. In the cross-trait LDSC expectation, which component is associated with LD-score variation, and which is approximately constant across LD score under the simple model?
  4. Why is a fitted positive LDSC intercept not, by itself, proof of either sample overlap or shared genetic effects?
  5. When might calculating an overlap count by rearranging the intercept formula be unreliable?
Clinical Biostats