Why HLA Hypersensitivity Studies Need Carefully Chosen Comparisons
Some drug hypersensitivity reactions are strongly associated with particular human leukocyte antigen (HLA) alleles. A large association, however, is not enough to establish how a study should be designed or what its results mean clinically. The central comparison is between people who developed a well-defined reaction after exposure to the implicated drug and people who had a sufficiently comparable opportunity to develop that reaction but did not.
The distinction matters because a person who never received the drug could not have developed a drug-induced reaction during that period. Comparing cases with such unexposed controls can make the groups differ in ways unrelated to susceptibility after exposure. It also does not directly answer the question of interest: among people who take this drug, how does risk differ by HLA genotype?
A study should specify the drug, the hypersensitivity phenotype, the exposure period, and the time window in which cases can occur. It should also specify the HLA allele at an appropriate level of resolution and the genotype contrast being estimated—for example, carriers of one named allele versus noncarriers. An association with an HLA allele is evidence of a statistical relationship in the studied setting; by itself, it does not establish that the allele is causal or that it predicts an individual outcome with high accuracy.
Define the Estimand Before Choosing Controls
The estimand is the quantity the study aims to describe or estimate. For a prospective exposed cohort, a natural estimand is the risk of the adjudicated hypersensitivity outcome among exposed carriers compared with exposed noncarriers over a stated follow-up period. A risk ratio compares these two risks; a risk difference subtracts them. Both require a denominator that represents the relevant exposed genotype group and adequate outcome follow-up.
A case-control study instead samples people according to outcome status: it selects cases and controls, then compares their allele frequencies. Provided the controls represent the exposure experience of the population that gave rise to the cases, this design can estimate an odds ratio for the allele–outcome association. The numbers of sampled cases and controls are set by the design, so the fraction of sampled participants who are cases is generally not the event prevalence in the exposed population.
Control selection is therefore not a minor implementation detail. Appropriate controls should have received the same implicated drug, or have a clearly justified exposure definition, and should have comparable opportunity to be observed through the relevant risk window. If controls are recruited from a group that is less likely to have received the drug, or is assessed over a shorter period, differences in exposure opportunity can distort the comparison. The outcome definition and case adjudication should be applied consistently, and the HLA genotype should be measured in a way that does not depend on whether someone is a case.
| Study question | Suitable design feature | What the result can estimate |
|---|---|---|
| How does risk differ by allele among people taking the drug? | Follow an exposed cohort with genotype-specific denominators and a defined outcome window. | Genotype-specific risks, risk ratio, and risk difference over that window. |
| Are cases more likely than comparable exposed controls to carry the allele? | Select cases and controls from the same relevant exposed population. | An allele–outcome odds ratio, subject to the sampling and comparability assumptions. |
| What fraction of carriers will experience the reaction? | Use a cohort or other data that identify the underlying event risk. | Positive predictive value for the specified population, exposure, and time period. |
Allele-Specific Risks and Case-Control Odds Ratios
Let \(A\) indicate carriage of the HLA allele of interest, and let \(Y=1\) indicate the defined hypersensitivity outcome during the risk window. In an exposed cohort, the risk among carriers is the number of carriers who become cases divided by the number of exposed carriers followed. The corresponding risk among noncarriers uses the noncarrier denominator. The risk ratio and risk difference are, respectively, the ratio and the difference of these risks.
Here \(a\) and \(c\) are the numbers of cases among exposed carriers and noncarriers, while \(b\) and \(d\) are the corresponding numbers without the outcome. A risk ratio above 1 indicates higher observed risk among carriers; a positive risk difference expresses the absolute excess risk over the specified period. These are estimates for the study population and exposure conditions, not automatically causal effects.
For a case-control study, arrange the sampled participants by outcome and allele status. The odds ratio is the odds of carrying the allele among cases divided by the odds of carrying it among controls; equivalently, it is the ratio of outcome odds in carriers to outcome odds in noncarriers in the source population under suitable sampling. For the table below, that estimate is \(ad/(bc)\).
| Allele carrier | Noncarrier | |
|---|---|---|
| Cases | \(a\) | \(c\) |
| Exposed controls without outcome | \(b\) | \(d\) |
The case-control odds ratio has an important sampling justification. Suppose cases are sampled with probability \(s_1\) and controls with probability \(s_0\), and within each outcome group these probabilities do not depend on allele status. Among carriers, the sampled case-to-control odds equal the population case-to-control odds multiplied by \(s_1/s_0\). The same factor applies among noncarriers. Dividing the two sampled odds cancels that factor, leaving the population odds ratio. If control selection depends on HLA status within outcome group, or controls do not represent the exposed population that generated the cases, this cancellation need not hold.
An odds ratio is not a risk ratio. When the outcome is uncommon in both genotype groups, the odds ratio may be numerically close to the risk ratio, but they are distinct quantities. A case-control sample also does not supply the population denominators needed to calculate absolute risks unless the sampling design and external information support that calculation.
Worked Example: Estimating an Allele Association With Exposed Controls
Consider hypothetical data from a case-control study of a defined hypersensitivity reaction after one drug. The 60 cases received the drug and met the prespecified phenotype definition. The 600 controls also received the drug, completed the relevant observation window, and did not meet the case definition. HLA testing identifies 36 carriers among cases and 90 carriers among controls. Assume selection of controls represents the exposed population that produced the cases, and that sampling within case and control groups does not depend on allele status.
There are \(60-36=24\) noncarriers among cases and \(600-90=510\) noncarriers among controls. The cross-product odds ratio is
The estimated odds of the defined reaction are 8.50 times as high for carriers as for noncarriers in the source exposed population, under the case-control sampling assumptions. This is not a statement that a carrier's risk is 8.50 times as high, nor does it provide the carrier's absolute probability of reaction.
For an approximate 95% confidence interval, the standard error of the log odds ratio is calculated from the four table counts:
Thus the interval on the log scale is \(\log(8.50)\pm1.96(0.2873)=2.140\pm0.563\), or approximately \((1.577,2.703)\). Exponentiating the endpoints gives an approximate 95% confidence interval of \((4.84,14.93)\) for the odds ratio. Under the model and sampling assumptions, this procedure has 95% coverage over repeated studies; it does not mean there is a 95% probability that this fixed parameter lies in this particular interval.
Worked Cohort Risk Estimates
Worked Example: Comparing Absolute Risks by Allele Status
In a separate hypothetical cohort, 10,000 people initiate the same drug and are followed for 12 months, with complete ascertainment of the defined hypersensitivity outcome. Five hundred are carriers of the HLA allele and 9,500 are noncarriers. There are 80 cases among carriers and 120 among noncarriers. Assume the exposure definition and follow-up window are applied consistently, and each person is counted once.
The estimated 12-month risks are
The estimated risk ratio is
The estimated risk difference is
In this hypothetical cohort, carriers had an estimated 12-month risk about 12.7 times the noncarrier risk, a difference of 14.74 percentage points. The cohort's overall event risk was \(200/10{,}000=2.00\%\). Because the study followed exposed people with genotype-specific denominators, it can estimate these risks for its specified population and period. The estimates still depend on adequate outcome ascertainment and may not transport unchanged to a population with different prescribing or follow-up conditions.
For comparison, the odds ratio from these cohort counts is \((80\times9{,}380)/(420\times120)=14.89\). It is larger than the risk ratio of 12.67, illustrating why the two measures should not be substituted for one another.
Low Prevalence and Predictive Value
A strong relative association can coexist with a low probability that any one carrier will experience the reaction. The probability of the outcome among carriers is the positive predictive value (PPV) of a carrier result for that outcome in the specified exposed population and time window. PPV depends not only on the relative association, but also on the event prevalence and how common the allele is. It is not a fixed property of the allele across all clinical settings.
This distinction is consequential for HLA-associated hypersensitivity because some serious reactions remain uncommon even in the genotype group at elevated risk. Many carriers may therefore tolerate the drug, while some cases still occur among noncarriers. A genotype can be strongly associated with risk yet be an imperfect predictor of who will develop the reaction. Absolute risks and predictive values are particularly important when considering the possible consequences of screening, avoiding a drug, or choosing an alternative treatment.
Worked Example: A Large Risk Ratio With Low Positive Predictive Value
Suppose a hypothetical cohort follows 100,000 exposed people through the same risk window. Five percent, or 5,000 people, carry the allele. There are 70 cases among carriers and 130 among 95,000 noncarriers. Assume all participants have adequate follow-up and the case definition is applied consistently.
The event prevalence is \(200/100{,}000=0.200\%\). The carrier risk is \(70/5{,}000=1.40\%\), and the noncarrier risk is \(130/95{,}000=0.1368\%\). Their ratio is
Despite this estimated risk ratio of about 10.2, the PPV of carrier status in this cohort is \(70/5{,}000=1.40\%\). In other words, 98.6% of carriers did not have the reaction during the observed window. Conversely, 130 of the 200 cases were noncarriers, so absence of the allele did not identify every person who would develop the outcome. These figures illustrate why a relative association alone cannot answer a clinical prediction question.
Interpretation, Limitations, and a Common Mistake
A defensible interpretation names the drug, the phenotype, the HLA contrast, the outcome measure, and the period or population to which the estimate applies. It also distinguishes a statistical association from a causal claim. HLA alleles can be correlated with other genetic variation, and allele frequencies and clinical settings can differ among populations. Study design and analysis should therefore consider whether the groups are comparable and whether the estimate is appropriate for the population where it will be used. An association alone does not establish a biological mechanism or prove that changing prescribing on the basis of genotype improves outcomes.
The most common design mistake is to use unexposed people as controls for a drug-specific reaction, then interpret the allele contrast as susceptibility among drug users. Those controls were not at risk of the reaction through the implicated exposure during that period. A related interpretation mistake is to report a case-control odds ratio as an absolute risk or as a risk ratio. The correction is to align controls with the exposed source population and use a cohort or valid external risk information when absolute risk is required.
Low event prevalence also has practical consequences for study planning: obtaining enough well-defined cases can be difficult, and an association estimate can be imprecise when the informative case count is limited. The present tutorial's central lesson is to keep the design question distinct from the prediction question. A case-control comparison can efficiently assess an allele association; estimating who will experience an event requires the relevant denominators and baseline risks.
Check Your Understanding
Answer each question in the context of a drug-specific hypersensitivity study.
- Why is a person who never received the implicated drug generally unsuitable as a control for estimating susceptibility to a reaction after that drug?
- A case-control study reports an odds ratio of 6. What additional information is needed before estimating the probability that an exposed carrier will develop the reaction?
- Explain why a large risk ratio can coexist with a low positive predictive value when the outcome is uncommon.
- In a cohort, what denominators are used to estimate risk among HLA carriers and noncarriers?
- State one reason why an HLA association, even if strong, does not by itself establish causation or guarantee useful individual prediction.