fe0627da07b9f989bc642d824af381968163d9a8
lrnassar
Wed Jul 8 13:40:04 2026 -0700
varFreqs: move ALFA to Global reference panels; add shared-sample bias caveat to description pages. refs #36642
Per Max's post-release feedback on the ticket: ALFA is not a European
cohort, it is an NCBI aggregator of dbGaP studies from many sources.
Moved from Europe to Global reference panels in the newsarch entry.
Also adds a shared-sample bias caveat to the Pooled allele frequency
section of both varFreqsAffected.html and varFreqsBackground.html. The
background page cites the concrete overlaps (1000 Genomes in both gnomAD
HGDP+1kG and HRC; HGDP/SGDP; AllOfUs/TOPMed; ALFA aggregating dbGaP
studies used elsewhere). The affected page cites the SPARK WGS cohort
being a subset of SPARK WES. Both explain that pooled AN is inflated and
pooled AF is skewed toward the frequency in the shared subset, and point
users to the per-cohort AC/AF/AN fields for unbiased single-cohort
numbers.
diff --git src/hg/makeDb/trackDb/human/varFreqsBackground.html src/hg/makeDb/trackDb/human/varFreqsBackground.html
index 598b693e6ba..56fdfc599b7 100644
--- src/hg/makeDb/trackDb/human/varFreqsBackground.html
+++ src/hg/makeDb/trackDb/human/varFreqsBackground.html
@@ -48,30 +48,40 @@
unaffected/control arms: backgroundAF = sum(AC) / sum(AN), where
backgroundAC sums the allele counts and backgroundAN sums the allele
numbers across each cohort/arm that provides both AC and AF (the per-arm AN is derived as
round(AC / AF)). Two cohorts that publish only AF (ABraOM, ALFA) are still
pooled by assigning them an assumed allele number, set as a default_an in the
build configuration; their per-arm AC is then derived as round(AF × default_an).
Cohorts that publish
only AC with no default_an set (currently MGRB and the GREGoR unaffected and
unknown arms), and cohorts that contribute only through per-population AC/AF (currently
AllOfUs), are listed in backgroundSources but do not contribute to the pool
numerator or denominator; their data remain visible in the per-database and per-population
AC/AF columns. The pooled rate is preferred over a max-across-cohorts statistic so a small
cohort with a high local AF (for example AllOfUs Oceanian) cannot dominate the displayed
frequency.
+The pooled rate also inherits a shared-sample bias: several source cohorts overlap +in the individuals they include. For example, 1000 Genomes samples appear in both gnomAD +HGDP+1kG and HRC; HGDP and SGDP overlap; AllOfUs and TOPMed share participants; ALFA +aggregates dbGaP studies used elsewhere in the pool. Where a variant sits in a shared +sample, both its AC and AN are counted more than once, so pooled AN is inflated and +pooled AF is skewed toward the frequency in the shared subset. Treat the pooled rate as +a cross-cohort summary rather than an unbiased population estimate; the per-cohort +AC/AF/AN fields on each variant give the single-cohort numbers. +
Alongside the pooled rate, the mouseover lists the top 3 contributing
background sources ranked by their own per-source AF, formatted as
Source (AF). This surfaces population cohorts where a variant
is specifically enriched, even when the pooled rate is small; the
East-Asian founder allele
rs4986893,
for example, ranks ToMMo Japan and KOVA Korea at the top while the pooled
rate across all contributing sources sits much lower. For disease cohorts
that ship a phenotype split (SPARK, SFARI WGS, SCHEMA, GREGoR), the
displayed AF is the unaffected-arm AF and the label includes the arm (for
example SPARK non-ASD, SCHEMA ctrl); for
population cohorts, the label is the cohort name and the AF is the unified