HLA-B*57 Allele Is Associated with Concomitant Anti-tuberculosis and Antiretroviral Drugs Induced Liver Toxicity in Ethiopians

Drug-induced liver injury (DILI) is a known adverse effect of both anti-tuberculosis (anti-TB) and antiretroviral (ARV) drugs. Recent studies highlight the implications of genetic predispositions to DILI. We performed a case-control study to identify Human Leukocyte Antigen-B (HLA-B) variant alleles associated with anti-TB and ARV co-treatment induced liver toxicity in Ethiopian TB and HIV co-infected patients. A total of 495 newly diagnosed TB and HIV co-infected patients were enrolled and received rifampicin based anti-TB and efavirenz based ARV therapy. Change in liver enzyme level from baseline was monitored 1st, 2nd, 4th, 8th, 12th, and 24th weeks after treatment initiation to identify patients who developed DILI (cases) and those who did not (treatment tolerants). Genomic DNA from 46 cases and 46 sex and age matched treatment tolerants were genotyped for HLA-B variant alleles using Olerup SSP^®HLA-B DNA Typing Kits. The proportion of HLA-B*57 allele carriers in DILI cases (37.0%), particularly in those who developed cholestatic type of DILI (44.8%) was significantly higher compared with those who tolerated the treatment (2.2%). The HLA-B*57 allele frequency was significantly higher in cases (25%) than treatment tolerants (1.1%). In a multivariate logistic analysis, the proportion of patients carrying HLA-B*57 (P = 0.002) and HLA-B*14 (P = 0.014) alleles were significantly higher in DILI cases compared with treatment tolerants. HLA-B*57 was significantly associated with cholestatic (P = 0.001) and mixed (P = 0.017) types of liver toxicity, and mild-to-moderate severity (P = 0.001). Of all HLA-B*57 alleles detected, HLA-B*57:03 accounted 58.3% and HLA-B*57:02 accounted 41.7%. HLA-B*57:01 was not detected. The variant allele frequencies of HLA-B*57:03 (15.2 vs. 0%) and HLA-B*57:02 (9.8 vs. 1.1%) were significantly higher in the DILI cases than treatment tolerants (P < 0.03). We conclude that HLA-B*57 alleles (B*57:03 and B*57:02) confer susceptibility to the development of anti-TB and ARV drugs co-treatment induced liver toxicity, which is mainly of cholestatic type. The possible association of HLA-B*14 with anti-TB and ARV drugs co-treatment induced liver toxicity requires further investigations.

Introduction

Tuberculosis (TB) is the most common opportunistic infection associated with human immunodeficiency virus (HIV) infection, and co-treatment of the two diseases is recommended (Harries et al., 2009). However, simultaneous treatment of TB and HIV infections is challenging due to drug interactions and overlapping toxicities (Cohen and Meintjes, 2010). Antiretroviral (ARV) and anti-tuberculosis (anti-TB) Drugs-induced liver injury (DILI) is a common adverse event, which can be fatal if therapy is not interrupted or changed on time (Devarbhavi et al., 2013; Naidoo et al., 2015; Shamanna et al., 2016). TB-HIV co-infected patients on anti-TB and ARV co-therapy are at a higher risk of developing DILI than TB or HIV only infected patients receiving monotherapy (Yimer et al., 2011, 2014; Mugusi et al., 2012). A recent study in TB/HIV patients on anti-TB and antiretroviral therapy (ART) with high levels of immune activation demonstrated impaired isoniazid clearance, implicating the need for exploring immune response and the risk of DILI (Vinnard et al., 2016). Up to 32% of HIV patients on ART discontinue their treatment or switch therapy mainly due to DILI (Bica et al., 2001), and genetic predisposition contributes partly (Lubomirov et al., 2011). Treatment interruption may increase the risk for developing of multidrug-resistant TB (MDR-TB) and HIV/AIDS (Hirpa et al., 2013). Thus, identifying genetic markers for drug-induced liver toxicity is valuable to identify high-risk patients and to introduce appropriate measures.

Both HIV and TB remain a major problem and co-infection is common in most Sub-Saharan African (SSA) countries including Ethiopia, the second most densely populated country in Africa with an estimated population size of 100 million. Ethiopia is listed among the top 20 high-TB burden countries globally, and one of the high MDR-TB burden countries (Biadglegne et al., 2014; World-Health-Organization, 2016). The rate of new HIV infection in Ethiopia is declining with the estimated number of people living with HIV being 769, 600 in 2014 (World-Health-Organization¹). The scale of ART is increasing in the country as part achieving the UNAID/WHO “90-90-90” target: to diagnose 90% of all HIV positive people, provide ART for 90% of those diagnosed and achieve viral suppression for 90% of those treated, by 2020 (UNAIDS, 2014). ART and anti-TB drug-induced liver toxicity is a common problem in Ethiopia causing treatment discontinuation and hence MDR-TB (Hirpa et al., 2013; Yimer et al., 2014).

Genetic variations in HLA gene is implicated with susceptibility to T-cell mediated adverse events to a wide range of pharmaceuticals making it a candidate gene relevant to pharmacogenetic studies (Barbarino et al., 2015). HLA alleles that are reported to be association with increased risk of idiosyncratic DILI include: HLA-DQB1*02:01 and DQB1*05 to anti-TB drugs (Sharma et al., 2002; Chen et al., 2015), and HLA-B*58:01 and DRB1*01:02 to nevirapine-containing ARV regimens (Phillips et al., 2013). HLA-B*57:01 and A*33:03 variant alleles were also reported as genetic markers for idiosyncratic liver injury induced by flucloxacillin (Daly et al., 2009) and ticlopidine (Hirata et al., 2008) respectively. Susceptibility to amoxicillin-clavulanate-induced liver injury is influenced by multiple HLA class I and II alleles (Lucena et al., 2011). A large genome-wide association study found a strong association of amoxicillin-clavulanate induced liver injury with HLA-A*02:01, HLA-DQB1*06:02, and DRB1*15:01 variant alleles (Lucena et al., 2011). Genetic screening for HLA-B*57:01 and subsequent treatment modifications have been shown to reduce incidence of life-threatening hypersensitivity to abacavir in HIV/AIDS patients carrying the allele (Hughes et al., 2008; Mallal et al., 2008) and HLA-B*15:02 to carbamazepine in Southeast Asian carriers (Amstutz et al., 2014).

Most previous reports investigating genetic risk factors for anti-TB and ARV drugs-induced liver toxicity focused on drug metabolizing enzymes and transporter proteins (Lee et al., 2010; Yimer et al., 2011, 2012). Previously we reported the association of high efavirenz plasma concentration and CYP2B6*6 genotype with DILI in TB-HIV patients (Yimer et al., 2011, 2012; Mugusi et al., 2012). However, only a few studies have explored the association of HLA genes with anti-TB or ARV drugs-induced liver toxicity. Therefore, in this study, we aimed to investigate the possible associations between HLA-B alleles, and anti-TB and ARV drugs co-treatment induced liver injury in TB and HIV co-infected patients in Ethiopia.

Methods

Study Design and Participants

Using a case-control comparative study design, we analyzed data from newly diagnosed TB and HIV co-infected patients, who were enrolled and followed up prospectively to identify the incidence, the pattern, and severity of anti-TB and ARV drugs-induced liver toxicity in Ethiopian patients (Yimer et al., 2014). In brief, 495 TB and HIV co-infected patients with CD4 count ≤ 200 cells/mm³ were recruited from three health institutions: Kazanchis and Beletshachew health centers, and Tikur Anbessa Specialized Hospital in Addis Ababa, Ethiopia, from June 2007 to June 2012. The inclusion criteria were TB and HIV co-infected men and non-pregnant women with age 18 years old and above. Patients were excluded if they had a history of prior treatment for TB/HIV or known pre-existing liver disease.

The study protocol was approved by the Institutional Review Board of College of Health Sciences, Addis Ababa University, the National Research Ethics Review Committee of Ethiopia, and Ethical Review Board of Karolinska Institutet, Sweden. Written informed consent was obtained from all the study participants in accordance with the Declaration of Helsinki.

Drug Treatment

All the study participants received first line ARV drugs containing efavirenz and lamivudine with tenofovir, zidovudine, or stavudine. A short-course anti-TB regimen consisting of rifampicin, isoniazid, pyrazinamide, and ethambutol for the first 2 months followed by rifampicin and isoniazid for the next 4 months was given. The patients did not receive other known hepatotoxic drugs concurrently, except co-trimoxazole prophylaxis that was given for TB and HIV co-infected patients according to the National Treatment Guideline. Change in liver enzymes levels from baseline was monitored on the 1st, 2nd, 4th, 8th, 12th, and 24th weeks after initiation of treatment.

Case Definitions, Severity Grade, and Pattern of Liver Toxicity

The criteria set by the International DILI Expert Working Group were used for DILI case definitions and pattern of liver injury determination (Aithal et al., 2011). The upper limit of normal (ULN) for liver enzymes used for the study population were alanine aminotransferase (ALT 33 U/L for male; 29 U/L for female), aspartate aminotransferase (AST, 41 U/L), alkaline phosphatase (ALP, 128 U/L), and 1.0 mg/dL for total bilirubin (Yimer et al., 2014). All cases recruited met at least one of the following criteria: (1) ALT ≥ 5xULN, (2) ALP ≥ 2xULN, or (3) ALT ≥ 3xULN along with total bilirubin (T Bil) ≥ 2xULN. Treatment tolerants (controls) were individuals who were on anti-TB and ARV drugs co-treatment but did not fulfill the case definitions for DILI and had not presented with clinical signs and symptoms consistent with DILI in the follow-up period.

The pattern of liver toxicity was defined using R-value, where R = (ALT/ULN)/(ALP/ULN). Cases were categorized as hepatocellular (R ≥ 5), cholestatic (R ≤ 2), or mixed (2 < R < 5) pattern of DILI. Clinical severity grading was determined by employing the highest measured values for each of the biochemical parameters (Yimer et al., 2014). Patients with grades one and two severities were grouped together into a “mild-to-moderate” group and those with grades three and four into a “severe” group. Causality assessment for DILI was performed using Roussel Uclaf Causality Assessment Method (RUCAM; Danan and Benichou, 1993).

Among the 495 TB and HIV co-infected patients involved in the initial cohort, 120 experienced DILI in the follow-up period (Yimer et al., 2014). Of these, 80 cases and 275 treatment tolerants had adequate DNA available for further analysis. After excluding patients that had abnormal liver biochemistry prior to starting treatment, or patients who had serological test positive for either hepatitis B virus surface antigen or anti-hepatitis C virus antibody, 46 cases and 46 treatment tolerants that have complete clinical data and matched with respect to gender and age in a 1:1 ratio were used for the current study.

HLA-B Genotyping

Genomic DNA was isolated from peripheral blood using QIAamp DNA Maxi Kit (QIAGEN GmbH. Hilden. Germany). We first screened for HLA-B variant alleles using a low-resolution (two digits) genotyping. HLA-B genotyping was performed using low-resolution Olerup SSP^®HLA-B Typing Kit (Olerup SSP AB, Franzengatan 5, SE-112 51 Stockholm, Sweden). Allele-specific polymerase chain reaction (PCR), using sequence-specific primers was done according to the protocol and recommendations of the manufacturer. The amplified PCR products were analyzed using 2% agarose gel, and the HLA-B allele types were determined using HELMBERG-SCORE software. Low-resolution typing results were recorded with the 2-digit code to ensure a uniform level of HLA resolution for the alleles.

As a next step, high resolution (four digits) typing were done for HLA-B variant alleles that showed significant association with DILI based on the low-resolution genotyping data. HLA-B*57 exhibited a significant association with DILI, and high-resolution subtyping was performed for all HLA-B*57 allele carriers using Olerup SSP^®HLA-B*57:01 Typing Kit (Olerup SSP AB, Franzengatan5, SE-112 51 Stockholm, Sweden) according to the protocol and recommendations of the manufacturer.

Statistical Analysis

Continuous variables were presented by a mean and standard deviation, and categorical variables as numbers and percentages. Univariate logistic regression analysis was used to identify potential independent risk factors for anti-TB and ARV drugs co-treatment induced liver toxicity. Variables with P < 0.1 in the univariate analysis were included in a multivariate logistic analysis. The strength of the associations was estimated by calculating the odds ratio (OR) and 95% confidence interval (CI). Fisher's exact test was used for HLA-B alleles with <5 expected cell count in a 2 × 2 table. To reduce bias in estimating the OR, Haldane's modification was employed (Haldane, 1956) i.e., whenever a zero-count cell was encountered, 0.5 was added to all cells in the 2 × 2 table. P < 0.05 were considered statistically significant. The corrected P-values (Pc) were adjusted by using Bonferroni's correction for multiple comparisons (18 tests) to account for the number of HLA-B alleles observed in the study participants. The statistical analysis was performed using IBM SPSS Statistics for Windows, Version 20.0 (IBM Corp., Armonk, NY, USA).

Results

A total of 92 TB and HIV co-infected patients on anti-TB and ARV drugs co-treatment were involved in this study; 46 treatment induced liver toxicity cases and 46 sex and age matched treatment tolerants. The demographics and clinical characteristics of the study participants are described in Table 1.

TABLE 1

Table 1. Demographics and clinical characteristics of the study participants.

In the univariate analysis, there were statistically significant differences in the CD4 count and Karnofsky score between DILI cases and treatment tolerants (P < 0.05). In a multivariate logistic analysis, baseline CD4 count remained as significant predictors of anti-TB and ARV drugs co-treatment induced liver injury. There were no statistically significant differences in the baseline liver enzyme levels and type of ARV regimens used between the DILI cases and treatment tolerants. More than half of the DILI cases developed the cholestatic type of liver toxicity, and 85% of the cases had mild-to-moderate severity of liver toxicity. All of the cases had a minimum score of three (“possible”) in RUCAM scoring system for DILI.

HLA-B genotype result from the low resolution typing for each study participant is presented in Supplementary Table 1. Comparison of HLA-B allele carriers' proportions between patients who developed DILI vs. treatment tolerants is presented in Table 2. A total of 18 HLA-B variant alleles were detected (Table 2). In the univariate analysis, the proportion of HLA-B*57 and HLA-B*14 alleles carriers were significantly higher in DILI cases than treatment tolerants. Association of being a carrier of HLA-B*57 with increased risk for DILI remained significant after correction for multiple testing. On the other hand, HLA-B*41 was negatively associated with DILI. The multivariate logistic analysis retained HLA-B*57 and HLA-B*14 as significant predictors of concomitant anti-TB and ARV drugs induced liver injury. For HLA-B*57, the association maintained after correcting for multiple comparisons (Pc = 0.036).

TABLE 2

Table 2. Comparison of proportion of HLA-B allele carriers between patients who developed ARV and anti-tuberculosis drugs induced liver injury (cases) and patients who did not (treatment tolerants).

Associations between HLA-B*57, HLA-B*14, and HLA-B*41 with the pattern and severity of liver toxicity was evaluated by comparing the proportion of allele carriers between DILI cases and treatment tolerants (Table 3). Compared with treatment tolerants, the proportion of HLA-B*57 allele carriers were significantly higher in cholestatic and mixed types of DILI. In addition, the proportion of HLA-B*14 allele carriers were higher in the mixed type of DILI cases than the treatment tolerants. There were no statistically significant differences in the hepatocellular type of DILI. Compared with treatment tolerants, the proportion of patients carrying HLA-B*57 allele was significantly higher in the mild-to-moderate DILI group, and those carrying HLA-B*14 in the severe DILI group were over-represented. There was a statistically significant difference in the proportion of HLA-B*41 allele carriers in the mild-to-moderate DILI group compared with treatment tolerants.

TABLE 3

Table 3. Association of HLA-B alleles with the pattern and severity of drug-induced liver injury.

HLA-B variant alleles that showed significant association with DILI from the low-resolution genotyping were further subjected to high resolution (four digit) testing to identify the sub-variant allele. The result from high resolution genotyping was interpreted in conjunction with the results from the prior low-resolution typing. HLA-B*57 was the only variant allele that exhibited a significant association with DILI after correcting for multiple testing. Thus, high-resolution typing was done for all subjects who were genotyped as carriers of HLA-B*57 in the low-resolution typing (see Table 4).

TABLE 4

Table 4. High resolution genotyping data for all HLA-B^*57 allele carriers stratified by DILI type.

Of all HLA-B*57 alleles identified, HLA-B*57:03 accounted 58.3% and HLA-B*57:02 accounted 41.7%. HLA-B*57:01 was not detected. Of the HLA-B*57 allele carrier cases, 10 (58.8%) had HLA-B*57:03 (four homozygous) and the rest seven (41.2%) had HLA-B*57:02 (two homozygous). There was only one HLA-B*57 allele carrier (heterozygous for HLA-B*57:02) out of the 46 treatment-tolerants. The proportion of HLA-B*57:03 and B*57:02 allele carriers were significantly higher in the DILI cases than treatment tolerants [P = 0.01, OR = 26.8 (1.5–47.2) and P = 0.03, OR = 8.1 (1.0–68.6)], respectively.

Comparisons of HLA-B*57, *14, and *41 (from the low-resolution genotyping) as well as HLA-B*57:02 and B*57:03 (from the high-resolution genotyping) allele frequencies between DILI cases and treatment tolerants is presented in Table 5. The overall allele frequency of HLA-B*57 was higher (13.0%) than that of the other HLA-B alleles (7.6% for B*14 and 7.1% for B*41). The allele frequency of HLA-B*57 was higher in DILI cases (25.0%) compared to treatment tolerants (1.1%). The HLA-B*57:03 and B*57:02 allele frequencies in DILI cases (15.2 and 9.8%) were higher than the treatment tolerants (0 and 1.1%), respectively.

TABLE 5

Table 5. Comparison of HLA-B allele frequencies distribution between patients who developed ARV and anti-tuberculosis drugs induced liver toxicity (cases) and patients who did not (treatment tolerants).

Discussion

In the present study, we investigated the association of HLA-B variant alleles with risk for concomitant anti-TB and ARV drugs induced liver toxicity. The proportion of HLA-B*57 allele carriers in Ethiopian patients who developed anti-TB and ARV drugs induced liver toxicity (37.0%), particularly in those who developed cholestatic type of liver toxicity (44.8%) was significantly higher compared with those who tolerated the treatment (2.2%). The proportion of HLA-B*14 allele carriers who developed DILI (21.7%) was also significantly higher compared with the treatment tolerants (6.5%). These indicate that HLA-B*57 and B*14 allele carriers might be at a higher risk of developing anti-TB and ARV drugs co-treatment induced liver toxicity. Accordingly, these variant alleles might play important roles in the pathogenesis of immune-mediated liver toxicity during anti-TB and ARV drugs co-treatment. The HLA-B*57 and B*14 molecules may function as endogenous antigen presenting molecules for the drugs/metabolites to HLA-restricted cytotoxic T-cell activation (Pichler, 2002). To our knowledge, this is the first report to investigate the association of HLA-B*57, B*14, and B*41 alleles with anti-TB and ARV drugs co-treatment induced liver toxicity.

The HLA-B*57 variant allele, which was observed in significantly higher proportion among DILI cases than treatment tolerants, had a high specificity (97.8%) and positive predictive value (94.4%). Hence, HLA-B*57 is likely to be an important predictor for anti-TB and ARV drugs co-treatment induced liver injury. However, there could be additional yet unidentified genetic markers and non-genetic risk factors involved in the pathogenesis of DILI. The matched case-control design used in this study minimizes effects of potential confounders and may increase power to identify genetic associations. Although this limits us from exploring associations of the matching variables such as sex, body mass index, Karnofsky score, CD4 count, and HIV viral load which were independently and significantly associated with the risk of developing DILI (Yimer et al., 2011, 2014). Association of CD4 cell counts and Karnofsky score as risk factors for DILI were also found in this study, although the others were not significant.

Studies suggest that a particular HLA-B allele may exert a protective effect against certain adverse drug reactions as evidenced by lower allele carrier frequencies in cases compared with treatment tolerants. HLA-DQA1*01:02 was identified as a protective variant for anti-TB drugs induced hepatotoxicity (Sharma et al., 2002). HLA-B*40:01 and HLA-B*07:02 were also identified as protective variant alleles for carbamazepine-induced severe cutaneous adverse reactions (Alfirevic et al., 2006; Hung et al., 2006). In our study, statistically significant lower allele carrier rate of HLA-B*41 was noted in the DILI cases compared with the treatment tolerants (6.5 vs. 21.7%), but this effect did not reach statistical significance after correcting for multiple comparisons. Further, analysis is required to clarify the role of HLA-B*41 in the prevention of development of anti-TB and ARV drugs co-treatment induced liver injury.

Our result indicated a positive association of HLA-B*57 allele with mild-to-moderate liver injury and the HLA-B*14 allele with severe liver injury. On the other hand, HLA-B*41 allele was negatively associated with mild-to-moderate liver injury. Accordingly, the association of HLA-B alleles with anti-TB and ARV drugs co-treatment induced liver injury may seem to depend on the severity of liver injury. The HLA-B*57 allele may be critical for the initiation of the immune response to cause DILI and the HLA-B*14 allele for the progression to severe degree of liver injury. On the other hand, the HLA-B*41 allele seems to play a role in the prevention of development of mild-to-moderate liver injury due to anti-TB and ARV drugs co-treatment. These findings warrant further investigation in a larger DILI case samples for each severity grade of liver injury.

DILI can be hepatocellular (predominant rise in ALT), cholestatic (predominant rise in ALP), or mixed type liver injury (Hussaini and Farrington, 2014). Recently, we conducted a prospective observational study to evaluate the incidence, type, severity, and predisposing risk factors for DILI in a large well-defined TB and/or HIV patient cohort receiving either anti-TB drugs alone, ARV drugs alone or concomitant anti-TB and ARV therapy (Yimer et al., 2014). We found rates of hepatocellular DILI being highest among patients treated with anti-TB drugs alone than patients treated with ARV drugs alone or co-treated with anti-TB drugs. On the other hand, the rates of cholestatic DILI was highest among patients treated with efavirenz based-ARV drugs alone than patients treated with anti-TB drugs alone or with ARV drugs (Yimer et al., 2014). DILI due to anti-TB drugs in TB mono-infected patients is known to be more of hepatocellular type. In the present study, most of the TB-HIV co-infected patients treated with concomitant anti-TB and efavirenz based-ARV drugs developed cholestatic DILI cases. Apparently, there is a significant contribution from efavirenz based-ARV drugs toward developing cholestatic type of DILI. Indeed, the significant association of HLA-B*57 variant allele with cholestatic type DILI identified in the present study might reflect for ARV-drugs induced hepatotoxicity. However, the role of HLA variant alleles for predisposition to anti-TB DILI cannot be ruled out (Sharma et al., 2002; Chen et al., 2015). Further, studies are necessary to investigate the association of HLA allele carrier status with anti-TB drugs alone as well as ARV drugs alone-induced liver injury.

Major histocompatibility complex (MHC) class I and class II-mediated immunological reactions are implicated in DILI, particularly in the cholestatic type that involves damage to the biliary system (Andrade et al., 2004; Daly, 2010). In line with this, carrier status of HLA-B*57 was significantly higher in patients who presented with the cholestatic type of DILI (44.8%) and mixed type (26.7%) compared with those who tolerated the treatment (2.2%). None of the patients who developed hepatocellular DILI were carriers of HLA-B*57 variant allele.

HLA-B*57 allele is associated with long-term non-progressive chronic HIV-1 infection by restricting cytotoxic T-lymphocyte response (Goulder et al., 2000). HLA-B*57:01 and B*57:03 are the most prevalent HLA-B*57 subtypes in Caucasian and African populations, respectively (Pelak et al., 2010; Apps et al., 2013). The HLA-B*57:01 and B*57:03 alleles are protective against HIV disease progression, and appear to present identical Gag epitopes (Payne et al., 2014). HLA-B*57:01, B*57:02, and B*57:03 share more than 90% sequence homology and as such have peptide-binding repertoires which substantially overlap (Illing et al., 2012; Ogese et al., 2017). Although in ART naïve patients HLA-B*57 (B*57:01 in Europe and US, B*57:03 in black Africans) confers protective effect against HIV-1 disease progression to AIDS (Costello et al., 1999; Migueles et al., 2000; López-Larrea et al., 2005; Frater et al., 2007), it may exert contradictory effect on treatment outcome when the disease course is altered by ARV therapy (Dold et al., 2015). Previous studies reported the association of HLA-B*57 with increased all causes of mortality (Dold et al., 2015) and reduced virological responses during ARV therapy (Kuniholm et al., 2011). HLA-B*57 allele is also known to be associated with immune-mediated drug-induced hypersensitivity reactions. Carriers of HLA-B*57:01 allele are at higher risk of developing abacavir-induced hypersensitivity reactions (Hetherington et al., 2002), whereas HLA-B*57:03 is associated with spondylarthropathies (López-Larrea et al., 2005). Indeed, genetic screening for HLA-B*57:01 variant allele has been shown to reduce drug toxicity and subsequently led to a labeled recommendation of routine screening before treatment initiation (Hughes et al., 2008; Mallal et al., 2008).

The frequency and subtypes of HLA-B*57 variant alleles display wide inter-ethnic variability globally ranging from 0 to 22.5% (http://www.allelefrequencies.net/). The overall frequency of HLA-B*57 in our study population from Ethiopia is 13% which is relatively high. HLA-B*57:01 occurs in Asians and Caucasians (up to 5%). The HLA-B*57:03 and B*57:02 variant alleles commonly occur in black population reaching up to 3 and 7% allele frequencies, respectively. Interestingly HLA-B*57:01 was not detected, and it may be rare or absent in Ethiopians similar to other black Africans where the allele frequency is <1%. The overall HLA-B*57:03 and B*57:02 allele frequencies in our TB/HIV co-infected study population was 7.6 and 5.4%, respectively, although the frequencies in healthy Ethiopians is yet unknown. Interestingly, HLA-B*57:03 and B*57:02 allele frequencies in DILI cases (15.2 and 9.8%) were significantly higher than the allele frequencies in the treatment-tolerants (0 and 1.1%), respectively.

There were some limitations in this study. First, as DILI is a rare event, it was not easy to get large number of cases (four years were required to collect the DILI cases in this study), which subsequently resulted in a small number of samples for sub-group analysis. The second limitation is that as drug combinations are the current treatment protocols for TB and HIV co-infections, and hence we cannot link the risk variant allele to a specific drug or class of drug(s). Since first line anti-TB and HIV treatment regimen consists of combination therapy, it is not possible to study individual drug-induced liver toxicity in TB and HIV co-infected patients for ethical reasons. However, our study represents an important first step in applying HLA-B typing to identify genetic variants for anti-TB and ARV drugs co-treatment induced liver injury.

Identification of genetic risk factors for anti-TB, and ARV drugs co-treatment induced liver injury is essential for patient safety. The HLA risk alleles predisposing to immune-mediated anti-TB and ARV drugs induced liver toxicity in black African population are not well investigated. A common problem encountered in HLA genotyping is inability to determine the variant alleles accurately using a simple genotyping procedure. This is mainly due to the extensive genetic diversity in HLA gene locus. Accurate allele-level HLA typing using the current methods requires high workload, cost and time. Because of extreme genetic variation of the HLA locus, pharmacogenetic testing for routine clinical practice is increasingly challenging in resource limited settings. However, the recent development of second-generation sequencing methods provides the possibility of sequencing a single DNA strand in isolation. Establishing a straightforward and affordable genotyping method for accurate HLA typing to identify patients at risk of developing drug-induced adverse events may lay the ground for the future application of pharmacogenetic testing in clinical practice for globalized personalized medicine.

In conclusion, HLA-B variant alleles may play important roles in determining the risk and severity of concomitant anti-TB and ARV drugs induced liver toxicity. HLA-B*57 variant alleles (HLA-B*57:03 and HLA-B*57:02) are risk factors to develop anti-TB and ARV drugs co-treatment induced liver injury, mainly of cholestatic type and mild DILI cases. The possible risk association of HLA-B*14 allele with severe DILI and the protective association of HLA-B*41 require further investigations. Additional studies are necessary to understand the roles of the identified HLA-B variant alleles in the pathogenesis of anti-TB and ARV drugs co-treatment induced liver toxicity.

Author Contributions

EA, JK, and EM conceived and designed the study; EA, GY, AH, EM collected the data; EA, JK, and ZP performed the experiment and analyzed the data; EA and ZP wrote this paper. All authors revised/edited the manuscript and approved for submission.

Funding

The study was financially supported by grants from European and Developing Countries Clinical Trials Partnership (NL) (CG_TA.05.40204_005) and from Swedish Research Council (Vetenskapsrådet, Grant number: 2015-03295). This work was also supported partly by the NIH/Fogarty International Center Global Infectious Diseases grant D43TW009127.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to express our heartfelt gratitude to all the study participants.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fphar.2017.00090/full#supplementary-material

Footnotes

1. ^World Health Organization. Ethiopia. HIV/AIDS. Available online at: http://www.afro.who.int/en/ethiopia/country-programmes/topics/4480-hivaids.html (Accessed January 29, 2017).

References

Aithal, G. P., Watkins, P. B., Andrade, R. J., Larrey, D., Molokhia, M., Takikawa, H., et al. (2011). Case definition and phenotype standardization in drug-induced liver injury. Clin. Pharmacol. Ther. 89, 806–815. doi: 10.1038/clpt.2011.58

PubMed Abstract | CrossRef Full Text | Google Scholar

Alfirevic, A., Jorgensen, A. L., Williamson, P. R., Chadwick, D. W., Park, B. K., and Pirmohamed, M. (2006). HLA-B locus in Caucasian patients with carbamazepine hypersensitivity. Pharmacogenomics 7, 813–818. doi: 10.2217/14622416.7.6.813

PubMed Abstract | CrossRef Full Text | Google Scholar

Amstutz, U., Shear, N. H., Rieder, M. J., Hwang, S., Fung, V., Nakamura, H., et al. (2014). Recommendations for HLA-B*15:02 and HLA-A*31:01 genetic testing to reduce the risk of carbamazepine-induced hypersensitivity reactions. Epilepsia 55, 496–506. doi: 10.1111/epi.12564

PubMed Abstract | CrossRef Full Text | Google Scholar

Andrade, R. J., Lucena, M. I., Alonso, A., García-Cortes, M., García-Ruiz, E., Benitez, R., et al. (2004). HLA class II genotype influences the type of liver injury in drug-induced idiosyncratic liver disease. Hepatology 39, 1603–1612. doi: 10.1002/hep.20215

PubMed Abstract | CrossRef Full Text | Google Scholar

Apps, R., Qi, Y., Carlson, J. M., Chen, H., Gao, X., Thomas, R., et al. (2013). Influence of HLA-C expression level on HIV control. Science 340, 87–91. doi: 10.1126/science.1232685

PubMed Abstract | CrossRef Full Text | Google Scholar

Barbarino, J. M., Kroetz, D. L., Klein, T. E., and Altman, R. B. (2015). PharmGKB summary: very important pharmacogene information for human leukocyte antigen B. Pharmacogenet. Genomics 25, 205–221. doi: 10.1097/FPC.0000000000000118

PubMed Abstract | CrossRef Full Text | Google Scholar

Biadglegne, F., Sack, U., and Rodloff, A. C. (2014). Multidrug-resistant tuberculosis in Ethiopia: efforts to expand diagnostic services, treatment and care. Antimicrob. Resist. Infect. Control 3:31. doi: 10.1186/2047-2994-3-31

PubMed Abstract | CrossRef Full Text | Google Scholar

Bica, I., McGovern, B., Dhar, R., Stone, D., McGowan, K., Scheib, R., et al. (2001). Increasing mortality due to end-stage liver disease in patients with human immunodeficiency virus infection. Clin. Infect. Dis. 32, 492–497. doi: 10.1086/318501

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, R., Zhang, Y., Tang, S., Lv, X., Wu, S., Sun, F., et al. (2015). The association between HLA-DQB1 polymorphism and antituberculosis drug-induced liver injury: a Case-Control Study. J. Clin. Pharm. Ther. 40, 110–115. doi: 10.1111/jcpt.12211

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, K., and Meintjes, G. (2010). Management of individuals requiring antiretroviral therapy and TB treatment. Curr. Opin. HIV AIDS 5, 61–69. doi: 10.1097/COH.0b013e3283339309

PubMed Abstract | CrossRef Full Text

Costello, C., Tang, J., Rivers, C., Karita, E., Meizen-Derr, J., Allen, S., et al. (1999). HLA-B*5703 independently associated with slower HIV-1 disease progression in Rwandan women. AIDS 13, 1990–1991. doi: 10.1097/00002030-199910010-00031

PubMed Abstract | CrossRef Full Text | Google Scholar

Daly, A. K. (2010). Drug-induced liver injury: past, present and future. Pharmacogenomics 11, 607–611. doi: 10.2217/pgs.10.24

PubMed Abstract | CrossRef Full Text | Google Scholar

Daly, A. K., Donaldson, P. T., Bhatnagar, P., Shen, Y., Pe'er, I., Floratos, A., et al. (2009). HLA-B*5701 genotype is a major determinant of drug-induced liver injury due to flucloxacillin. Nat. Genet. 41, 816–819. doi: 10.1038/ng.379

PubMed Abstract | CrossRef Full Text | Google Scholar

Danan, G., and Benichou, C. (1993). Causality assessment of adverse reactions to drugs–I. A novel method based on the conclusions of international consensus meetings: application to drug-induced liver injuries. J. Clin. Epidemiol. 46, 1323–1330. doi: 10.1016/0895-4356(93)90101-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Devarbhavi, H., Singh, R., Patil, M., Sheth, K., Adarsh, C. K., and Balaraju, G. (2013). Outcome and determinants of mortality in 269 patients with combination anti-tuberculosis drug-induced liver injury. J. Gastroenterol. Hepatol. 28, 161–167. doi: 10.1111/j.1440-1746.2012.07279.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dold, L., Ahlenstiel, G., Althausen, E., Luda, C., Schwarze-Zander, C., Boesecke, C., et al. (2015). Survival and HLA-B*57 in HIV/HCV co-infected patients on highly active antiretroviral therapy (HAART). PLoS ONE 10:e0134158. doi: 10.1371/journal.pone.0134158

PubMed Abstract | CrossRef Full Text | Google Scholar

Frater, A. J., Brown, H., Oxenius, A., Günthard, H. F., Hirschel, B., Robinson, N., et al. (2007). Effective T-cell responses select human immunodeficiency virus mutants and slow disease progression. J. Virol. 81, 6742–6751. doi: 10.1128/JVI.00022-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Goulder, P. J., Tang, Y., Pelton, S. I., and Walker, B. D. (2000). HLA-B57-restricted cytotoxic T-lymphocyte activity in a single infected subject toward two optimal epitopes, one of which is entirely contained within the other. J. Virol. 74, 5291–5299. doi: 10.1128/JVI.74.11.5291-5299.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Haldane, J. B. (1956). The estimation and significance of the logarithm of a ratio of frequencies. Ann. Hum. Genet. 20, 309–311. doi: 10.1111/j.1469-1809.1955.tb01285.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Harries, A. D., Zachariah, R., and Lawn, S. D. (2009). Providing HIV care for co-infected tuberculosis patients: a perspective from sub-Saharan Africa. Int. J. Tuberc. Lung Dis. 13, 6–16.

PubMed Abstract | Google Scholar

Hetherington, S., Hughes, A. R., Mosteller, M., Shortino, D., Baker, K. L., Spreen, W., et al. (2002). Genetic variations in HLA-B region and hypersensitivity reactions to abacavir. Lancet 359, 1121–1122. doi: 10.1016/S0140-6736(02)08158-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirata, K., Takagi, H., Yamamoto, M., Matsumoto, T., Nishiya, T., Mori, K., et al. (2008). Ticlopidine-induced hepatotoxicity is associated with specific human leukocyte antigen genomic subtypes in Japanese patients: a preliminary case-control study. Pharmacogenomics J. 8, 29–33. doi: 10.1038/sj.tpj.6500442

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirpa, S., Medhin, G., Girma, B., Melese, M., Mekonen, A., Suarez, P., et al. (2013). Determinants of multidrug-resistant tuberculosis in patients who underwent first-line treatment in Addis Ababa: a case control study. BMC Public Health 13:782. doi: 10.1186/1471-2458-13-782

PubMed Abstract | CrossRef Full Text | Google Scholar

Hughes, S., Hughes, A., Brothers, C., Spreen, W., Thorborn, D., and Team, C. N. A. S. (2008). PREDICT-1 (CNA106030): the first powered, prospective trial of pharmacogenetic screening to reduce drug adverse events. Pharm. Stat. 7, 121–129. doi: 10.1002/pst.286

PubMed Abstract | CrossRef Full Text | Google Scholar

Hung, S. I., Chung, W. H., Jee, S. H., Chen, W. C., Chang, Y. T., Lee, W. R., et al. (2006). Genetic susceptibility to carbamazepine-induced cutaneous adverse drug reactions. Pharmacogenet. Genomics 16, 297–306. doi: 10.1097/01.fpc.0000199500.46842.4a

PubMed Abstract | CrossRef Full Text | Google Scholar

Hussaini, S. H., and Farrington, E. A. (2014). Idiosyncratic drug-induced liver injury: an update on the 2007 overview. Expert Opin. Drug Saf. 13, 67–81. doi: 10.1517/14740338.2013.828032

PubMed Abstract | CrossRef Full Text | Google Scholar

Illing, P. T., Vivian, J. P., Dudek, N. L., Kostenko, L., Chen, Z., Bharadwaj, M., et al. (2012). Immune self-reactivity triggered by drug-modified HLA-peptide repertoire. Nature 486, 554–558. doi: 10.1038/nature11147

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuniholm, M. H., Gao, X., Xue, X., Kovacs, A., Anastos, K., Marti, D., et al. (2011). Human leukocyte antigen genotype and risk of HIV disease progression before and after initiation of antiretroviral therapy. J. Virol. 85, 10826–10833. doi: 10.1128/JVI.00804-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, S. W., Chung, L. S., Huang, H. H., Chuang, T. Y., Liou, Y. H., and Wu, L. S. (2010). NAT2 and CYP2E1 polymorphisms and susceptibility to first-line anti-tuberculosis drug-induced hepatitis. Int. J. Tuberc. Lung Dis. 14, 622–626.

PubMed Abstract | Google Scholar

López-Larrea, C., Njobvu, P. D., González, S., Blanco-Gelaz, M. A., Martínez-Borra, J., and López-Vázquez, A. (2005). The HLA-B*5703 allele confers susceptibility to the development of spondylarthropathies in Zambian human immunodeficiency virus-infected patients with slow progression to acquired immunodeficiency syndrome. Arthritis Rheum. 52, 275–279. doi: 10.1002/art.20722

PubMed Abstract | CrossRef Full Text | Google Scholar

Lubomirov, R., Colombo, S., di Iulio, J., Ledergerber, B., Martinez, R., Cavassini, M., et al. (2011). Association of pharmacogenetic markers with premature discontinuation of first-line anti-HIV therapy: an observational cohort study. J. Infect. Dis. 203, 246–257. doi: 10.1093/infdis/jiq043

PubMed Abstract | CrossRef Full Text | Google Scholar

Lucena, M. I., Molokhia, M., Shen, Y., Urban, T. J., Aithal, G. P., and Andrade, R. J., et al (2011). Susceptibility to amoxicillin-clavulanate-induced liver injury is influenced by multiple HLA class I and II alleles. Gastroenterology 141, 338–347. doi: 10.1053/j.gastro.2011.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Mallal, S., Phillips, E., Carosi, G., Molina, J. M., Workman, C., and Tomazic, J., et al (2008). HLA-B*5701 screening for hypersensitivity to abacavir. N.Engl. J. Med. 358, 568–579. doi: 10.1056/NEJMoa0706135

PubMed Abstract | CrossRef Full Text | Google Scholar

Migueles, S. A., Sabbaghian, M. S., Shupert, W. L., Bettinotti, M. P., Marincola, F. M., Martino, L., et al. (2000). HLA B*5701 is highly associated with restriction of virus replication in a subgroup of HIV-infected long term nonprogressors. Proc. Natl. Acad. Sci. U.S.A. 97, 2709–2714. doi: 10.1073/pnas.050567397

PubMed Abstract | CrossRef Full Text | Google Scholar

Mugusi, S., Ngaimisi, E., Janabi, M., Minzi, O., Bakari, M., Riedel, K. D., et al. (2012). Liver enzyme abnormalities and associated risk factors in HIV patients on efavirenz-based HAART with or without tuberculosis co-infection in tanzania. PLoS ONE 7:e40180. doi: 10.1371/journal.pone.0040180

PubMed Abstract | CrossRef Full Text | Google Scholar

Naidoo, S., Evans, D., Jong, E., Mellet, K., and Berhanu, R. (2015). Outcomes of TB/HIV co-infected patients presenting with antituberculosis drug-induced liver injury. S.Afr. Med. J. 105, 393–396. doi: 10.7196/SAMJ.8217

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogese, M. O., Ahmed, S., Alferivic, A., Betts, C. J., Dickinson, A., Faulkner, L., et al. (2017). New approaches to investigate drug-induced hypersensitivity. Chem. Res. Toxicol. 30, 239–259. doi: 10.1021/acs.chemrestox.6b00333

PubMed Abstract | CrossRef Full Text | Google Scholar

Payne, R. P., Branch, S., Kløverpris, H., Matthews, P. C., Koofhethile, C. K., Strong, T., et al. (2014). Differential escape patterns within the dominant HLA-B*57:03-restricted HIV Gag epitope reflect distinct clade-specific functional constraints. J. Virol. 88, 4668–4678. doi: 10.1128/JVI.03303-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Pelak, K., Goldstein, D. B., Walley, N. M., Fellay, J., Ge, D., Shianna, K. V., et al. (2010). Host determinants of HIV-1 control in African Americans. J. Infect. Dis. 201, 1141–1149. doi: 10.1086/651382

PubMed Abstract | CrossRef Full Text | Google Scholar

Phillips, E., Bartlett, J. A., Sanne, I., Lederman, M. M., Hinkle, J., Rousseau, F., et al. (2013). Associations between HLA-DRB1*0102, HLA-B*5801, and hepatotoxicity during initiation of nevirapine-containing regimens in South Africa. J. Acquir. Immune Defic. Syndr. 62, e55–e57. doi: 10.1097/QAI.0b013e31827ca50f

PubMed Abstract | CrossRef Full Text | Google Scholar

Pichler, W. J. (2002). Modes of presentation of chemical neoantigens to the immune system. Toxicology 181–182, 49–54. doi: 10.1016/S0300-483X(02)00254-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Shamanna, S. B., Naik, R. R., and Hamide, A. (2016). Causes of liver disease and its outcome in HIV-infected individuals. Indian J. Gastroenterol. 35, 310–314. doi: 10.1007/s12664-016-0676-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, S. K., Balamurugan, A., Saha, P. K., Pandey, R. M., and Mehra, N. K. (2002). Evaluation of clinical and immunogenetic risk factors for the development of hepatotoxicity during antituberculosis treatment. Am. J. Respir. Crit. Care Med. 166, 916–919. doi: 10.1164/rccm.2108091

PubMed Abstract | CrossRef Full Text | Google Scholar

UNAIDS (2014). 90-90-90 An Ambitious Treatment Target to Help End the AIDS Epidemic. Geneva. Available online at: http://www.unaids.org/en/resources/documents/2014/90-90-90 (Accessed January 29, 2017).

Vinnard, C., Ravimohan, S., Tamuhla, N., Ivaturi, V., Pasipanodya, J., Srivastava, S., et al. (2016). Isoniazid clearance is impaired among human immunodeficiency virus/tuberculosis patients with high levels of immune activation. Br. J. Clin. Pharmacol. doi: 10.1111/bcp.13172. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

World-Health-Organization (2016). Global Tuberculosis Report 2016. Geneva. Available online at: http://apps.who.int/iris/bitstream/10665/250441/1/9789241565394-eng.pdf?ua=1 (Accessed January 29, 2017).

Yimer, G., Amogne, W., Habtewold, A., Makonnen, E., Ueda, N., Suda, A., et al. (2012). High plasma efavirenz level and CYP2B6*6 are associated with efavirenz-based HAART-induced liver injury in the treatment of naive HIV patients from Ethiopia: a prospective cohort study. Pharmacogenomics J. 12, 499–506. doi: 10.1038/tpj.2011.34

PubMed Abstract | CrossRef Full Text | Google Scholar

Yimer, G., Gry, M., Amogne, W., Makonnen, E., Habtewold, A., Petros, Z., et al. (2014). Evaluation of patterns of liver toxicity in patients on antiretroviral and anti-tuberculosis drugs: a prospective four arm observational study in ethiopian patients. PLoS ONE 9:e94271. doi: 10.1371/journal.pone.0094271

PubMed Abstract | CrossRef Full Text | Google Scholar

Yimer, G., Ueda, N., Habtewold, A., Amogne, W., Suda, A., Riedel, K. D., et al. (2011). Pharmacogenetic & pharmacokinetic biomarker for efavirenz based ARV and rifampicin based anti-TB drug induced liver injury in TB-HIV infected patients. PLoS ONE 6:e27810. doi: 10.1371/journal.pone.0027810

PubMed Abstract | CrossRef Full Text | Google Scholar