Tutorials › Statistical Genetics › Reference Databases for Rare-Variant Interpretation

Rare-Variant Data and Inference · Tutorial 401 of 989

Reference Databases for Rare-Variant Interpretation

Evaluate what a reference database can—and cannot—tell you about a variant’s frequency in the population relevant to an interpretation.

Intermediate 10 min read

What You'll Learn

  • Distinguish a recorded zero allele count from a variant that was not assessed or reported.
  • Compare allele counts using the number of callable chromosomes as the denominator.
  • Explain how ancestry representation affects the relevance and precision of a frequency estimate.
  • Use a zero-count bound to judge what an absence can support.
  • Recognize how sequencing, quality control, and variant representation affect database evidence.
  • Combine database evidence cautiously rather than treating one resource as a definitive frequency answer.

Why a Reference Database Needs Context

A reference database can help assess how often a genetic variant has been observed in sampled people. But its value is not determined by the database’s total number of participants alone. The relevant questions include how many chromosomes were actually assessed at the site, whether the sampled people represent the population of interest, and whether the data are sufficiently reliable for the type of variant being considered.

In particular, “not listed” can have several meanings. A variant may have been assessed and observed zero times; the site may not have been adequately assessed; or the database may not report variants that fail its inclusion or quality rules. Those situations provide different evidence. An absent record is not, by itself, proof that a variant is exceptionally rare.

Definition: A reference database is a collection of genetic data from sampled individuals that can be used to describe observed variation and estimate allele frequencies in the represented samples. For a variant, the allele count (AC) is the number of observed copies of that allele, and the allele number (AN) is the number of chromosome copies successfully assessed at the site. When the denominator is appropriate, the observed allele frequency is \( \hat{p}=\mathrm{AC}/\mathrm{AN} \). A database frequency describes the sampled data; how well it represents a target population depends on sampling, ancestry representation, and data quality.

For a diploid autosomal site, a fully assessed individual can contribute two chromosome copies, but missing or unusable data can reduce the number assessed. Thus, participant count and AN are not interchangeable. Frequency estimates should be considered within the relevant ancestry group when possible, not only in a pooled sample whose composition may differ from the population of interest.

Read the Denominator Before the Frequency

The observed frequency is meaningful only alongside its denominator. A count of one allele copy among 100 assessed chromosomes gives a different estimate, and much more sampling uncertainty, than a count of one among 100,000. The key quantity is the number of relevant chromosome copies successfully assessed, not just the headline number of database participants.

Population-specific frequency matters for the same reason. If an allele occurs in one ancestry group but is uncommon in a much larger group, pooling those groups can make the pooled frequency a poor summary for either group. Ancestry categories in databases are approximations, and they may not capture all relevant population structure. A frequency estimate should therefore be matched as closely as possible to the population about which an interpretation is being made.

Worked Example: Interpreting a Small Allele Count

Suppose a hypothetical reference database reports a variant with AC = 3 and AN = 12,000 in an ancestry group of interest. Calculate the observed frequency and describe what the result means.

Calculation. The observed frequency is the allele count divided by the assessed allele number:

$$ \hat{p}=\frac{\mathrm{AC}}{\mathrm{AN}}=\frac{3}{12{,}000}=0.00025 $$

The variant was observed in 3 of the 12,000 chromosome copies assessed, for an observed frequency of 0.00025, or 0.025% (about 1 in 4,000 chromosome copies in these sampled data). This does not mean that exactly 1 in 4,000 people in every population carries the variant. The estimate is based on a small count, refers to chromosome copies rather than people, and applies most directly to the sampled ancestry group and database methods.

What a Zero Count Can Establish

A reported zero allele count can still be informative if the database provides a credible denominator for the relevant site and population. Under a simple sampling model, suppose \(N\) independent chromosome copies are assessed and each has probability \(p\) of carrying the allele. The number of observed copies follows a binomial model. The probability of seeing zero copies is \( (1-p)^N \). A zero count is more surprising when \(N\) is large than when \(N\) is small, if the true frequency is held fixed.

A useful way to express the limit of a zero count is a one-sided 95% upper bound. It is the frequency \(p_U\) for which the probability of observing zero copies is 0.05:

$$ (1-p_U)^N=0.05 \qquad\Longrightarrow\qquad p_U=1-0.05^{1/N} $$

For large \(N\), this is approximately \(3/N\). The bound is not an estimate that the frequency equals \(3/N\); it indicates that, under the model, frequencies substantially above that value would make a zero count unlikely. This calculation assumes independent sampled chromosome copies and a common frequency \(p\) among them. It does not fix a mismatch between the database sample and the target population, or compensate for a site that was not adequately assessed.

Worked Example: Comparing Zero Counts at Different Sample Sizes

Two hypothetical ancestry-specific database summaries report no observed copies of a variant. The first has AN = 20,000; the second has AN = 400. Use the approximate one-sided 95% upper bound \(3/N\) to compare the evidence.

For AN = 20,000:

$$ \frac{3}{20{,}000}=0.00015 $$

The approximate upper bound is 0.00015, or 0.015%. Under the stated sampling model, a zero count in this many assessed chromosome copies is less compatible with frequencies much above that level.

For AN = 400:

$$ \frac{3}{400}=0.0075 $$

The approximate upper bound is 0.0075, or 0.75%. A zero count in this smaller sample is compatible with a substantially higher frequency than the first zero count. Therefore, the two absences do not provide equally strong evidence about rarity. Neither bound should be applied to a different population without justification, and both rely on the assumed sampling model.

Ancestry Representation and Pooled Summaries

A large total database can still offer limited evidence for a particular ancestry group if relatively few people from that group were included or if relatively few chromosome copies were assessable. A pooled frequency can conceal this imbalance. For interpretation, examine both the overall summary and the ancestry-stratified counts and denominators, while recognizing that ancestry labels are imperfect proxies for genetic and population history.

Worked Example: Why a Pooled Frequency Can Mislead

Consider a hypothetical database with two ancestry groups. In group A, AC = 2 and AN = 8,000. In group B, AC = 0 and AN = 500. Calculate the observed frequencies and the pooled frequency, then assess what the zero in group B supports.

Group A:

$$ \hat{p}_A=\frac{2}{8{,}000}=0.00025 $$

Group B:

$$ \hat{p}_B=\frac{0}{500}=0 $$

The pooled count is \(2+0=2\), and the pooled allele number is \(8{,}000+500=8{,}500\). Thus:

$$ \hat{p}_{\mathrm{pooled}}=\frac{2}{8{,}500}\approx0.0002353 $$

The pooled frequency is about 0.000235, but it should not replace the group-specific information. In group B, a zero count among 500 assessed chromosome copies gives an approximate 95% upper bound of \(3/500=0.006\), or 0.6%. That is far above group A’s observed frequency. The data therefore do not establish that the allele is absent, or equally rare, in group B; the smaller denominator leaves considerable uncertainty. The pooled value mostly reflects the much larger group A sample.

Assess Data Quality and Database Scope

A database frequency is useful only when the allele and its denominator are credible. Before treating a reported count or absence as evidence, consider whether the database’s scope and data-generation procedures fit the question. Quality is not a single yes-or-no property: performance can differ by genomic region, variant type, sequencing methods, and filtering rules.

  • Site assessment: Determine whether the site was assessed in the relevant samples. An AN of zero or a missing denominator is not evidence for a zero frequency. A reported zero count is informative only when an adequate number of copies was actually assessed.
  • Variant type and region: Some variant classes or difficult genomic regions may be less reliably represented by a resource’s methods. A database designed primarily for one data type may not provide equally strong evidence for every other type.
  • Quality and inclusion rules: Resources may exclude low-confidence observations or omit variants that do not meet reporting thresholds. A missing entry can therefore reflect filtering or reporting policy, not necessarily a verified absence in all sampled individuals.
  • Consistent variant representation: Confirm that the database record refers to the same genomic position and alternate allele as the variant being interpreted. Differences in genome build, allele alignment, or representation of complex variants can make matching records difficult.
  • Sample composition: Consider how participants were selected and whether their ancestry and other characteristics are relevant. A database sample is not automatically a random, representative sample of every population.

These considerations do not require treating every database result as unusable. Rather, they determine how much weight the result can bear. A well-supported count in a relevant group can inform frequency assessment; a missing record without a reliable site-specific denominator usually supports a much weaker conclusion.

Worked Example: Distinguishing Zero From Not Assessed

Suppose a hypothetical database summary shows two entries for a variant site in a target ancestry group. Entry A reports AC = 0 and AN = 620. Entry B has no variant record and supplies no site-specific AN. Use the approximate zero-count bound to describe what can be concluded from Entry A, and compare it with Entry B.

For Entry A, the reported denominator shows that 620 chromosome copies were assessed. Its approximate one-sided 95% upper bound is:

$$ \frac{3}{620}\approx0.004839 $$

This corresponds to about 0.484%. Under the zero-count sampling assumptions, the observation is evidence against frequencies much higher than approximately 0.00484 in the represented sample, but it does not establish that the frequency is exactly zero.

For Entry B, there is no reported variant and no site-specific allele number. Without evidence that the site was assessed and reportable in the relevant samples, there is no valid denominator for a zero-count bound. The missing record could reflect lack of assessment, filtering, or reporting conventions. It therefore cannot be interpreted like Entry A’s documented zero count.

Putting Database Evidence Together

A careful interpretation treats database information as one component of evidence, not as a verdict. The practical sequence is to identify the exact variant and relevant population, check the count and denominator, judge the precision of any zero count, and evaluate whether data quality and database scope support the comparison. If sources disagree, differences in ancestry composition, sample size, site assessment, or representation may explain some of the discrepancy; the largest database total does not automatically settle the question.

1
Match the variant.
Confirm that position, allele, and representation correspond to the variant under review.
2
Match the population.
Inspect ancestry-specific counts where available and consider whether the sampled group is relevant to the interpretation.
3
Inspect AC and AN.
Calculate AC/AN for an observed frequency. For a zero count, first verify that an adequate site-specific AN is reported.
4
Judge scope and quality.
Consider whether the database’s methods, variant coverage, and reporting rules make the result reliable for this variant.
5
State the conclusion narrowly.
Describe what was observed in the database sample and avoid claiming that an unobserved allele is absent from the population.

Common Mistake: Treating “Not Listed” as “Exceptionally Rare”

The incorrect reasoning is: “The database has no record of this variant, so the variant must be exceptionally rare.” This conclusion skips two questions: was the relevant site actually assessed, and how many relevant chromosome copies were assessed? Even a genuine zero count can be weak evidence when the denominator is small or the target ancestry is underrepresented. And when there is no site-specific denominator, a missing entry is not a measured zero.

The correction is to separate absence of a record from zero observed copies among assessed chromosomes. For the latter, report the ancestry group, AC, AN, and the resulting frequency or uncertainty. For the former, do not infer a frequency unless the database documentation supports that inference. In either case, describe the limits of the database evidence and avoid converting a sample observation into a claim of population-wide absence.

Key takeaway: A reference database supports a frequency claim only to the extent that the relevant allele was consistently represented, adequately assessed, and observed in a sample relevant to the population of interest. “Not listed” alone does not establish exceptional rarity.

Check Your Understanding

Use the database counts, denominators, and sampling assumptions carefully.

  1. A database reports AC = 4 and AN = 16,000 in an ancestry group. What is the observed allele frequency, and what does it describe?
  2. Why does a zero count among 300 assessed chromosome copies provide weaker evidence about rarity than a zero count among 30,000?
  3. A variant is absent from a database’s variant list, but the site-specific allele number is unavailable. What additional information is needed before interpreting that absence?
  4. A pooled database frequency is low, but the ancestry group of interest has a much smaller denominator. Why should the pooled result not be treated as the group-specific frequency?
  5. Name two data-quality or database-scope issues that could make an absent record uninformative about the population frequency.
Clinical Biostats