An ancestry coverage audit of 108 pharmacogenomic recommendations

Incoprea / ResearchResearch blog ↗

Working research · Version 2 · 30 October 2026

An ancestry coverage audit of 108 pharmacogenomic recommendations

108actionable gene–drug pairs in the snapshot, across 21 genes and 27 CPIC guidelines
46pairs (42.6%) have no computed ancestry-specific phenotype frequency in the audited API field
6 of 10genes with phenotype frequencies cover fewer than nine ancestry groups
3gene–ancestry cells have an allele-record support proxy below 100

Preprint draft · not peer reviewed

Pharmacogenomic guidance turns a genotype into a prescribing decision. The frequencies that make that guidance operational come from particular study populations. This audit asks a narrow, checkable question of one public database: for how many recommendations does the CPIC API actually return a computed phenotype frequency, and for which ancestry groups?

A descriptive database audit, not medical advice. It records which fields are populated in a dated snapshot of a public API. It is not a systematic review, it does not assess whether any guideline is clinically adequate, and it does not estimate anyone’s individual risk.

1. Forty-six of 108 pairs lack a computed phenotype frequency

We retrieved the CPIC public API on 30 October 2026 and selected every pair at evidence level A, A/B or B with usedforrecommendation = Yes: 108 pairs, 21 genes, 27 guidelines. For each, we asked whether the gene carried a value in CPIC’s computed phenotype-frequency field.

For 62 of the 108 pairs, a phenotype frequency was present for at least one ancestry group. Forty-five pairs had allele-frequency records but no computed phenotype-frequency value, and one pair — ivacaftor and CFTR — had neither.

These categories describe which fields are populated in this snapshot. They do not measure the completeness of published evidence. CPIC recommendations can key on phenotype, on an activity score, or on allele status, and an absent phenotype field is not evidence that a guideline is unsupported or that a frequency has never been measured.

Two bar charts. Left: 62 pairs have phenotype frequencies, 45 have allele frequencies only, 1 has neither. Right: pairs with phenotype frequencies by ancestry group range from 52 for the American group to 62 for European, East Asian and Central/South Asian groups.
Figure 1. Which frequency layer the audited API returns for each pair, and the number of pairs with a phenotype frequency by ancestry group. Both panels count API fields, not the state of the literature.

Where phenotype frequencies were present, they were often incomplete. Six of the ten genes with any phenotype frequency covered fewer than nine ancestry groups:

GeneGroups with a phenotype frequencyMissing
NUDT154 / 9African American, Sub-Saharan African, Near Eastern, Oceanian, American
NAT25 / 9African American, Near Eastern, Oceanian, American
CYP3A57 / 9Oceanian, American
TPMT7 / 9Oceanian, American
CYP2B68 / 9American
SLCO1B18 / 9Latino

This is an observation about the snapshot. It does not show that NUDT15 has never been measured in those groups, and it does not support ancestry-based dosing decisions.

2. Where phenotype frequencies exist, they differ across groups

For genes with phenotype frequencies we summed the categories that CPIC’s own ehrpriority field marks for clinical attention. We call this the CPIC high-priority proportion. It is a clinical-attention label, not a calibrated probability of harm, and it aggregates categories that can move in different directions — an ultrarapid and a poor metaboliser are both flagged, for opposite reasons.

Sixty-one pairs had enough data to rank. The median difference between the lowest and highest ancestry group was 16.1 percentage points, and 25 pairs exceeded 20 points. Those 25 pair-level results are not 25 independent biological contrasts: they repeat five gene-level estimates — CYP2C19 (14 pairs), SLCO1B1 (7), CYP2B6 (2), ABCG2 (1) and CYP3A5 (1).

Heat map of the CPIC high-priority proportion for nine genes across nine ancestry groups, ranging from about 1 to 96 percent. Grey cells have no published frequency data; orange borders mark cells whose allele-record support proxy is below 100.
Figure 2. CPIC high-priority proportion by gene and ancestry group. Grey cells have no phenotype frequency in the snapshot. Orange borders mark cells whose support proxy is below 100; that proxy is described in section 3.
GeneLowest groupHighest groupGap (points)
CYP2C19American · 37.2%Oceanian · 96.5%59.3
CYP3A5European · 14.3%Sub-Saharan African · 72.9%58.6
ABCG2Sub-Saharan African · 1.2%East Asian · 52.0%50.8
CYP2B6East Asian · 37.0%Oceanian · 87.0%50.0
SLCO1B1Oceanian · 0.0%*American · 42.2%*42.2
NUDT15European · 1.4%East Asian · 20.8%19.4
CYP2C9East Asian · 5.6%Near Eastern · 21.7%16.1
CYP2D6Sub-Saharan African · 6.7%Oceanian · 17.7%11.0
TPMTCentral/South Asian · 3.7%African American · 13.5%9.8

*The SLCO1B1 value is a composite: the American 42.2% is 36.48% Decreased Function plus 5.76% Poor Function, not Decreased Function alone. Both SLCO1B1 cells carry a support proxy below 100 — see section 3. CYP2C19, not SLCO1B1, is the largest contrast in the table.

3. The support proxy, and what it cannot tell us

For each gene × ancestry cell we recorded the largest subjectcount among the allele-level records contributing to it. This is a sample-support proxy. It is not a demonstrated denominator for the phenotype estimate, not a count of unique people, and not an effective sample size.

Three of 75 cells had a proxy below 100, all for SLCO1B1: African American/Afro-Caribbean 38, Oceanian 28, American 64. The Oceanian cell is 0.0% and the American cell 42.2%.

It is tempting to call that contrast an artefact. This audit cannot support that. We did not reconstruct a denominator, propagate uncertainty, or examine the contributing cohorts, so we cannot say whether the contrast reflects sampling variation, cohort composition, modelling assumptions or a population difference. It requires cautious interpretation and further validation — a conclusion that applies to every cell here, not only the small ones.

4. HLA carrier frequencies

For genes where the actionable unit is a single allele, the allele table is the relevant measure. We estimated carrier frequency under Hardy–Weinberg as 1 − (1 − p)².

HLA-B*57:01, the abacavir hypersensitivity marker, ranges from about 1.2% of African American/Afro-Caribbean and Sub-Saharan African groups to 13.2% of Central/South Asian groups. HLA-B*15:02 ranges from approximately 0% to 8.9%, HLA-B*58:01 from 1.5% to 11.7%, and HLA-A*31:01 from 1.6% to 12.0%.

Four bar charts showing Hardy-Weinberg carrier frequency estimates by ancestry for HLA-B*57:01 (about 1 to 13 percent), HLA-B*15:02 (0 to 9 percent), HLA-B*58:01 (1.5 to 12 percent) and HLA-A*31:01 (1.6 to 12 percent).
Figure 3. Estimated population carrier frequencies under Hardy–Weinberg. These are prevalence estimates for a genotype in a population, not the risk of a reaction in a person whose genotype is already known.

The research map

Each row is an actionable gene, each column an ancestry group, and each cell records what the snapshot contains. Dark means a computed phenotype frequency. Mid-blue means allele records only. Pale means neither.

Grid of 21 actionable genes against nine ancestry groups. Four genes have phenotype frequencies in all groups. NUDT15 is missing in five groups, NAT2 in four. Ten genes have allele-level records only, shown in mid-blue. CFTR has neither.
Figure 4. Coverage of the audited API fields. Pale cells are places where this snapshot provides no frequency of either kind.

The coverage ledger

Every actionable gene, graded by what the audited fields contain. The tier is a summary of data availability only — it is not a judgement about evidence quality, clinical importance, or the soundness of any guideline.

GenePairsField coverageGroupsTier
ABCG21Phenotype9/9A
CYP2C1914Phenotype9/9A
CYP2C911Phenotype9/9A
CYP2D619Phenotype9/9A
CYP2B62Phenotype8/9B
CYP3A51Phenotype7/9B
TPMT3Phenotype7/9B
SLCO1B17Phenotype8/9B−
NAT21Phenotype5/9C
NUDT153Phenotype4/9C
CACNA1S7Allele only5/9D
CYP4F21Allele only8/9D
DPYD2Allele only6/9D
G6PD8Allele only5/9D
HLA-A1Allele only9/9D
HLA-B6Allele only9/9D
MT-RNR111Allele only5/9D
RYR17Allele only6/9D
UGT1A11Allele only8/9D
VKORC11Allele only8/9D
CFTR1Neither0/9E

Tier A: phenotype frequency in all nine groups. B: seven or eight. C: six or fewer. D: allele records only. E: neither. The ledger is downloadable as gene_reference.csv.

What this audit supports — and what it does not

It supports

  • The statement that 46 of the 108 selected pairs lack a computed ancestry-specific phenotype frequency in this snapshot, and that six of ten phenotyped genes cover fewer than nine groups.
  • The statement that where phenotype frequencies exist, the CPIC high-priority proportion varies across groups — with the caveat that 25 pair-level spreads repeat five gene-level estimates.
  • The statement that three gene–ancestry cells carry a support proxy below 100, and that their contrasts should be read with caution.

It does not support

  • Any claim that 43% of clinical recommendations lack evidence. The audit measured fields in one API, not the literature.
  • Any claim that a particular contrast is noise, an artefact, or a population difference. Uncertainty was not propagated.
  • Any claim about prescribing, testing policy, outcomes, cost-effectiveness or research priority. None of those were measured.
  • Any ancestry-specific treatment recommendation. Genotype prevalence is not harm, and this is not medical advice.

Limits

  • One source, one snapshot. CPIC only; DPWG and others publish overlapping guidance. The API can change between releases.
  • Fields, not evidence. An empty phenotype field does not mean guideline authors lacked data, nor that allele-level guidance is unsupported.
  • Coarse categories. CPIC’s ancestry groups are administrative groupings that pool heterogeneous populations and do not describe genetic ancestry.
  • Derived frequencies. CPIC computes phenotype frequencies from allele counts under modelling assumptions; we report them as published and did not propagate uncertainty.
  • Attention is not risk. ehrpriority flags a result for clinical attention. It is not a calibrated probability of harm, and the composite mixes directions of effect.
  • Repeated estimates. The 61 rankable pairs are not 61 independent comparisons; they rest on nine informative gene-level estimates.
  • No external replication. We did not cross-check against gnomAD or 1000 Genomes. This version contains no independent population-genomic validation.

Method

We queried api.cpicpgx.org/v1/ on 30 October 2026 and selected pairs at level A, A/B or B with usedforrecommendation = Yes. For each gene we recorded whether gene_result.frequency was populated, whether population_frequency_view contained allele records, and the largest subjectcount among those records. High-priority phenotypes are those whose ehrpriority contains “High Risk”; their frequencies were summed per group. HLA carrier frequencies use 1 − (1 − p)². The supplied analysis reports phenotype-frequency sums approximately one within each gene–ancestry cell — 71 of 75 within 0.001, with a minimum of 0.996. These source-reported checks test internal consistency, not correctness; the publisher did not independently reconstruct them from the raw API snapshot. Definitions and provenance notes accompany the processed-data download. Additional working pages are being assembled in the research library.

Sources

Primary resource: CPIC public API, snapshot date reported as 30 October 2026. See the official data-model documentation, CPIC content license and usage notes, and the related 2024 biogeographic-frequency study. No CPIC endorsement is implied.

Data

Cite this

Mikle T. An ancestry coverage audit of 108 pharmacogenomic recommendations.
Preprint draft, not peer reviewed. Version 2, 30 October 2026.
https://incoprea.com/research/pharmacogenomics-coverage/
Revision note, version 2 (30 October 2026). Corrected coverage wording so that the 43% figure refers to a missing computed phenotype-frequency field, not to absent ancestry data or absent evidence. Relabelled support_n as an allele-record support proxy rather than a phenotype denominator. Removed the conclusion that the SLCO1B1 contrast is an artefact or sampling noise; the audit does not propagate uncertainty. Corrected the SLCO1B1 42.2% figure as a composite of Decreased Function (36.48%) and Poor Function (5.76%). Renamed the aggregate metric from “high-risk” to “CPIC high-priority proportion”. Removed novelty and clinical-priority claims. Added the preprint status. Version 1 is archived and unchanged.

Research, analysis and drafting were prepared with DeepSeek Harness. Internal consistency of the supplied processed tables was checked during preparation. No independent clinical review has been performed. This descriptive database audit is not medical advice.