HLA-B*14:01 and HLA-B*35:01 are associated with trimethoprim-sulfamethoxazole induced liver injury

Abstract


Introduction
Trimethoprim-sulfamethoxazole (TMP-SMX) is a fixed combination synthetic antimicrobial that is commonly used to treat various bacterial infections and as a prophylaxis against opportunistic infections (1) .TMP-SMX use has been associated with life-threatening and presumed immune-mediated idiosyncratic adverse drug reactions such as Stevens-Johnson syndrome (SJS), toxic epidermal necrolysis (TEN), drug-reaction with eosinophilia and systemic symptoms (DRESS) , severe drug induced liver injury (DILI), and blood dyscrasias (2) .TMP-SMX is among the top 5 leading causes of DILI in the United States (3) .The risk factors for liver injury associated with TMP-SMX are largely unknown, other than that they are more common among African-Americans and those infected with human immunodeficiency virus (1,4) .The pattern of liver injury associated with TMP-SMX is typically cholestatic or mixed, although severe hepatocellular injury and acute liver failure have been reported as well (5) Immunogenetic factors have been increasingly shown to play an important role in the pathogenesis of severe immune mediated adverse drug reactions including DILI (6) .Previous

Study Participants:
Individuals aged ≥ 2 years of age with suspected DILI and meeting predefined eligibility criteria were enrolled into the DILIN Prospective (NCT00345930) and Retrospective (NCT00360646) studies.The designs of these two studies have been previously published (13, 14) .In brief, written informed consent were obtained from individuals with acute liver injury suspected due to prescription, over-the-counter, or herbal and dietary supplements prior to enroll at multiple centers throughout the United States.Our study followed the ethical guidelines of the 1975 Declaration of Helsinki and was approved by the Institutional Review Board of each study site.Upon recruitment, participants underwent systematic evaluation for competing etiologies and all cases were adjudicated for causality by structured expert consensus.The causal relationship between suspected medication(s) and the liver injury episode was categorized into definite (>95% likelihood), highly likely (76-95% likelihood), probable (51-75% likelihood), possible (26-50% likelihood), and unlikely (≤25% likelihood) (14).Cases included in this study were adjudicated to have definite, highly likely or probable DILI due to TMP-SMX, referred to as high confidence DILI cases.For this study's purposes, instead of self-reported race and ethnicity, we utilized the genetic ancestry inferred by genome-wide single nucleotide polymorphism (SNP) data from DILIN and 1000 genomes project (15) (Supplementary Materials).Of 1,916 DILI patients, 86 (4.5%) were suspected to have DILI due to TMP-SMX with 72 patients classified as high confidence cases.Based on genetic inferred ethnicity, these patients were then grouped to 51 European Americans, 10 African Americans, 6 Hispanics, 4 Asians, and 1 other (Supplementary materials Section I, Table S2).In this study, the 51 European American and 10 African American DILI cases were analyzed.Participants in this paper have been included in other publications from the DILIN (16) .

Study design and HLA data:
We utilized a case-control study design to identify HLA alleles that may play a role in the risk of developing TMP-SMX DILI in European Americans and African Americans, respectively.Due to the lack of control data in DILIN, considering the low incidence rate of DILI in the population, we used large population cohorts as controls to overcome the potential bias of including potential subjects with DILI events.With the approval from NCBI dbGaP, we obtained the GWAS data from eMERGE-I: Genome Wide Association Studies of Network Phenotypes (phs000360.v3.p1) and PAGE: The Charles Bronfman Institute for Personalized Medicine BioMe BioBank (phs000925.v1.p1).The eMERGE-1 GWAS dataset consisted of 17,045 subjects who were predominantly Americans with European ancestry, while the PAGE GWAS dataset consisted of 12,932 subjects with diverse racial backgrounds.Genotyping platforms were Illumina Human 660W Quad for eMERGE-1 and MEGA for PAGE.Similarly, we used genome-wide SNP data to infer ethnicity for each dataset.This led to 12,156 European Americans from eMERGE-I and 5,439 African Americans from PAGE to serve as controls.
For TMP-SMX DILI cases, HLA alleles were determined from the deep sequencing of HLA Class I and II genes using Illumina MiSeq (details in Supplementary Materials).For population controls, we utilized the GWAS SNPs to impute HLA alleles in HLA Class I and II genes for European Americans and African Americans, respectively.Specifically, fourdigit HLA alleles were imputed based on the SNPs in the MHC region in chromosome 6 using the HLA genotype imputation with attribute bagging (HIBAG) program (17) .All analyses were conducted for European Americans (51 cases vs. 12,156 controls) and African Americans (10 cases vs. 5,439 controls), respectively.

Statistical Analysis:
Descriptive statistics for demographic and patient characteristics were computed for 51 European American and 10 African American TMP-SMX DILI cases, respectively.We presented frequency (percentage) for categorical variables and median (interquartile range) for continuous variables.Allele frequencies (AF) and carriage frequencies (CF, proportion of patients carrying the allele) of all four-digit HLA alleles were computed for cases and controls, respectively, for each ethnic group.To evaluate the accuracy of imputed HLA alleles in controls, we compared the AF of imputed HLA alleles to population AF available from BeTheMatch (18) and AFND (19) for European Americans and African Americans, respectively.We set an arbitrary threshold to exclude imputed HLA alleles that have an absolute AF difference from population AF greater than 0.05.Since HLA alleles were determined by all variants in the coding sequence of the gene, this reduced variants to one marker with multiple alleles for each HLA gene (e.g., HLA-A*, HLA-B*, etc).Therefore, association tests were performed for a single allele at a time.Considering rare alleles will have low statistical power in single-allele association test, we excluded ultra-rare HLA alleles that have either allele count < 2 in TMP-SMX DILI cases or < 10 in controls.
For allelic association tests, each four-digit HLA allele was coded as 0, 1, or 2 based on the number of the allele an individual carried.The Fisher's exact test was used to test AF differences between cases and controls for each HLA allele.For the European American subset, multivariable logistic regression with Firth penalization for rare events (Firth logistic regression) was used to test the association between each HLA allele and DILI with covariate adjustment.For covariate selection, univariate Firth logistic regression was performed for age, sex, and PC1 to PC10 individually first.We selected the top three variables meeting p<0.05 as covariates.Furthermore, since eMERGE-1 released age data in decades, we categorized age to three levels (1: < 40 years; 2: 40-69 years; 3: >69 years).False discovery rate (FDR) by Benjamini-Hochberg (20) was computed to correct multiple testing.Significant HLA alleles were determined based on FDR < 0.15 from the multivariable Firth logistic regression analysis.To delineate the relationship among the significant HLA alleles, we conducted conditional analyses by adjusting the most significant HLA allele along with other covariates.For the African American dataset, due to 10 cases only, we did not pursue multivariable regression analysis.Instead, we conducted Fisher exact test only and highlighted alleles meeting p<0.05.
Several secondary analyses were performed to understand the top HLA alleles further.First, haplotype association was tested for two and three gene combinations using the score statistics in the haplo.statsprogram (21) .Second, using the NetMHCpan-2.8algorithm in MHCcluster (22), HLA alleles in the gene of interest were clustered based on their predicted peptide-binding specificity.We then tested the association between each allelic cluster and TMP-SMX DILI.That is, HLA alleles in the same cluster were coded as the same marker and then the same association analysis as described above were followed.Third, we checked if the HLA allele of interest is present in persons with high confidence DILI due to nonantibacterial sulfonamides enrolled into the DILIN studies.Fourth, to evaluate if the HLA allele of interest is specific to TMP-SMX subset, we compared CF of the HLA allele of interest among DILIN drug groups that had at least 10 high confidence DILI cases.Finally, we tested if the HLA allele of interest is associated with the any-cause DILI for the corresponding ethnicity.

Replication of top HLA alleles:
We were able to identify five European Americans with TMP-SMX and three patients with TMP induced liver injury enrolled in the iDILIC consortium to evaluate our findings in DILIN.Causality assessment was performed as previously reported (23) .The cases were previously analyzed as a part of the broad GWAS on susceptibility to DILI (24) .

Amino Acid association analysis:
We first obtained the AA sequence alignment of all available HLA alleles in the same HLA gene from the Immuno Polymorphism Database (25) .We then identified a set of polymorphic AA positions within the binding site.Using the HLA allele of interest as the reference allele, we converted the genotypes of the HLA allele to AA residue pairs of all polymorphic AA positions based on if the reference AA residue is present.In other words, assuming K polymorphic AA positions identified for a HLA gene, HLA genotype data of each subject is recoded as a vector of K markers with 0, 1, or 2 based on the number of reference AA residues present at the AA position.The same association tests, as described above for HLA alleles, were performed for each polymorphic AA position for European Americans and African Americans, respectively.

Molecular Docking:
Due to the lack of crystal structure of HLA-B*14:01, an atomic model of HLA-B*14:01 was generated based on the crystal structure of HLA-B*14:02, 99.6% identical, PDB code 3BVN (26) using SWISS-MODELLER (27) .The HLA-B*14:01 model and HLA-B*35:01 crystal structure, PDB 6BJ8 (28) , were used for molecular docking.AutoDock Tools was used for molecular docking of trimethoprim and sulfamethoxazole with AutoDock Vina (29) (e.g., add H, generate charges for each atom to be used in scoring).Scoring grids were 20 × 20 × 20 Å, and centered on 3 sites within the antigen binding clefts of HLA-B*14:01 and HLA-B*35:01, including sites corresponding to the Cα of the first peptide position (P1), Cα of the fifth peptide position (P5, the middle of the antigen binding cleft), and Cα of the terminal peptide position (P9).TMP and SMX were docked with exhaustiveness set to 10.The top 8 scoring orientations were output and compared.Sulfonamide was manually positioned in the modeled antigen binding cleft of HLA-B*14:01 with PyMol, which was used to generate molecular graphics (The PyMOL Molecular Graphics System, Version 1.8 Schrödinger, LLC.).

Results
From September 2004 through January 2019, the DILIN has enrolled a total of 2,256 participants with suspected DILI and 1,916 of them had undergone causality adjudication and had their HLA genes sequenced (Figure 1).There were 86 participants with suspected DILI in whom TMP-SMX was the primary implicated agent.Upon adjudication, 72 were deemed to have high-confidence DILI events (causality scores of definite, highly likely, or probable) with 51 European Americans and 10 African Americans.Table 1 describes the summary statistics for selected clinical characteristics of European Americans and African Americans, respectively.The median latency was longer in European Americans (25 days) than in African Americans (14.5 days).While the pattern of liver injury was evenly distributed among hepatocellular, cholestatic, and mixed categories in European Americans, it was predominantly hepatocellular in African Americans.Eleven patients (22%) of European Americans had peripheral eosinophilia but none appeared in African Americans.While a cutaneous rash was commonly reported (European Americans 47%, African Americans 40%), there were no reported instances of Stevens-Johnson syndrome in our TMP-SMX DILI patients.
We first examined the quality of imputed HLA alleles in controls: the eMERGE-1 and PAGE data.We identified four class II alleles in European Americans and seven class II alleles in African Americans with AF deviating from population AF more than 0.05 (Supplementary Table S3).These alleles were excluded from the association analysis.In addition, 114 ultrarare HLA alleles in European Americans and 130 in African Americans were excluded.This led to 85 alleles in European Americans and 26 alleles in African Americans remained for interrogation in this study.
Given the rarity of the HLA-B*14:01 allele in the European American general population (population AF=0.008), we did not expect any carrier in the five TMP-SMX DILI cases from iDILI.However, instead, we found that one TMP-SMX DILI case carried the HLA-B*14:01 allele, a higher AF than expected (AF=0.1).Interestingly, this case with HLA-B*14:01 also carried HLA-C*08:02.None of the five TMP-SMX iDILI cases carried HLA-A*34:02 or HLA-B*27:02.Further, among three TMP cases in iDILI, they did not carry HLA-B*14:01 or HLA-B*14:02.
In the African American subset, among the same 45 AA residues in the binding sites, there were 26 polymorphic AA residues in the available HLA-B alleles.The AA residues at positions 63 (Asn 63 : Asparagine) and 67 (Phe 67 : Phenylalanine) in HLA-B*35:01 reached the nominal significance threshold (Table 4).Phe 67 , present in HLA-B*35:01, HLA-B*08:01, and HLA-B*51:01, ranked at the top with a two-fold higher frequency in cases than in controls (frequency: 0.55 vs. 0.25, p=0.007).

Molecular docking of Sulfamethoxazole or Trimethoprim with B*14:01 and B*35:01:
For HLA-B*14:01, molecular docking predictions of non-covalent interactions (H bonds and van der Waal contacts) showed that neither SMX nor TMP were predicted to bind with high affinity to the antigen binding cleft of HLA-B*14:01.The estimated ΔG values were −6.7 kcal per mole for SMX and −6.4 kcal per mole for TMP (<−10 kcal per mole suggests high affinity binding, −6.5 to −9.0 moderate affinity, >−6.5 kcal per mole, low affinity).However, drug metabolites and peptides may form intermolecular contacts with HLA-B*14:01 by covalent interactions.It is notable that HLA-B*14:01 exhibited a chemically reactive side chain located at a site in the antigen binding cleft critical for peptide binding: a cysteine at position 67 (Cys 67 ), the residue detected in our amino acid analysis.As the cysteine sulfhydryl group is often in a reactive thiolate form in cells, free unpaired cysteine residues have the potential to form covalent bonds with reactive groups (e.g., with peptides that contain cysteine, or other molecules with reactive groups).SMX is known to be metabolized into reactive forms, such as sulfonamide, that have the potential to bind Cys 67 in the antigen binding cleft of HLA-B*14:01.These covalent bonds may partially explain the association between TMP-SMX and HLA-B*14:01, as metabolite binding at this site would result in altered peptide presentation to T cells.Based on the location of the reactive group in sulfonamide and the free sulfhydryl of Cys 67 in the HLA-B*14:01 model, SMX was accommodated in a site comprised of residues associated with TMP-SMX DILI (Cys 67 , Trp 97 and Thr 163 , Figures 3A and 3B).These data suggested functional interactions between HLA antigen binding clefts and reactive SMX metabolites.
For HLA-B*35:01, molecular docking predicted SMX to bind better than TMP (estimated ΔG = −7.6 kcal/mol for SMX, versus ΔG = −6.9kcal/mol for TMP) (Figures 3C and 3D).SMX was predicted to form hydrogen-bond interactions with residues in the antigen-binding clefts of HLA-B*35:01 that are identical in HLA-B*15:02, a previously identified HLA allele associated with TMP-SMX induced severe cutaneous adverse drug reactions such as SJS/TEN in a Thai population (12) .

Patient characteristics for HLA-B*14:01 and HLA-B*35:01 carriers and non-carriers:
We observed slightly longer latency of liver injury (median 29 vs. 21) and higher proportion of peripheral eosinophilia (40% vs. 16%) in TMP-SMX DILI patients with HLA-B*14:01 than those without the allele, but the differences were not statistically significant.No significant differences were found for liver biochemistries at onset and at peak, pattern of liver injury, recovery of liver biochemistries, severity and outcomes either (Supplementary Table S7).Similarly, the characteristics and the outcomes of African Americans with and without HLA-B*35:01 were not significantly different from each other.

Discussion
Our study found that HLA-B*14:01 was associated with TMP-SMX related DILI in European Americans with the risk being 9.2-fold higher among individuals carrying the allele than those without.The HLA-B*14:01 allele was about 5.5-fold more frequent in European American patients with TMP-SMX DILI than controls.This large frequency difference was also present in one of five TMP-SMX DILI validation cases carrying HLA-B*14:01 from the iDILIC cohort.For the African American subset, HLA-B*35:01 ranked at the top as a potential risk factor for TMP-SMX DILI.While only 10 African American TMP-SMX DILI patients are available in our dataset, half of them carried HLA-B*35:01, a much higher allele and carriage frequencies than in the control group.Although replication studies with a larger sample size are needed, our results are promising and consistent with previous genetic studies of drug-specific DILI that identified HLA Class I and II genes (8,30,31) as genetic risk factors for DILI.
In the European American subset, HLA-A*34:02, HLA-B*27:02, and HLA-C*08:02 were the other top alleles that showed a significant association with TMP-SMX DILI.However, we considered HLA-B*14:01 as the most promising risk factor for the following reasons.First, HLA-B*14:01 remained significant after correcting for multiple testing.Second, both HLA-A*34:02 and HLA-B*27:02 were present in only two TMP-SMX DILI patients.
Together with their low frequency in controls, the estimated effect sizes (OR) were unstable, as shown by the large 95% CI.Therefore, further replication is definitely needed.Third, while HLA-C*08:02 is the most common allele among these four top alleles, HLA-C*08:02 is known to be in strong linkage disequilibrium with B*14:01/B*14:02 (32) .In our data, all eight HLA-C*08:02 carriers carried either HLA-B*14:01 (five patients) or HLA-B*14:02 (three patients).Since HLA-B*14:02 did not show association with TMP-SMX DILI, the association signal of HLA-C*08:02 is more likely dependent on HLA-B*14:01.As we demonstrated in the conditional analysis, after adjusting for the effect of HLA-B*14:01, the significant effect of HLA-C*08:02 on TMP-SMX DILI was diminished.
Interestingly, within the TMP-SMX DILI subset, we observed an opposite phenomenon where HLA-B*14:01 had higher frequency than HLA-B*14:02 (AF=0.049 vs. 0.029).Further, the HLA-B*14:02 frequency was similar between cases and controls.This implied that the effect of HLA-B*14:01 was independent from HLA-B*14:02.To date, we have not found any reports linking HLA-B*14:01 to DILI in general or to drug-specific cases.Among all causal drugs in the DILIN cohort, HLA-B*14:01 was more prevalent in the TMP-SMX subset than in other drug-specific subsets or the subset caused by non-antibacterial sulfonamide like drugs.Therefore, if the effect HLA-B*14:01 on DILI risk we observed in this study is correct, it is specific to TMP-SMX only (and likely to sulfamethoxazole).
In the literature, the HLA-B*14:01 allele has been linked to protection against HIV progression and with lower plasma viral levels (pVL) in HIV patients (32,33) .Different AA residues at a single AA position can also a play role (protective or risk) in disease outcomes as illustrated in HIV related studies (33) .In the current study, the main identified risk allele, HLA-B*14:01, has an unpaired cysteine at position 67.This provides clues not only to a shared epitope with other HLA alleles that could in theory be implicated in risk of DILI but also a clue to the structure of the reactive metabolites and specific groups of SMX such as sulfonamide RS(O)NR' 2 that might be expected to covalently bind with host proteins.
In the smaller group of African American persons with TMP-SMX related DILI, HLA-B*35:01 appeared to be an important HLA risk allele.HLA-B*35:01 was previously reported to be associated with DILI caused by Polygonum multiflorum, a Chinese herbal medicine (34) , in Chinese patients.However, the relationship between HLA-B*35:01 and DILI with other drugs has not been reported.HLA-B*35:01 is known to share peptidebinding specificities with HLA-B*51:01 and HLA-B*53:01.Our allelic cluster analysis also showed that the cluster of these three alleles had stronger association with TMP-SMX DILI than HLA-B*35:01 alone.The two significant AA residues in positions 63 (Asn 63 ) and 67 (Phe 67 ) identified are also present in HLA-B*51:01 and HLA-B*53:01.Therefore, it is likely that these three alleles are equally important, but this conclusion deserves further study in large cohorts of African Americans with DILI.On the other hand, molecular docking led us to HLA-B*15:01, which shares nearly identical peptide binding preference as HLA-B*35:01.Since HLA-B*15:01 was reported to associate with SJS/TEN due to TMP-SMX in Thai population, the similarity in the antigen-binding cleft between these two alleles offer potential clues to explain the role of HLA-B*35:01 in TMP-SMX DILI in African Americans.
Despite these comprehensive analyses suggesting HLA-B*14:01 and HLA-B*35:01 as risk factors for TMP-SMX induced DILI in European Americans and African Americans, respectively, the limitations in this study need to be considered.First, the numbers of cases were limited, a common problem in most drug-specific genetic association studies.Using our European American subset for a post-hoc power calculation, given 51 cases and 12,156 controls and assuming AF of the risk allele at 0.009 (same as HLA-B*14:01) and complete LD between the risk allele and the testing allele, we have 80% power to detect genetic relative risk at 5.9 for an allele at the same AF.This calculation implies that our sample size was sufficient to detect HLA-B*14:01 but will require more samples to detect other HLA risk alleles with smaller effect sizes.In addition, because of this sample size limitation, several approaches were used to enhance the analysis quality, including using large population control datasets for comparison, Fisher exact tests to accommodate the low counts, and logistic regression models with Firth penalization to reduce the potential bias of parameter estimates.With these strategies, we expected a reduced chance of biased findings.Second, while HLA genotype data in DILI patients were obtained from advanced sequencing techniques, the HLA data in population controls (eMERGE-1 and PAGE) were obtained fromimputation of SNPs in the MHC region.This may raise concerns if HLA calls are inaccurate in the control group.The imputation quality by various programs have been compared before.Karnes et al. (35) reported 97.6% concordance rate for European Americans and 92.9% concordance rate for African Americans using HIBAG.Particularly, alleles in Class I genes had over 98% concordance rates.Therefore, with the large sample size of the control groups used in this study, the impact of minor imputation errors on the result should be small, particularly for the Class I genes.Third, binding interactions between SMX and TMP were predicted, rather than being experimentally determined.Moreover, predictions of SMX and TMP contact with HLA-B*14:01 were based on an atomic model (as opposed to a solved crystal structure).However, HLA-B*14:01 is nearly identical to the crystal structure of HLA-B*14:02, with only one difference in α1 and α2 domains, residues 1-180, that form the antigen binding cleft: (Ser 11 vs. Ala 11 ).Since the side chain at this position is not oriented toward peptide ligands (in class I HLA molecules, the side chain at position 11 is oriented towards the α1 helix), this polymorphic difference between HLA-B*14:01 and -B*14:02 is unlikely to influence reduce confidence in modeling the antigen binding cleft of HLA-B*14:01.Overall, the results from our molecular docking may to help us understanding interactions with SMX and its metabolites.
Given the rarity of liver injury due to TMP-SMZ and the relatively small effect size for the associated HLA alleles, our observations are not necessarily immediately clinically actionable from the perspective of pre-prescription screening to avoid a clinical event.However, our observations may help to both inform the immunopathogenesis of DILI related to trimethoprim-sulfamethoxazole and also help in the risk stratification and early diagnosis of such patients.Indeed, in future state where integrated health care systems (e.g., Geisenger) undertake whole exome or whole genome sequencing of all their patients, our genetic association findings may assist in the thorough evaluation of risk -benefit ratios and help promote safe prescription practices.

Supplementary Material
Refer to Web version on PubMed Central for supplementary material.

Acknowledgements:
We are extremely grateful to all patients and study coordinators participated in the DILIN study.We also thank NCBI/dbGaP at http://www.ncbi.nlm.nih.gov/gapapproved our application to obtain control datasets.The data of non-DILI control subjects were from eMERGE-I and PAGE GWAS studies from dbGaP.The eMERGE network is a consortium of five participating sites (Group Health Seattle, Marshfield Clinic, Mayo Clinic, Northwestern University, and Vanderbilt University) funded by the National Human Genome Research Institute (NHGRI).Funding supports for the cohort and genotyping of eMERGE samples are as the following: (1) Group Health Seattle: U01AG006781 and U01HG004438 with genotyping performed at Johns Hopkins University; (2) Mayo clinic: UOIHG004599 and U01HG004424 with genotyping performed at The Broad Institute; (3) Marshfield Clinic: U01HG004608 and U01HG004438 with genotyping performed at Johns Hopkins University; (4) Northwestern University: U01HG04603 and U01HG004424 with the genotyping performed at The Broad Institute; (5) Vanderbilt University: U01HG004603 and U01HG004424 with the genotyping performed at The Broad Institute.The eMERGE GWAS datasets were from dbGaP accession number phs000360.v3.p.PAGE dataset was from the Charles Bronfman Institute for Personalized Medicine (IPM) BioMe BioBank at the Icahn School of Medicine at Mount Sinai (New York).Phenotype data collection was supported by The Andrea and Charles Bronfman Philanthropies.Funding support for genotyping, which was performed at The Center for Inherited Disease Research (CIDR), was provided by the NIH (U01HG007417).The GWAS data were from dbGaP accession number phs000925.v1.p1.Disclosures: Chalasani and Watkins have consulting agreements and research grants from several pharmaceutical companies.Nicoletti is also an employee of Sema4.They all declare that their engagement does not have direct or significant relationship to this paper.Predicted H bond interactions are shown as black dashed lines.The position 67 (Phe 67 ), which we detected in our AA analysis, was not predicted to form contact with SMX by molecular docking.

Figure 1 :
Figure 1: Consort diagram to outline the sample sizes used in this study

Figure 2 :
Figure 2:HLA-B allele cluster from predicted peptide-binding specificity for European American (panel A) and African Americans (panel B).The color key was based on the distance between allele binding specificity.Red color (near 0) indicates high similarity on binding specificity between alleles, and white (near 1) indicates completely different binding specificity.Allele clusters tested for association with TMP-SMX DILI are marked by rectangles.

Table 2 :
HLA alleles associated with TMP-SMX DILI by ethnic group False discovery rate by Benjamini-Hochberg method to correct multiple testing of 85 alleles.

Table 4 :
Amino acids (AAs) showing significant association with TMP-SMX DILI by ethnic group

Exact test Firth logistic regression * Race/Reference HLA allele AA position (residue) B* alleles shared the same residue Case AA Freq Control AA Freq P OR (95% CI) Adjusted P
Hepatology.Author manuscript; available in PMC 2022 January 01.