This is a rigorously scoped, purely comparative re-analysis: the authors extracted variants for eight globin genes (HBA1, HBA2, HBZ, HBB, HBD, HBG1, HBG2, HBE1) from gnomAD v2.1.1 and v4.1 callsets, harmonized them by locus and alleles, applied ClinGen Hemoglobinopathy VCEP thresholds for BA1, BS1, BS2, BS2_Supporting and PM2_Supporting, and cross-referenced ClinVar (accessed 10 Oct 2025). Sample scale is large: v4.1 includes 730,947 exomes and 76,215 genomes vs 125,748 exomes and 15,708 genomes in v2.1.1 .
The core findings are directly supported by the reported numbers: ~4-fold variant increase but a near-halving of quality-pass proportion (76.6%β47.7%; joint callset 43.6%), PM2_Supporting rising 2.6-fold (81.5%β93.5%), and conflicting criterion combinations in <1% of shared variants .
Gene-level distribution (HBB 1,051; HBD 893; HBE1 879; HBA2 359; HBA1 316 quality-passing variants) plausibly reflects paralog callability rather than biology β the authors honestly attribute clustering patterns to short-read mapping limits in HBA1/HBA2 and HBG1/HBG2 . The 68.7% ancestry-exclusive variant figure is correctly framed as a sampling-depth artifact, not true population specificity β a stance consistent with gnomAD's own local-ancestry work showing that within-tract frequencies differ substantially from continental grpmax values .
Strengths: transparent methods (R 4.4.3, public scripts on GitHub), open data (gnomAD, ClinVar), declared absence of conflicts, EU/COST funding with no apparent product incentive, and careful avoidance of over-interpretation (descriptive ancestry exclusivity; no reclassification claims). Limitations the authors partly address: VCEP thresholds formally apply only to HBA1/HBA2/HBB, yet five additional genes were scored β the authors flag this but it remains a generalization risk. Concordance with ClinVar could be assessed for only 903/5,021 variants, and no globin-gene ClinVar record exceeds two stars because expert-panel review is absent, capping evidentiary resolution . Blind spots: no sensitivity analysis of how borderline variants near PopMax thresholds behave under resampling or confidence-interval uncertainty; the 'how to falsify' test (re-estimation in larger future releases flipping PM2_Supporting) is acknowledged implicitly but not quantified. Constraint-metric shifts (e.g., HBG2 LOEUF 0.83β0.38) are interpreted cautiously as partly technical, which is appropriate given known paralog artifacts in constraint estimation .
Verdict (moderate confidence): a methodologically sound, honestly caveated resource paper whose central conclusion β that population evidence is release- and callset-dependent and gnomAD version must be reported routinely β is well supported and practically actionable for diagnostic laboratories. What would change it: demonstration that VCEP threshold crossings are stable across future gnomAD releases despite the quality-pass decline, or evidence that joint-callset discordance (4.4% of shared variants) materially alters final classifications.
Know what changed, what holds up, and what remains uncertain. Every Friday. No ads.