Tutorials › Statistical Genetics › Design Choices That Undermine GWAS Reproducibility

Genome-Wide Association Study Design · Tutorial 357 of 989

Design Choices That Undermine GWAS Reproducibility

Learn to identify design mismatches that make GWAS results difficult to reproduce and to distinguish biological differences from differences in measurement, recruitment, precision and variant coverage.

Intermediate 10 min read

What You'll Learn

  • Define reproducibility in terms of comparable phenotypes, variants, effect estimates and study populations
  • Diagnose how phenotype-definition drift changes the association being estimated
  • Explain how selective recruitment can limit generalizability or introduce selection bias
  • Assess how sample size affects the precision of an association estimate
  • Identify ancestry and platform differences that complicate comparisons across studies
  • Use a design audit to distinguish a genuine inconsistency from a mismatch between studies

What Reproducibility Requires in a GWAS

A genome-wide association study (GWAS) result is not just a variant and a p-value. It is an estimate tied to a particular phenotype definition, recruited population, genetic variant or proxy, measurement platform, and analysis. If a later study changes one or more of these, its estimate may differ even when both studies have been conducted carefully. Conversely, repeating the same design does not guarantee identical estimates: estimates vary because each sample is finite.

In this tutorial, reproducibility means that a result can be evaluated in a new study with sufficiently comparable measurements and design, and that the estimates are reasonably consistent in direction and magnitude given their uncertainty. It does not mean that two studies must produce identical estimates or identical p-values. The first diagnostic question is therefore not simply “Did the result replicate?” but “Were the studies estimating and measuring sufficiently similar things?”

Definition: A GWAS finding is reproducible to the extent that an independent study using comparable phenotype definitions, sampling and analysis, and comparable genetic measurements produces evidence consistent with the original association. A meaningful comparison requires checking what each study measured and in whom, as well as which variant or genetic signal was analyzed.

Several design choices can make that comparison difficult. Phenotype drift changes the outcome definition. Selective recruitment can change who enters the study and, in some circumstances, create or distort associations. A small sample yields imprecise estimates. Ancestry underrepresentation limits the populations in which the result can be assessed and may change the performance of a marker. Platform limitations can leave a variant unmeasured or poorly represented. These issues can occur together.

Phenotype Drift: Different Definitions, Different Outcomes

A phenotype is an operational definition of the trait or disease used in an analysis. It may come from clinical records, laboratory measurements, medication records, questionnaires, or a combination. Two studies can use the same disease name while applying different rules for who counts as a case. They can also use nominally identical rules but differ in measurement timing, instruments, thresholds, or data availability.

These differences matter because an association estimate refers to the outcome as defined in that study. For example, an association with a carefully adjudicated diagnosis is not automatically the same estimand as an association with self-reported history of that diagnosis. A broader definition may include milder cases or people whose reports have not been clinically confirmed. A stricter definition may exclude some true cases while yielding a more specific case group. Neither choice is universally preferable; the key is to state the definition and assess its comparability.

Worked Example: Auditing a Changed Case Definition

Suppose a hypothetical discovery study defines a condition using a clinical-record algorithm requiring two diagnostic entries. A replication study instead asks participants whether a clinician has ever told them they have the condition. In a separate hypothetical review of 500 records, 100 people meet the clinical algorithm, 70 meet the self-report definition, and 50 meet both definitions.

The clinical-only group contains \(100-50=50\) people. The self-report-only group contains \(70-50=20\). Thus, \(50+20=70\) people meet exactly one of the two definitions, while the union of the definitions contains \(100+70-50=120\) people. The fraction of people meeting at least one definition who are classified differently is \(70/120=0.583\), or 58.3%.

This calculation does not show that either definition is wrong, nor does it predict the direction or size of a genetic association. It shows that the labels “case” and “control” do not select the same people under these rules. Before interpreting different GWAS estimates as a failure to reproduce, investigators should compare the case algorithms, timing, and any validation data, and consider analyzing harmonized definitions where possible.

Phenotype harmonization makes definitions sufficiently comparable for a planned analysis; it does not make every source of measurement error disappear. If studies cannot use identical data, they should document the differences, describe their likely consequences, and avoid presenting the results as though the outcomes were interchangeable.

Recruitment and Selection: Who Enters the Study?

Recruitment determines the population represented by the analyzed sample. A hospital-based study, a volunteer biobank, and a population-based cohort may enroll people with different health profiles, access to care, socioeconomic circumstances, and willingness to participate. A result in a recruited sample does not automatically describe every population from which participants might have come.

Selection is especially concerning when inclusion or participation depends on factors related to both genotype and the phenotype, or to causes of the phenotype. Conditioning on such selection can induce or distort an association. This is one form of selection bias. Recruitment alone does not prove that bias occurred: a case-control design, for example, can be appropriate when its sampling scheme and analysis match the intended comparison. The risk depends on how inclusion relates to the variables involved, not simply on whether the sample is clinic-based or volunteer-based.

Worked Example: Comparing Clinic and Population Recruitment

Consider two hypothetical studies of the same disease. One recruits cases from a specialist clinic and controls from patients attending that clinic for unrelated conditions. The other recruits cases and controls from a defined population. In a recruitment audit, 40% of the clinic controls and 10% of the population controls have a particular chronic condition that can affect clinic attendance.

The difference is \(0.40-0.10=0.30\), or 30 percentage points. This is evidence that the control groups differ in a measured feature relevant to recruitment; it is not, by itself, proof that the genetic association is biased. To assess the risk, investigators would ask whether the chronic condition, or factors associated with it, relate to the genotype or disease outcome, and whether selection into the clinic depends on those factors.

For a replication comparison, the studies should report their recruitment sources and eligibility rules, compare relevant characteristics, and state the population each estimate represents. If recruitment processes differ substantially, a difference in estimates may reflect different sampled populations or selection mechanisms rather than a contradiction about the same target population.

A careful interpretation separates two questions: whether an association is present in the analyzed sample, and whether it generalizes to another target population. The second question needs evidence about the recruitment process and population differences; it does not follow automatically from the first.

Sample Size: Inconsistent Estimates Can Be Noise

A small study may estimate an effect imprecisely even when its design and phenotype match a larger study. Under a simple additive linear model for a quantitative trait, let \(G\) be genotype dosage, coded as 0, 1, or 2 copies of an allele, and let \(\beta\) be the mean trait difference per additional allele. If dosage variance is \(\operatorname{Var}(G)\) and the residual standard deviation is \(\sigma\), an approximate standard error under a simple model is

$$ SE(\hat{\beta}) \approx \frac{\sigma}{\sqrt{n\operatorname{Var}(G)}}. $$

This approximation illustrates why larger samples usually yield more precise estimates when the phenotype, allele-frequency distribution, and residual variability are otherwise comparable. Under Hardy–Weinberg equilibrium for a biallelic variant with allele frequency \(f\), genotype variance is \(2f(1-f)\). Real analyses may include covariates and other complications, so this formula is a diagnostic illustration rather than a substitute for a design-specific calculation.

Worked Example: How Sample Size Changes Precision

Suppose two hypothetical studies analyze the same quantitative trait and variant, with \(f=0.20\), residual standard deviation \(\sigma=1\) trait unit, and comparable designs. The first has \(n=10{,}000\); the second has \(n=2{,}500\). Under the stated approximation, \(\operatorname{Var}(G)=2(0.20)(0.80)=0.32\).

For the first study, \(SE(\hat{\beta})\approx 1/\sqrt{10{,}000(0.32)}=1/\sqrt{3{,}200}=0.01768\) trait units per allele. For the second, \(SE(\hat{\beta})\approx 1/\sqrt{2{,}500(0.32)}=1/\sqrt{800}=0.03536\) trait units per allele. The second standard error is twice as large.

If the underlying per-allele effect were \(0.05\) trait units, the effect divided by the approximate standard error would be \(0.05/0.01768=2.83\) in the first study and \(0.05/0.03536=1.41\) in the second. These ratios illustrate the difference in precision; they do not guarantee that either study will observe that exact ratio or establish a result by themselves. In particular, a less precise estimate may differ substantially from the underlying effect through sampling variation.

A reproducibility assessment should compare estimates and their uncertainty, not just whether each study crosses a decision threshold. A small study can yield an estimate in the same direction but with wide uncertainty; it can also yield an estimate in the opposite direction by chance. Conversely, increasing sample size cannot repair a different phenotype definition, systematic selection bias, or a variant that was not measured comparably.

Ancestry Representation and Variant Comparability

Recruitment plans that substantially underrepresent some ancestry groups restrict the populations in which a study can directly evaluate its findings. “Ancestry” in genetic analyses refers to patterns of genetic similarity and shared population history; it should not be treated as interchangeable with a broad social or racial label. Both genetic ancestry and the social conditions that shape recruitment and health can matter, but they are distinct considerations.

A variant associated with a trait may itself be causal, or it may be a marker correlated with a causal variant through linkage disequilibrium. Allele frequencies and linkage disequilibrium patterns can differ among populations. Consequently, a marker can have different precision or represent a causal signal less well in a population not adequately represented in discovery. A difference in the marker’s association estimate does not automatically imply that the underlying biology differs.

Worked Example: Allele Frequency and Precision Across Samples

Suppose, hypothetically, that the same additive effect is present in two equally sized samples of \(n=10{,}000\), with the same residual standard deviation of 1 trait unit. The allele frequency is \(f=0.30\) in one sample and \(f=0.05\) in the other. Using \(\operatorname{Var}(G)=2f(1-f)\), the respective genotype variances are \(2(0.30)(0.70)=0.42\) and \(2(0.05)(0.95)=0.095\).

The approximate standard errors are \(1/\sqrt{10{,}000(0.42)}=0.01543\) and \(1/\sqrt{10{,}000(0.095)}=0.03245\) trait units per allele. Their ratio is \(0.03245/0.01543=2.10\), so the lower-frequency variant has about 2.10 times the standard error under these simplified assumptions.

This calculation concerns precision for the measured variant. It does not account for differences in linkage disequilibrium with a causal variant, which can further change the association of a marker. Therefore, a study that adds participants from an underrepresented ancestry group should not assume that a marker will have the same frequency, tagging performance, or precision there. It should assess variant comparability and report which populations its results directly support.

Platform Limitations and an Integrated Design Audit

Genotyping arrays, exome sequencing, whole-genome sequencing, and imputation do not provide identical coverage of genetic variation. A variant may be directly measured in one study, inferred with varying quality in another, or absent from the available data. A nearby marker can serve as a proxy only to the extent that it tracks the relevant signal in that population. Thus, “the same locus” does not always mean that studies tested the same variant or an equally informative proxy.

For example, if a discovery study tests a directly genotyped variant but a replication study lacks that variant and uses a proxy, investigators should document the proxy, its imputation or genotyping quality, and its relationship to the discovery variant in the study population. Low-quality measurement can reduce information and make estimates less stable; systematic measurement problems can also affect validity. The precise consequence depends on the platform and data processing, so it should not be inferred from the word “imputed” alone.

A practical audit proceeds by checking comparability in a deliberate order. This does not replace statistical analysis; it helps identify whether two studies are addressing the same question with comparable information.

1
Compare the phenotype.
Check case and control rules, measurement source, timing, thresholds, and missing-data handling. Record any definition changes that could alter who is classified or the trait scale.
2
Compare recruitment.
Identify the source population, eligibility criteria, and participation process. Ask whether the two samples represent comparable target populations and whether selection could relate to genotype and outcome or their causes.
3
Compare precision and ancestry representation.
Review sample size and estimate uncertainty, then check which ancestry groups are represented. Consider allele-frequency differences and whether the tested marker tags the same signal.
4
Compare genetic measurements.
Verify that the same variant, allele coding, and sufficiently comparable genotyping or imputation information were used. A different marker or allele orientation can make an apparent mismatch misleading.

This audit matters because “not reproduced” can describe several different situations: a genuinely inconsistent estimate, an imprecise estimate, a changed phenotype, a selected sample, a different marker, or a population in which the discovery marker is less informative. These explanations are not interchangeable, and none should be invoked automatically to excuse a disagreement. The study reports and data should determine which are plausible.

Common Mistake

Common mistake: “The replication study did not produce the same p-value, so the original association was false.” This reasoning treats p-values as if they were stable properties of a variant. They depend on the estimate, its precision, and the analysis. The correction is to compare effect estimates and uncertainty, then check phenotype, recruitment, ancestry representation, and variant measurement. A different p-value alone does not establish either replication or non-replication.

The reverse mistake is equally important: attributing every discrepancy to design differences without evidence. Reproducibility is strengthened by transparent definitions, representative recruitment, adequate precision for the intended comparison, and comparable genetic measurements. If a study differs from its predecessor, report the difference and assess its consequences rather than assuming that it explains the result.

Key takeaway: Before interpreting a GWAS result as reproduced or not reproduced, establish whether both studies measured comparable phenotypes, sampled comparable populations, had enough precision for the intended comparison, and analyzed comparable genetic signals. Design differences can explain difficulty in comparison, but they must be evaluated rather than presumed.

Check Your Understanding

Use the design-audit principles to answer the following questions.

  1. Why can two studies using the same disease name still be analyzing different phenotypes?
  2. What recruitment information would help you assess whether a clinic-based study and a population-based study are comparable?
  3. In the approximate quantitative-trait model, what happens to the standard error when sample size is reduced by a factor of four, with other quantities unchanged?
  4. Why might a marker’s association differ across ancestry groups even if the effect of a causal variant were similar?
  5. What information should be checked when a replication study uses an imputed proxy rather than the discovery variant?
Clinical Biostats