recast-jats-builde082d6219abmjopenBMJ OpenbmjopenBMJ OpenBMJ OpenBMJ Open2044-60552044-6055British Medical Journal Publishing Groupbmjopen-2018-02794010.1136/bmjopen-2018-027940/bmjopen/9/8/e027940.atomRespiratory medicineResearchOpen accessRespiratory medicineSpecial collectionsBMJ OpenRespiratory medicineAssociation of genetic polymorphisms of CYP2E1, NAT2, GST and SLCO1B1 with the risk of anti-tuberculosis drug-induced liver injury: a systematic review and meta-analysisYangSeungwon1HwangSe Jung2ParkJung Yun3ChungEun Kyoung24LeeJangik I35 Department of Pharmacy and Yonsei Institute of Pharmaceutical Science, Yonsei University, Incheon, Republic of Korea Department of Pharmacy, College of Pharmacy, Kyung Hee University, Seoul, Republic of Korea College of Pharmacy, Seoul National University, Seoul, Republic of Korea Department of Pharmacy, Kyung Hee University Hospital at Gangdong, Seoul, Republic of Korea College of Pharmacy and Research Institute of Pharmaceutical Sciences, Seoul National University, Seoul, Republic of Korea Dr Eun Kyoung Chung; cekchung@khu.ac.kr and Professor Jangik I Lee; jangik.lee@snu.ac.kr 8201918201957201957201998e0279401911201811062019© Author(s) (or their employer(s)) 2019. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.2019This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.Objectives

The objective of this study was to investigate the association between genetic polymorphisms of N-acetyltransferase 2 (NAT2), cytochrome P450 2E1 (CYP2E1), glutathione S-transferase (GST) and solute carrier organic anion transporter family member 1B1 (SLCO1B1) and the risk of anti-tuberculosis drug-induced liver injury (ATDILI).

Design

Systematic review and meta-analysis.

Data sources

PubMed, Embase, Web of Science and Cochrane Reviews databases were searched through April 2019.

Eligibility criteria

We included case-control or cohort studies investigating an association between NAT2, CYP2E1, GST or SLCO1B1 polymorphisms and the ATDILI risk in patients with tuberculosis.

Data extraction and synthesis

Three authors screened articles, extracted data and assessed study quality. The strength of association was evaluated for each gene using the pooled OR with a 95% CI based on the fixed-effects or random-effects model. Sensitivity analysis was performed to confirm the reliability and robustness of the results.

Results

Fifty-four studies were included in this analysis (n=26 for CYP2E1, n=35 for NAT2, n=19 for GST, n=4 for SLCO1B1). The risk of ATDILI was significantly increased with the following genotypes: CYP2E1 RsaI/PstI c1/c1 (OR=1.39, 95% CI 1.06 to 1.83), NAT2 slow acetylator (OR=3.30, 95% CI 2.65 to 4.11) and GSTM1 null (OR=1.30, 95% CI 1.12 to 1.52). No significant association with ATDILI was found for the genetic polymorphisms of CYP2E1 DraI, GSTT1, GSTM1/GSTT1, SLCO1B1 388A>G and SLCO1B1 521T>C (p>0.05).

Conclusions

ATDILI is more likely to occur in patients with NAT2 slow acetylator genotype, CYP2E1 RsaI/PstI c1/c1 genotype and GSTM1 null genotype. Close monitoring may be warranted for patients with these genotypes.

anti-tuberculosis drug-induced liver injurygenetic polymorphismsmeta-analysisdrug-metabolizing enzymedrug transportertuberculosisspecial-featureunlockedspecial-propertycontains-inline-supplementary-materialStrengths and limitations of this study

This is the first meta-analysis to evaluate the association between the risk of anti-tuberculosis drug-induced liver injury (ATDILI) and solute carrier organic anion transporter family member 1B1 (SLCO1B1) in patients with tuberculosis.

We included most updated studies with the large sample sizes to better clarify the association of genetic polymorphisms with the risk of ATDILI.

The effect of anti-tuberculosis drug dosages on the risk of ATDILI could not be accounted for in this study due to the lack of drug dosing information in the majority of the included studies.

Introduction

Tuberculosis is a rampant infectious disease caused by Mycobacterium tuberculosis. It poses a major public health threat globally with approximately 1.3 million deaths and 10 million new cases in 2017.1 The mainstay of first-line tuberculosis treatment is a four-drug combination regimen of isoniazid (INH), rifampicin (RIF), pyrazinamide (PZA) and ethambutol (EMB) during the first 2 months, followed by INH and RIF for additional 4 months.2 3 The currently recommended therapy for tuberculosis is highly effective, resulting in high cure rates if patients are adherent to therapy.4 However, treatment adherence is often suboptimal in patients receiving the combination anti-tuberculosis therapy due to many adverse drug reactions, some of which are considered detrimental.5 One of the common adverse drug reactions associated with anti-tuberculosis medications is anti-tuberculosis drug-induced liver injury (ATDILI) affecting 2% to 28% of patients with tuberculosis.6 ATDILI could be potentially serious and fatal, resulting in the treatment interruption and ultimately, treatment failure.7 8

Recently, increasing evidence suggests an association between the risk of ATDILI and genetic polymorphisms of drug-metabolising enzymes (DMEs) and drug transporters.9 10 Altered enzyme activity due to polymorphic genotypes of various DMEs including cytochrome P450 2E1 (CYP2E1), N-acetyltransferase 2 (NAT2) and glutathione S-transferase (GST) can result in the accumulation of toxic substances in the liver, leading to the development of ATDILI.11 However, conflicting results have been reported regarding the association between the risk of ATDILI and genetic polymorphisms of various DMEs in patients with tuberculosis.9 12 13 In addition to DMEs, drug transporters have been emerging as a key determinant of the pharmacokinetics and pharmacodynamics of a drug.14 Among various drug transporters, organic anion transporting polypeptide 1B1 (OATP1B1), encoded by SLCO1B1, is the major influx transporter responsible for hepatic uptake of RIF.15 Although several studies have previously examined the association between SLCO1B1 polymorphisms and the risk of ATDILI, conflicting results have been reported regarding the effect of SLCO1B1 polymorphisms on ATDILI risk. Therefore, an updated meta-analysis has been warranted to confirm the association between the ATDILI risk and genetic polymorphisms of DMEs. In our preliminary literature search, several polymorphic genes, including many DMEs, transporters and other genes such as those involved in the immune system, were identified to have an association with the risk of ATDILI. Among these, sufficient, published information was available to confirm the effect of CYP2E1, NAT2, GST, and SLCO1B1 genetic polymorphisms on the ATDILI risk through meta-analysis.

Objectives

The objective of this meta-analysis was to evaluate the association between the risk of ATDILI and genetic polymorphisms of CYP2E1, NAT2, GST and SLCO1B1 in patients with tuberculosis.

Methods

This study was in compliance with the Meta-analysis Of Observational Studies in Epidemiology checklist for reporting the study design, search strategy, methods, results and conclusions (online supplementary table S1). Three authors (SY, JYP and SJH) independently conducted a literature search, study selection, quality assessment and data extraction. Any discrepancies were adjudicated by corresponding authors (JIL and EKC).

SP110.1136/bmjopen-2018-027940.supp1

Search strategy

Electronic databases of PubMed, Embase, Web of Science and Cochrane Reviews were systematically searched from their inception to April 2019 to identify relevant studies evaluating the association of NAT2, CYP2E1, GST and SLCO1B1 polymorphisms with ATDILI risk. A comprehensive literature search was conducted using a combination of the following keywords and Medical Subject Heading terms: (‘genetic polymorphism’ or ‘NAT2’ or ‘CYP2E1’ or ‘GST’ or ‘SLCO1B1’ or ‘drug-metabolizing enzymes’ or ‘drug transporter’) AND (‘anti-tuberculosis agents drug-induced liver injuries’ or ‘hepatotoxicity’). The detailed search strategies for each electronic database used in this analysis are presented in online supplementary table S2. The reference lists in the selected reviews and meta-analyses were reviewed to ensure the inclusion of all relevant evidence in this analysis.

Study selection

Studies were considered eligible if they met all of the following inclusion criteria: (1) studies in patients with tuberculosis receiving anti-tuberculosis drug regimen; (2) studies with the control group of patients with tuberculosis, tolerant of anti-tuberculosis medications; (3) studies evaluating the association between the occurrence of ATDILI and genetic polymorphisms of CYP2E1, NAT2, GST and SLCO1B1 388A>G and 521T>C and (4) case-control or cohort observational studies. Excluded studies were as follows: (1) studies available only in the form of abstracts or meeting posters, (2) review or meta-analysis articles, (3) studies providing insufficient data necessary for the statistical data analysis, (4) studies in non-English language, (5) non-human studies including animal and in vitro studies, (6) studies with unpublished data, (7) studies providing insufficient information on genotyping methods and (8) healthy controls.

Quality assessment and data extraction

The quality of included studies was assessed using the revised Little’s recommendation based on the following criteria16 17: (1) scientific design, (2) definite inclusion of study population, (3) explicit information on study population, (4) explicit diagnostic criteria of ATDILI, (5) genetic detection method, (6) appropriate statistical analysis and (7) logical discussion of study bias. Studies with an overall score of ≥4 (range 0 to 7) were considered high quality and retained in the analysis.

The following data were extracted from each study using a standardised extraction form: (1) name of the first author, (2) year of publication, (3) the polymorphic gene(s) and genotype(s) under investigation, (4) ethnicity, (5) sample size, (6) mean or median age, (7) sex distribution, (8) anti-tuberculosis drug regimens, (9) diagnostic criteria of ATDILI, (10) genotyping methods and (11) the number of cases and controls for each polymorphic genotype.

Statistical analysis

The genotypes were analysed based on the following proposed genetic risk model: (1) NAT2 (slow acetylator vs intermediate and fast acetylator), (2) CYP2E1 (c1/c1 vs c1/c2 and c2/c2 for the RsaI/PstI polymorphism, D/D vs D/C and C/C for the DraI polymorphism), (3) GSTM1 (null vs non-null), (4) GSTT1 (null vs non-null), (5) GSTM1/GSTT1 (dual-null vs one-null or non-null) and (6) SLCO1B1 388A>G and 521T>C polymorphisms. The genetic risk models for NAT2, CYP2E1, GSTM1, GSTT1 and GSTM1/GSTT1 have been studied in previous studies.9 18 19 Based on these previous studies, the most clinically significant and plausible model for each polymorphic gene was selected. Due to the relative paucity of data suggesting the most clinically relevant genetic model for SLCO1B1 388A>G and 521T>C polymorphisms, all three genomic models including dominant, recessive and additive models were evaluated. The Mantel-Haenszel or DerSimonian-Laird method based on fixed-effects or random-effects models, respectively, were used depending on the presence of heterogeneity.20 21 The random-effects model was used in the presence of significant heterogeneity; otherwise, the fixed-effects model was used to estimate the total effect of a polymorphic gene genotype on the risk of ATDILI. Heterogeneity of study outcomes among included studies was evaluated using Cochran’s Q test (Q) and quantified using Higgin’s I2 test. Significant heterogeneity was defined as the I2 score of >40% accompanied by p<0.10 from the Cochran’s Q test.22 The strength of the association between the genetic polymorphisms and the risk of ATDILI was estimated using pooled ORs with the corresponding 95%CIs. The statistical significance of an OR was defined as p<0.05 from the Z test.

Subgroup analysis was performed based on ethnicity, anti-tuberculosis drug regimen used and the type of study design. Sensitivity analysis was conducted to assess the robustness of the results and to identify the source of heterogeneity using the leave-one-out method. In each analysis, one study was deleted, and with the one study left out, the meta-analysis was performed; this process was repeated until every study had been deleted from our included study pool for each tested polymorphic gene. Publication bias was evaluated with a symmetrical funnel plot. Statistical analyses were performed using Review Manager Software V.5.3 (Cochrane Collaboration, London, UK).

Patient and public involvement

Patients and public were not involved in the design of this study.

ResultsStudy selection and characteristics

Overall, 388 articles were identified through electronic database search and three articles through manual search by reviewing the reference lists of retrieved articles. After removing 99 duplicates, 289 articles were screened for relevance based on the title and abstract. Among them, 72 relevant articles were assessed for eligibility through full-text evaluations. Finally, a total of 54 articles which met the inclusion criteria were included in our analysis (figure 1). Among the 54 studies, 26 studies were included for CYP2E1, 35 studies for NAT2, 19 studies for GST (19 for GSTM1, 17 for GSTT1 and 11 for GSTM1/GSTT1) and four studies for SLCO1B1 388A>G and 521T>C.

F1

Study selection process flowchart. ATDILI, anti-tuberculosis drug-induced liver injury; CYP2E1, cytochrome P450 2E1; GSTM1, glutathione S-transferase Mu 1; GSTT1, glutathione S-transferase Theta 1; NAT2, N-acetyltransferase 2; SLCO1B1, solute carrier organic anion transporter family member 1B1; TB, tuberculosis.

Table 1 summarises the characteristics of the included studies. Across the included studies, large variability in study population was observed in terms of ethnicity (Chinese, Japanese, Korean, Indian, Taiwanese, Brazilian, Caucasian, Iranian, Tunisian and Turkish), age (mean or median age ranging from 27 to 70 years) and sex (the proportion of males ranging from 13% to 90%). Patients in our included studies received either monotherapy with INH or RIF or a combination therapy including a four-drug regimen of INH, RIF, PZA and EMB for the treatment of tuberculosis. ATDILI was defined as an elevated serum alanine aminotransferase concentration by 1.5-fold to 5-fold or greater above the upper limit of normal depending on the study. The quality score of the included studies was 6 or greater based on the revised Little’s recommendation (table 1, online supplementary table S3).16 17 Genotype distribution and genotyping method used in the included studies are summarised for each polymorphic gene in online supplementary table S4 to S7. Funnel plots for CYP2E1, NAT2, GST and SLCO1B1 are provided in online supplementary figure S8. None of the funnel plots showed an asymmetric inverted funnel shape, indicating the absence of potential publication bias.

T1

Characteristics of the studies included in the meta-analysis (n=54 studies)

Last name of the first author, yearPolymorphic geneStudy designEthnicitySample size (case/control)Age (years) (case/control)*Male (%) (case/control)Anti-TB drug regimen administeredDiagnostic criteria of ATDILIQuality score†
Feng et al, 201423 CYP2E1 Case-controlChinese173/17348.8/48.668.0/68.0INH, RIF, PZAALT>3 × ULN6
Kim et al, 200948 CYP2E1 Case-controlKorean67/15942.1/42.865.7/65.4INH, RIF, PZA, EMBALT>2 × ULN7
Singh et al, 201449 CYP2E1 CohortIndian50/135NA/NANA/NANAALT>2 ×  ULN7
Tang et al, 201350 CYP2E1 Case-controlChinese89/35643.7/43.673.0/73.0INH, RIF, PZA, EMBALT>2 × ULN7
Ben Mahmoud et al, 201251 NAT2 CohortTunisian14/5242.4/42.142.8/48.1INH, RIF containing regimenALT>2 × ULN7
Bozok Cetintaş et al, 200852 NAT2 Case-controlTurkish30/7039.8/37.350.0/72.8INH, RIF, PZA, EMBALT>3 × ULN6
Higuchi et al, 200753 NAT2 CohortJapanese18/8260.8/64.750.0/57.3INH, RIF containing regimenALT>2 × ULN7
Ho et al, 201354 NAT2 CohortTaiwanese20/328NA/NANA/NAINH, RIF, PZA, EMBALT>5 × ULN6
Huang et al, 200255 NAT2 CohortTaiwanese33/19173.3/63.787.9/88.5INH, RIF, PZA, EMBALT>2 × ULN7
Khalili et al, 201156 NAT2 Case-controlIranian14/36NA/NANA/NAINH, RIF, PZA, EMBALT>3 × ULN6
Leiro-Fernandez et al, 201157 NAT2 Case-controlCaucasian50/6734.0/30.5‡54.0/56.7INH, RIF, PZAALT>3 × ULN7
Lv et al, 201224 NAT2 Case-controlChinese89/35642.0/42.0‡73.0/73.0INH, RIF, PZA, EMBALT>2 ×  ULN7
Ng et al, 201458 NAT2 Case-controlMixed26/10148.3/NA38.5/NAINH containing regimenALT>5 × ULN7
Ohno et al, 200059 NAT2 CohortJapanese14/63NA/NANA/NAINH, RIFALT>1.5 × ULN7
Possuelo et al, 200860 NAT2 CohortBrazilian14/24038.9/36.550.0/66.9INH, RIF, PZAALT>3 × ULN7
Rana et al, 201261 NAT2 CohortIndian50/20145.3/43.876.0/57.2INH, RIF, PZA, EMBALT>5 × ULN7
Shimizu et al, 200662 NAT2 Case-controlJapanese10/3260.5/64.970.0/46.9INH, RIFALT>2 × ULN6
Yuliwulandari et al, 201663 NAT2 Case-controlIndonesian50/191NA/NANA/NANAALT>2 × ULN7
Wattanapokayakit et al, 201625 NAT2 Case-controlThailand53/8551.4/50.258.5/60.0INH containing regimenALT>2 × ULN7
Chatterjee et al, 201064 GSTM1, GSTT1 Case-controlIndian51/10037.2/33.249.0/63.0INH, RIF, PZAALT>3 × ULN7
Gupta et al, 201365 GSTM1, GSTT1 CohortIndian50/24637.0/36.5‡48.0/56.5INH, RIF, PZA, EMBALT>2 × ULN7
Huang et al, 200766 GSTM1, GSTT1 Case-controlTaiwanese63/6362.0/NANA/NANAALT>5 × ULN6
Kim et al, 201067 GSTM1, GSTT1 Case-controlKorean57/19047.3/42.459.6/67.9INH, RIF, PZA, EMBALT>3 × ULN7
Leiro et al, 200868 GSTM1, GSTT1 Case-controlCaucasian35/6034.0/31.0‡40.0/41.7INH, RIF, PZAALT>3 × ULN7
Liu et al, 201469 GSTM1, GSTT1 Case-controlChinese20/14335.9/61.260.0/59.4INH containing regimenALT>2 × ULN7
Monteiro et al, 201270 GSTM1, GSTT1 CohortBrazilian59/11837.0/38.0‡76.0/61.0NAALT>2 × ULN7
Rana et al, 201371 GSTM1, GSTT1 CohortIndian30/22043.6/42.360.0/64.5INH, RIFALT>5 × ULN6
Roy et al, 200172 GSTM1, GSTT1 Case-controlIndian33/33NA/NANA/NAINH, RIF, PZA, EMBALT>2 × ULN7
Chen et al, 201542 SLCO1B1 Case-controlChinese89/35643.7/43.673.0/73.0INH, RIF, PZA, EMBALT>2 ×  ULN7
Kim et al, 201210 SLCO1B1 Case-controlKorean67/15943.0/42.865.7/65.4INH, RIF, PZA, EMBALT>2 × ULN7
Li et al, 201228 SLCO1B1 Case-controlChinese118/15540.5/39.348.3/54.8RIFALT>3 × ULN7
An et al, 201273 NAT2, CYP2E1 Case-controlChinese101/10736.0/33.4‡55.0/70.0INH, RIF, PZA, EMBALT>2 × ULN7
Bose et al, 201174 NAT2, CYP2E1 CohortIndian41/17738.0/36.043.9/47.4INH, RIF, PZA, EMBALT>2 × ULN7
Chamorro et al, 201375 NAT2, CYP2E1 CohortMixed (South American)47/12829.0/27.041.3/64.8INH, RIF, PZA, EMBALT>3 × ULN7
Cho et al, 200776 NAT2, CYP2E1 CohortKorean18/11451.2/46.766.7/55.3INH, RIF, PZA, EMBALT>2 × ULN7
Gupta et al, 201327 NAT2, CYP2E1 Case-controlIndian50/16537.0/38.048.0/60.0INH, RIF, PZA, EMBALT>2 × ULN7
Huang et al, 200377 NAT2, CYP2E1 CohortTaiwanese49/26970.0/59.0‡18.4/14.9INH, RIF, PZA, EMBALT>2 × ULN7
Lee et al, 201078 NAT2, CYP2E1 CohortTaiwanese45/9558.4/54.960.0/66.3INH, RIF, PZA, EMBALT>2 × ULN7
Mishra et al, 201379 NAT2, CYP2E1 Case-controlIndian33/17338.0/NA52.0/NAINH, RIF, PZA, EMBALT>3 × ULN7
Santos et al, 201380 NAT2, CYP2E1 Case-controlBrazilian18/25247.7/45.656.0/49.0INH, RIFALT>3 × ULN7
Vuilleumier et al, 200681 NAT2, CYP2E1 Case-controlMixed8/6327–35: 2/22§ >36: 5/18§38.0/51.0INHAST or ALT>4 ×  ULN7
Yamada et al, 200982 NAT2, CYP2E1 Case-controlMixed23/147NA/NA13.0/42.9INHALT>2 × ULN7
Zaverucha-do-Valle et al, 201483 NAT2, CYP2E1 CohortBrazilian50/79<40: 28/43§ >40: 20/36§60.4/72.2INH, RIF, PZAALT>2 × ULN6
Sharma et al, 201484 CYP2E1, GSTM1 CohortIndian105/18535.2/27.655.7/72.1INH, RIF, PZA, EMBALT>5 × ULN7
Wang et al, 201085 CYP2E1, GSTM1 Case-controlChinese104/11148.6/44.767.3/67.6INH, RIF, PZA, EMBALT>2 × ULN7
Tang et al, 201286 CYP2E1, GSTM1, GSTT1 Case-controlChinese89/35643.7/43.673.0/73.0INH, RIF, PZA, EMBALT>2 × ULN7
Yimer et al, 201145 NAT2, SLCO1B1 CohortEthiopian41/160NA/NANA/NAINH, RIF, PZA, EMBALT>2 × ULN6
Brito et al, 201487 NAT2, CYP2E1, GSTM1, GSTT1 CohortBrazilian15/23038.1/36.846.7/NAINH, RIF, PZAALT>3 × ULN7
Forestiero et al, 201388 NAT2, CYP2E1, GSTM1, GSTT1 CohortBrazilian59/40NA/NA49.2/60.0INH, RIF, PZAALT>2.5 × ULN6
Rana et al, 201489 NAT2, CYP2E1, GSTM1, GSTT1 CohortIndian55/24543.6/42.360.0/62.0INH, RIF, PZA, EMBALT>5 × ULN6
Singla et al, 201426 NAT2, CYP2E1, GSTM1, GSTT1 Case-controlIndian17/39148.2/32.764.7/61.4INH, RIF, PZA, EMB, STMALT>2 ×  ULN7
Sotsuka et al, 201190 NAT2, CYP2E1, GSTM1, GSTT1 Case-controlJapanese20/9254.9/50.490.0/73.9INH, RIF, PZA, EMB or STMALT>3 × ULN7
Teixeira et al, 201191 NAT2, CYP2E1, GSTM1, GSTT1 Case-controlBrazilian26/14147.6/43.061.5/52.5INH containing regimenALT>3 × ULN7
Xiang et al, 201492 NAT2, CYP2E1, GSTM1, GSTT1 CohortChinese89/215537.0/44.567.4/55.7INH, RIF, PZA, EMBALT>2 × ULN7

*Mean unless otherwise stated.

†Detailed scoring for each quality assessment criterion based on the revised Little’s recommendation in online supplementary table S2.

‡Median age.

§Number of individuals in the age ranges.

ALT, alanine aminotransferase; AST, aspartate aminotransferase; ATDILI, anti-tuberculosis drug-induced liver injury; CYP2E1, cytochrome P450 2E1; EMB, ethambutol; GSTM1, glutathione S-transferase Mu 1; GSTT1, glutathione S-transferase Theta 1; INH, isoniazid; NA, not available; NAT2, N-acetyltransferase 2; PZA, pyrazinamide; RIF, rifampicin; SLCO1B1, solute carrier organic anion transporter family member 1B1 (encoding organic anion transporting polypeptide 1B1 (OATP1B1)); STM, streptomycin; TB, tuberculosis; ULN, upper limit of normal.

CYP2E1

For the CYP2E1 RsaI/PstI polymorphism, 24 studies with 1293 cases and 5450 controls were included in our primary analysis. Using the random-effects model, the pooled estimates of all included studies (n=24) showed a significant association between the risk of ATDILI and the CYP2E1 RsaI/PstI polymorphism (OR for the c1/c1 genotype=1.39, 95% CI 1.06 to 1.83, p=0.02; I2=60%, P heterogeneity  <0.0001) (figure 2A). In the subgroup analysis based on ethnicity and anti-tuberculosis drug regimens, the risk of ATDILI was significantly increased for the CYP2E1 RsaI/PstI c1/c1 genotype in East Asian patients (OR=1.62, 95% CI 1.26 to 2.36, p=0.01; I2=69%, P heterogeneity =0.0006) and in patients receiving a combination of anti-tuberculosis medications (OR=1.35, 95% CI 1.01 to 1.79, p<0.00001; I2=61%, P heterogeneity =0.0002) (online supplementary table S9). No significant association was observed between the risk of ATDILI and the CYP2E1 RsaI/PstI c1/c1 genotype when evaluating studies with the same study design only (ie, either case-control studies or cohort studies) (online supplementary table S9).

F2

Risk of anti-tuberculosis drug-induced liver injury in patients with the CYP2E1 (A) RsaI/PstI c1/c1 genotype compared with c1/c2+c2/c2 genotypes and (B) DraI D/D genotype compared with D/C+C/C genotypes. CYP2E1, cytochrome P450 2E1.

In our primary analysis for the CYP2E1 DraI polymorphism with six studies including 233 cases and 1272 controls, no significant association was observed using the fixed-effects model between the risk of ATDILI and the DraI polymorphism (OR for the D/D genotype=0.93, 95% CI 0.68 to 1.27, p=0.64; I2=0%, P heterogeneity =0.51) (figure 2B).

NAT2

Overall, 35 studies with 1323 cases and 7319 controls were included in our primary analysis for the NAT2 polymorphism. Using the random-effects model, the pooled estimates of all included studies (n=35) showed a significant association between the risk of ATDILI and the NAT2 polymorphism (OR for the slow acetylator genotype=3.30, 95% CI 2.65 to 4.11, p<0.00001; I2=54%, P heterogeneity  <0.0001) (figure 3). In the subgroup analysis based on ethnicity, anti-tuberculosis drug regimens used, and study design, the risk of ATDILI was significantly increased in slow acetylators compared with fast or intermediate acetylators in all subgroups (online supplementary table S10).

F3

Risk of anti-tuberculosis drug-induced liver injury in patients with the NAT2 slow acetylator genotype compared with those with the intermediate/fast acetylator genotypes. NAT2, N-acetyltransferase 2.

GST

For the GSTM1 polymorphism, a total of 19 studies with 977 cases and 5119 controls were included in our primary analysis. Using the fixed-effects model, the pooled estimates of all included studies (n=19) showed a significant association between the risk of ATDILI and the GSTM1 polymorphism (OR for the GSTM1 null genotype=1.30, 95% CI 1.12 to 1.52, p=0.0007; I2=33%, P heterogeneity =0.08) (figure 4A). When studies were stratified for ethnicity, the risk of ATDILI was significantly increased for the GSTM1 null genotype in Indians (OR=1.68, 95% CI 1.30 to 2.19, p<0.0001; I2=36%, P heterogeneity =0.15) (online supplementary table S11). In the subgroup analyses by study design, the estimated OR (95% CI, p value; I2, P heterogeneity ) for the GSTM1 null genotype relative to the non-null genotype was 1.41 (1.04 to 1.93, p=0.03; I2=44%, P heterogeneity =0.08) in cohort studies and 1.25 (1.01 to 1.55, p=0.20; I2=29%, P heterogeneity =0.17) in case-control studies, respectively (online supplementary table S11).

F4

Risk of anti-tuberculosis drug-induced liver injury in patients with (A) the GSTM1 null genotype compared with the non-null genotype, (B) the GSTT1 null genotype compared with the non-null genotype and (C) the GSTM1/GSTT1 dual-null genotype compared with the one-null and non-null genotypes. GSTM1, glutathione S-transferase Mu 1; GSTT1, glutathione S-transferase Theta 1.

For the GSTT1 and GSTM1/GSTT1 polymorphisms, 17 studies (768 cases, 4823 controls) and 11 studies (547 cases, 4233 controls) were included in our primary analyses, respectively. The risk of ATDILI was not significantly associated with the GSTT1 polymorphism (OR for the null genotype=1.03, 95% CI 0.85 to 1.25, p=0.76; I2=16%, P heterogeneity =0.26) or the GSTM1/GSTT1 polymorphism (OR for the dual-null genotype=1.05, 95% CI 0.67 to 1.62, p=0.84; I2=59%, P heterogeneity =0.006) (figure 4B and C). When studies were stratified for ethnicity, anti-tuberculosis drug regimens used, and study design, no subgroups showed significant association between the risk of ATDILI and the GSTT1 and the GSTM1/GSTT1 polymorphisms (online supplementary table S11).

SLCO1B1

For the SLCO1B1 388A>G polymorphism, four studies with 302 cases and 913 controls were included in our primary analysis. Using the dominant, recessive or additive genomic model, no significant association was observed between the risk of ATDILI and the SLCO1B1 388A>G polymorphism (table 2). For the SLCO1B1 521T>C polymorphism, four studies with 314 cases and 912 controls were included in our primary analysis. No significant association was found between the ATDILI risk and the SLCO1B1 521T>C polymorphism under the dominant, recessive or additive genetic model (table 2). Due to the lack of significant association between the risk of ATDILI and the tested SLCO1B1 genetic polymorphisms in our primary meta-analysis, subgroup analyses were not performed for these genetic polymorphisms.

T2

Association between the SLCO1B1 polymorphisms and the risk of anti-tuberculosis drug-induced liver injury

PolymorphismGenetic modelNumber of studiesOR (95% CI)P valueI2, %P heterogeneity Model of meta-analysis
SLCO1B1 388A>G (rs2306283)Dominant modelAA+AG vs GG41.00 (0.76 to 1.31)1.0000.73Fixed
Recessive modelAA vs AG+GG41.45 (0.93 to 2.25)0.1000.84Fixed
Additive modelAA vs GG41.36 (0.85 to 2.15)0.2000.98Fixed
SLCO1B1 521T>C (rs4149056)Dominant modelCC+TC vs TT40.74 (0.43 to 1.28)0.28660.03Random
Recessive modelCC vs TC+TT41.21 (0.40 to 3.64)0.7300.57Fixed
Additive modelCC vs TT41.27 (0.42 to 3.84)0.6700.61Fixed

SLCO1B1, solute carrier organic anion transporter family member 1B1.

Sensitivity analysis

Our primary analysis results showed significantly high heterogeneity between studies for CYP2E1 RsaI/PstI (I2=60%, p<0.0001), NAT2 (I2=54%, p<0.0001), GSTM1/GSTT1 (I2=59%, p=0.006) and SLCO1B1 521T>C (dominant genetic model: I2=66%, p=0.03) polymorphisms. This high heterogeneity between studies may be due to substantial differences in ethnicity, anti-tuberculosis drug regimen, the genotyping method used, study design and diagnostic criteria of ATDILI among the included studies (table 1). Through the sensitivity analyses, outlier studies were identified as the major source of heterogeneity. After removing these outlier studies, heterogeneity was substantially reduced (I2=60% to 42% for CYP2E1 RsaI/PstI,23 I2=54% to 34% for NAT2 24 25, I2 =59% to 0% for GSTM1/GSTT126 27 and I2=66% to 0% for SLCO1B1 521T>C dominant genetic model28). The overall results for the association between the risk of ATDILI and these genetic polymorphisms after excluding the outlier studies stayed the same as those from our primary analysis results.

Discussion

In this study, we conducted a large-scale meta-analysis evaluating the association between the risk of ATDILI and genetic polymorphisms of SLCO1B1 as well as various DMEs including CYP2E1, NAT2 and GST to provide more updated, comprehensive and compelling evidence. Compared with previous meta-analyses, our present study included a larger number of studies, which may sufficiently increase the statistical power compared with individual studies. However, a limited number of studies for the SLCO1B1 genetic polymorphisms were included (n=4). Consistently with previous studies, our current study suggested a significantly increased risk of ATDILI in patients with the CYP2E1 RsaI/PstI c1/c1 genotype (OR=1.39, 95% CI 1.06 to 1.83), the NAT2 slow acetylator genotype (OR=3.30, 95% CI 2.65 to 4.11) and the GSTM1 null genotype (OR=1.30, 95% CI 1.12 to 1.52).9 12 29 Among these genotypes, the largest increase in the risk of ATDILI was shown in patients with the NAT2 slow acetylator genotype. In contrast, no significant association was observed between the risk of ATDILI and the genetic polymorphisms of CYP2E1 DraI, GSTT1, GSTM1/GSTT1, SLCO1B1 388A>G and SLCO1B1 521T>C. Caution needs to be exercised when interpreting this study finding because the lack of significant association between these polymorphisms and the risk of ATDILI might be due to small sample sizes and the low frequency of ATDILI reported in patients with these genetic polymorphisms.

When evaluating the impact of the CYP2E1 RsaI/PstI and DraI genetic polymorphisms on the risk of ATDILI in our study, patients with the RsaI/PstI c1/c1 genotype were 1.39-times more likely to develop ATDILI. Similarly, in a previous meta-analysis by Deng and colleagues, the risk of ATDILI was 1.4-times higher in patients with the RsaI/PstI c1/c1 genotype compared with other genotypes.30 In the liver, INH is metabolised by NAT2 to acetylisoniazid which is consequently oxidised by CYP2E1 to reactive hepatotoxic intermediates.31 32 The increased inducibility or greater activity of CYP2E1 in patients with the CYP2E1 RsaI/PstI c1/c1 genotype may result in the production of more intermediate hepatotoxins, ultimately leading to the increased risk of ATDILI.31 32 Our subgroup analysis showed a significantly increased risk of ATDILI in the CYP2E1 RsaI/PstI c1/c1 genotype carriers of East Asian ethnicity (S9 Table), suggesting a potential gene-ethnicity interaction.33 A previous study identified age, female sex, white race, non-Hispanic ethnicity, lower body mass index, elevated plasma aspartate transaminase concentrations at baseline and 9 months of daily INH use as risk factors for ATDILI.34 Considering their race, ethnicity and relatively lower body mass index compared with other ethnicities, East Asians may be at an increased risk of ATDILI. As the CYP2E1 RsaI/PstI c1 allele frequency is relatively low in this population (79.8% vs 88.5% to 99.8% in other ethnicities), the ethnicity itself might play an important role in developing hepatotoxicity through gene-ethnicity interaction.35 Furthermore, the relatively high frequency of c2 allele in this population might serve as a good control to estimate the effect of c1 allele on the risk of ATDILI; the rarity of this minor allele in other ethnicities could make it difficult to evaluate the association between the ATDILI risk and this genetic polymorphism.35 In addition to ethnicity, combination anti-tuberculosis therapy was shown to significantly increase the risk of ATDILI in patients with the CYP2E1 RsaI/PstI c1/c1 genotype (S9 Table). This is consistent with previous study findings because hepatotoxicity commonly occurs with anti-tuberculosis drugs such as INH and RIF and thus, use of more than one hepatotoxic anti-tuberculosis medication increases the risk of ATDILI.7

Similar to previous studies, our current study suggested a significantly increased risk of ATDILI in patients with the NAT2 slow acetylator genotype compared with those with intermediate/fast acetylator genotypes.9 29 The risk of ATDILI in slow acetylators remained significantly increased in all tested subgroups regardless of ethnicity and the anti-tuberculosis drug regimen used (S10 Table). The frequencies of NAT2 slow acetylator alleles are highly variable between ethnic groups, ranging from 32% in Koreans to 76% in Caucasians.36 Despite this large inter-ethnic variability in the NAT2 polymorphic allele frequency, the NAT2 slow acetylator genotype consistently and significantly increased the risk of ATDILI across all ethnicities, suggesting the critical role of NAT2 polymorphism in the development of ATDILI. In addition, the increased risk of ATDILI in slow acetylators receiving INH monotherapy or combination therapy further highlights the importance of the NAT2 polymorphism in the development of INH-induced hepatotoxicity. The clearance of INH is slower in slow acetylators compared with rapid or intermediate acetylators, resulting in the accumulation of INH in these patients.37 38 This high level of INH may increase the risk of ATDILI in patients with tuberculosis carrying NAT2 slow acetylator genotype due to immune-mediated liver injury through the binding of INH to liver proteins.39 Therefore, clinicians should closely monitor patients with tuberculosis carrying the NAT2 slow acetylator genotype for hepatotoxicity when INH-based treatment is administered to these patients.

According to previous studies, GST enzymes, particularly those coded by GSTM1 and GSTT1 loci, are associated with the risk of drug-induced hepatotoxicity.9 40 Similar to previous studies, our current study demonstrated a significantly increased risk of ATDILI in individuals with the GSTM1 null genotype compared with those with the non-null genotype; however, the risk of ATDILI was not affected by the GSTT1 or GSTM1/GSTT1 genetic polymorphisms. GSTs are important enzymes to detoxify various xenobiotics and play an essential role in INH metabolism by eliminating acetyldiazene ketene acetylonium ion, which is a possibly hepatotoxic free radical metabolite of INH, from the body through GSTM1. This may account for the significant association of the ATDILI risk with the GSTM1 genotype, but not with the GSTT1 or GSTM1/GSTT1 genotypes.9 40 Our subgroup analysis showed a significantly increased risk of ATDILI in the GSTM1 null genotype carriers of Indian ethnicity; although not statistically significant, the risk of ATDILI was relatively high in the East Asian population with the GSTM1 null genotype (online supplementary table S11). Considering the substantial difference in the GSTM1 null allele frequency between Indians (29.6%) and East Asians (52.1%), a potential gene-ethnicity interaction may exist based on their race, ethnicity and body size as aforementioned.34 41 Other characteristics than the GSTM1 polymorphism in these ethnicities may play a more important role in the development of ATDILI. In addition, when studies were stratified by study design, the risk of ATDILI was significantly increased in patients with the GSTM1 null genotype for cohort studies only, but not for case-control studies, probably due to a relatively larger sample size with cohort studies compared with case-control studies.

SLCO1B1 encodes OATP1B1 which is a major influx drug transporter responsible for the hepatic uptake of various endogenous and exogenous substances including RIF.42 Previous studies showed significantly altered systemic exposure of RIF in carriers of the SLCO1B1 polymorphism.43 44 To our knowledge, only four studies have been conducted to examine the association between the ATDILI risk and the SLCO1B1 genetic polymorphisms.10 28 42 45 Various single nucleotide polymorphisms of SLCO1B1 were evaluated in these studies; however, SLCO1B1 388A>G (rs2306283) and 521T>C (rs4149056) were the only polymorphisms assessed in common.10 28 42 45 Therefore, to maximise the sample size in our current meta-analysis, we examined the association between the risk of ATDILI and the polymorphic genotypes of SLCO1B1 388A>G and 521T>C. Similar to each of the included studies, we did not find significant difference in the risk of ATDILI among patients with different SLCO1B1 388A>G and 521T>C genotypes.

There are limitations to this study. First, due to the lack of information regarding other patient characteristics potentially associated with ATDILI, our estimated ORs were not adjusted based on the potential risk factors such as age, anti-tuberculosis drug dosages, alcohol consumption, cigarette smoking and other lifestyle characteristics.7 46 Second, our literature search limited to the articles published in English may lead to language bias. Third, a specific causative agent of ATDILI could not be identified in our analysis because most patients in our included studies received a combination regimen of anti-tuberculosis drugs. Fourth, only the limited number of polymorphic genotypes were assessed for the association with the risk of ATDILI, particularly for SLCO1B1. In addition, only one genetic model was used for CYP2E1, NAT2 and GST when evaluating the association between genetic polymorphisms of these genes and the risk of ATDILI. Although we acknowledge dominant, recessive and additive genomic models can be used for two alleles, it could not be applied to our meta-analysis because we compared patients with different genotype-based phenotype, that is, slow acetylator versus fast/intermediate acetylator and null vs. non-null GSTs. Multiple allelic variants or allele subgroups may represent the same phenotype (eg, NAT2*5B, *6A, and *7B all represent slow acetylator genotypes), and the genetic model selection can be varied depending on the specific allelic variant.47 Therefore, the genetic models used in previous original and meta-analysis studies were adopted for these polymorphic genes in our current study.9 18 19 Future studies are needed to comprehensively and adequately address the relationship between the ATDILI risk and various genetic polymorphisms by using different genetic risk models and including more polymorphic genotypes.

In conclusion, the risk of ATDILI during tuberculosis therapy was significantly increased in patients with tuberculosis carrying NAT2 slow acetylator, CYP2E1 RsaI/PstI c1/c1, or GSTM1 null genotypes. Screening for these genetic polymorphisms, particularly for the NAT2 slow acetylator genotype, may be of great clinical benefit to identify patients at high risk for ATDILI and minimise the risk of ATDILI. Future studies are pertinent to develop dose and/or treatment adjustment strategies, to evaluate the feasibility and cost-effectiveness of the genetic screening test, and to assess the effect of more genetic polymorphisms on the risk of ATDILI for optimal prevention and management of ATDILI.

This study was supported by Research Resettlement Fund for the new faculty of Seoul National University. The authors thank Ms Kyeoul Jeong for her assistance in data management.

EKC and JIL contributed equally.

SY devised and designed the study. SY, JYP and SJH conducted the literature search, performed data extraction and analysis and interpreted the data. SY, EKC and JIL prepared and reviewed the manuscript. All authors reviewed, amended and approved the submitted manuscript.

Research Resettlement Fund for the new faculty of Seoul National University provided funding to JIL. The funder had no role in data collection, data analysis, translation and drafting of the manuscript.

None declared.

Not commissioned; externally peer reviewed.

No additional data are available.

Not required.

References WHO. Global tuberculosis report Geneva. Switzerland: WHO, 2014. WHO. Treatment of tuberculosis guidelines. 4th edition. Geneva, Switzerland: WHO, 2010. Blumberg HM , Burman WJ , Chaisson RE , et al . American Thoracic Society, Centers for Disease Control and Prevention and the Infectious Diseases Society. American Thoracic Society/Centers for Disease Control and Prevention/Infectious Diseases Society of America: treatment of tuberculosis. Am J Respir Crit Care Med 2003;167:603.doi:10.1164/rccm.167.4.603 Zumla A , Raviglione M , Hafner R , et al . Tuberculosis. N Engl J Med Overseas Ed 2013;368:745–55.doi:10.1056/NEJMra1200894 Lange C , Abubakar I , Alffenaar JW , et al . TBNET. Management of patients with multidrug-resistant/extensively drug-resistant tuberculosis in Europe: a TBNET consensus statement. Eur Respir J 2014;44:23–63.doi:10.1183/09031936.00188313 Tostmann A , Boeree MJ , Aarnoutse RE , et al . Antituberculosis drug-induced hepatotoxicity: concise up-to-date review. J Gastroenterol Hepatol 2008;23:192–202.doi:10.1111/j.1440-1746.2007.05207.x Saukkonen JJ , Cohn DL , Jasmer RM , et al . ATS (American Thoracic Society) Hepatotoxicity of Antituberculosis Therapy Subcommittee. An official ATS statement: hepatotoxicity of antituberculosis therapy. Am J Respir Crit Care Med 2006;174:935–52.doi:10.1164/rccm.200510-1666ST Yee D , Valiquette C , Pelletier M , et al . Incidence of serious side effects from first-line antituberculosis drugs among patients treated for active tuberculosis. Am J Respir Crit Care Med 2003;167:1472–7.doi:10.1164/rccm.200206-626OC Cai Y , Yi J , Zhou C , et al . Pharmacogenetic study of drug-metabolising enzyme polymorphisms on the risk of anti-tuberculosis drug-induced liver injury: a meta-analysis. PLoS One 2012;7:e47769.doi:10.1371/journal.pone.0047769 Kim SH , Kim SH , Lee JH , et al . Polymorphisms in drug transporter genes (ABCB1, SLCO1B1 and ABCC2) and hepatitis induced by antituberculosis drugs. Tuberculosis 2012;92:100–4.doi:10.1016/j.tube.2011.09.007 Roy PD , Majumder M , Roy B . Pharmacogenomics of anti-TB drugs-related hepatotoxicity. Pharmacogenomics 2008;9:311–21.doi:10.2217/14622416.9.3.311 Sun F , Chen Y , Xiang Y , et al . Drug-metabolising enzyme polymorphisms and predisposition to anti-tuberculosis drug-induced liver injury: a meta-analysis. Int J Tuberc Lung Dis 2008;12:994–1002. Chan SL , Chua APG , Aminkeng F , et al . Association and clinical utility of NAT2 in the prediction of isoniazid-induced liver injury in Singaporean patients. PLoS One 2017;12:e0186200.doi:10.1371/journal.pone.0186200 Kerb R . Implications of genetic polymorphisms in drug transporters for pharmacotherapy. Cancer Lett 2006;234:4–33.doi:10.1016/j.canlet.2005.06.051 Cui Y , König J , Leier I , et al . Hepatic uptake of bilirubin and its conjugates by the human organic anion transporter SLC21A6. J Biol Chem 2001;276:9626–30.doi:10.1074/jbc.M004968200 Little J , Bradley L , Bray MS , et al . Reporting, appraising, and integrating data on genotype prevalence and gene-disease associations. Am J Epidemiol 2002;156:300–10.doi:10.1093/oxfordjournals.aje.a000179 Tang N , Deng R , Wang Y , et al . GSTM1 and GSTT1 null polymorphisms and susceptibility to anti-tuberculosis drug-induced liver injury: a meta-analysis. Int J Tuberc Lung Dis 2013;17:17–25.doi:10.5588/ijtld.12.0447 Wang FJ , Wang Y , Niu T , et al . Update meta-analysis of the CYP2E1 RsaI/PstI and DraI polymorphisms and risk of antituberculosis drug-induced hepatotoxicity: evidence from 26 studies. J Clin Pharm Ther 2016;41:334–40.doi:10.1111/jcpt.12388 Li C , Long J , Hu X , et al . GSTM1 and GSTT1 genetic polymorphisms and risk of anti-tuberculosis drug-induced hepatotoxicity: an updated meta-analysis. Eur J Clin Microbiol Infect Dis 2013;32:859–68.doi:10.1007/s10096-013-1831-y Mantel N , Haenszel W . Statistical aspects of the analysis of data from retrospective studies of disease. J Natl Cancer Inst 1959;22:719–48. DerSimonian R , Laird N . Meta-analysis in clinical trials. Control Clin Trials 1986;7:177–88.doi:10.1016/0197-2456(86)90046-2 Higgins JP , Thompson SG . Quantifying heterogeneity in a meta-analysis. Stat Med 2002;21:1539–58.doi:10.1002/sim.1186 Feng FM , Guo M , Chen Y , et al . Genetic polymorphisms in metabolic enzymes and susceptibility to anti-tuberculosis drug-induced hepatic injury. Genet Mol Res 2014;13:9463–71.doi:10.4238/2014.November.11.11 Lv X , Tang S , Xia Y , et al .  NAT2 genetic polymorphisms and anti-tuberculosis drug-induced hepatotoxicity in Chinese community population. Ann Hepatol 2012;11:700–7.doi:10.1016/S1665-2681(19)31446-2 Wattanapokayakit S , Mushiroda T , Yanai H , et al . NAT2 slow acetylator associated with anti-tuberculosis drug-induced liver injury in Thai patients. Int J Tuberc Lung Dis 2016;20:1364–9.doi:10.5588/ijtld.15.0310 Singla N , Gupta D , Birbian N , et al . Association of NAT2, GST and CYP2E1 polymorphisms and anti-tuberculosis drug-induced hepatotoxicity. Tuberculosis 2014;94:293–8.doi:10.1016/j.tube.2014.02.003 Gupta VH , Amarapurkar DN , Singh M , et al . Association of N-acetyltransferase 2 and cytochrome P450 2E1 gene polymorphisms with antituberculosis drug-induced hepatotoxicity in Western India. J Gastroenterol Hepatol 2013;28:1368–74.doi:10.1111/jgh.12194 Li LM , Chen L , Deng GH , et al . SLCO1B1 *15 haplotype is associated with rifampin-induced liver injury. Mol Med Rep 2012;6:75–82.doi:10.3892/mmr.2012.900 Wang P-Y , Xie S-Y , Hao Q , et al . NAT2 polymorphisms and susceptibility to anti-tuberculosis drug-induced liver injury: a meta-analysis [Review article]. The International Journal of Tuberculosis and Lung Disease 2012;16:589–95.doi:10.5588/ijtld.11.0377 Deng R , Yang T , Wang Y , et al . CYP2E1 RsaI/PstI polymorphism and risk of anti-tuberculosis drug-induced liver injury: a meta-analysis. Int J Tuberc Lung Dis 2012;16:1574–81.doi:10.5588/ijtld.12.0304 Mitchell JR , Zimmerman HJ , Ishak KG , et al . Isoniazid liver injury: clinical spectrum, pathology, and probable pathogenesis. Ann Intern Med 1976;84:181–92.doi:10.7326/0003-4819-84-2-181 Ryan DE , Ramanathan L , Iida S , et al . Characterization of a major form of rat hepatic microsomal cytochrome P-450 induced by isoniazid. J Biol Chem 1985;260:6385–93. Ortega VE , Meyers DA . Pharmacogenetics: implications of race and ethnicity on defining genetic profiles for personalized medicine. J Allergy Clin Immunol 2014;133:16–26.doi:10.1016/j.jaci.2013.10.040 Bliven-Sizemore EE , Sterling TR , Shang N , et al . TB Trials Consortium. Three months of weekly rifapentine plus isoniazid is less hepatotoxic than nine months of daily isoniazid for LTBI. Int J Tuberc Lung Dis 2015;19:1039–44. i-v.doi:10.5588/ijtld.14.0829 Zhou Y , Ingelman-Sundberg M , Lauschke VM . Worldwide Distribution of Cytochrome P450 Alleles: A Meta-analysis of Population-scale Sequencing Projects. Clin Pharmacol Ther 2017;102:688–700.doi:10.1002/cpt.690 Lin HJ , Han CY , Lin BK , et al . Ethnic distribution of slow acetylator mutations in the polymorphic N-acetyltransferase (NAT2) gene. Pharmacogenetics 1994;4:125–34.doi:10.1097/00008571-199406000-00003 Seng KY , Hee KH , Soon GH , et al . Population pharmacokinetic analysis of isoniazid, acetylisoniazid, and isonicotinic acid in healthy volunteers. Antimicrob Agents Chemother 2015;59:6791–9.doi:10.1128/AAC.01244-15 Wang P , Pradhan K , Zhong XB , et al . Isoniazid metabolism and hepatotoxicity. Acta Pharm Sin B 2016;6:384–92.doi:10.1016/j.apsb.2016.07.014 Metushi IG , Nakagawa T , Uetrecht J . Direct oxidation and covalent binding of isoniazid to rodent liver and human hepatic microsomes: humans are more like mice than rats. Chem Res Toxicol 2012;25:2567–76.doi:10.1021/tx300341r Ramachandran G , Swaminathan S . Role of pharmacogenomics in the treatment of tuberculosis: a review. Pharmgenomics Pers Med 2012;5:89–98.doi:10.2147/PGPM.S15454 Kasthurinaidu SP , Ramasamy T , Ayyavoo J , et al . GST M1-T1 null allele frequency patterns in geographically assorted human populations: a phylogenetic approach. PLoS One 2015;10:e0118660.doi:10.1371/journal.pone.0118660 Chen R , Wang J , Tang S , et al . Association of polymorphisms in drug transporter genes (SLCO1B1 and SLC10A1) and anti-tuberculosis drug-induced hepatotoxicity in a Chinese cohort. Tuberculosis 2015;95:68–74.doi:10.1016/j.tube.2014.11.004 Weiner M , Peloquin C , Burman W , et al . Effects of tuberculosis, race, and human gene SLCO1B1 polymorphisms on rifampin concentrations. Antimicrob Agents Chemother 2010;54:4192–200.doi:10.1128/AAC.00353-10 Chigutsa E , Visser ME , Swart EC , et al . The SLCO1B1 rs4149032 polymorphism is highly prevalent in South Africans and is associated with reduced rifampin concentrations: dosing implications. Antimicrob Agents Chemother 2011;55:4122–7.doi:10.1128/AAC.01833-10 Yimer G , Ueda N , Habtewold A , et al . Pharmacogenetic & pharmacokinetic biomarker for efavirenz based ARV and rifampicin based anti-TB drug induced liver injury in TB-HIV infected patients. PLoS One 2011;6:e27810.doi:10.1371/journal.pone.0027810 Chen R , Wang J , Zhang Y , et al . Key factors of susceptibility to anti-tuberculosis drug-induced hepatotoxicity. Arch Toxicol 2015;89:883–97.doi:10.1007/s00204-015-1473-1 Khan S , Mandal RK , Elasbali AM , et al . Pharmacogenetic association between NAT2 gene polymorphisms and isoniazid induced hepatotoxicity: trial sequence meta-analysis as evidence. Biosci Rep 2019;39:BSR20180845.doi:10.1042/BSR20180845 Kim SH , Kim SH , Bahn JW , et al . Genetic polymorphisms of drug-metabolizing enzymes and anti-TB drug-induced hepatitis. Pharmacogenomics 2009;10:1767–79.doi:10.2217/pgs.09.100 Singh M , Gupta VH , Amarapurkar DN , et al . Association of genetic variants with anti-tuberculosis drug induced hepatotoxicity: a high resolution melting analysis. Infect Genet Evol 2014;23:42–8.doi:10.1016/j.meegid.2014.01.027 Tang S , Lv X , Zhang Y , et al . Cytochrome P450 2E1 gene polymorphisms/haplotypes and anti-tuberculosis drug-induced hepatitis in a Chinese cohort. PLoS One 2013;8:e57526.doi:10.1371/journal.pone.0057526 Ben Mahmoud L , Ghozzi H , Kamoun A , et al . Polymorphism of the N-acetyltransferase 2 gene as a susceptibility risk factor for antituberculosis drug-induced hepatotoxicity in Tunisian patients with tuberculosis. Pathol Biol 2012;60:324–30.doi:10.1016/j.patbio.2011.07.001 Bozok Cetintaş V , Erer OF , Kosova B , et al . Determining the relation between N-acetyltransferase-2 acetylator phenotype and antituberculosis drug induced hepatitis by molecular biologic tests. Tuberk Toraks 2008;56:81–6. Higuchi N , Tahara N , Yanagihara K , et al . NAT2 6A, a haplotype of the N-acetyltransferase 2 gene, is an important biomarker for risk of anti-tuberculosis drug-induced hepatotoxicity in Japanese patients with tuberculosis. World J Gastroenterol 2007;13:6003–8.doi:10.3748/wjg.v13.45.6003 Ho HT , Wang TH , Hsiong CH , et al . The NAT2 tag SNP rs1495741 correlates with the susceptibility of antituberculosis drug-induced hepatotoxicity. Pharmacogenet Genomics 2013;23:200–7.doi:10.1097/FPC.0b013e32835e95e1 Huang YS , Chern HD , Su WJ , et al . Polymorphism of the N-acetyltransferase 2 gene as a susceptibility risk factor for antituberculosis drug-induced hepatitis. Hepatology 2002;35:883–9.doi:10.1053/jhep.2002.32102 Khalili H , Fouladdel S , Sistanizad M , et al . Association of N-acetyltransferase-2 genotypes and anti-tuberculosis induced liver injury; first case-controlled study from Iran. Curr Drug Saf 2011;6:17–22.doi:10.2174/157488611794479946 Leiro-Fernandez V , Valverde D , Vázquez-Gallardo R , et al . N-acetyltransferase 2 polymorphisms and risk of anti-tuberculosis drug-induced hepatotoxicity in Caucasians. Int J Tuberc Lung Dis 2011;15:1403–8.doi:10.5588/ijtld.10.0648 Ng CS , Hasnat A , Al Maruf A , et al . N-acetyltransferase 2 (NAT2) genotype as a risk factor for development of drug-induced liver injury relating to antituberculosis drug treatment in a mixed-ethnicity patient group. Eur J Clin Pharmacol 2014;70:1079–86.doi:10.1007/s00228-014-1703-0 Ohno M , Yamaguchi I , Yamamoto I , et al . Slow N-acetyltransferase 2 genotype affects the incidence of isoniazid and rifampicin-induced hepatotoxicity. Int J Tuberc Lung Dis 2000;4:256–61. Possuelo LG , Castelan JA , de Brito TC , et al . Association of slow N-acetyltransferase 2 profile and anti-TB drug-induced hepatotoxicity in patients from Southern Brazil. Eur J Clin Pharmacol 2008;64:673–81.doi:10.1007/s00228-008-0484-8 Rana SV , Ola RP , Sharma SK , et al . Comparison between acetylator phenotype and genotype polymorphism of n-acetyltransferase-2 in tuberculosis patients. Hepatol Int 2012;6:397–402.doi:10.1007/s12072-011-9309-4 Shimizu Y , Dobashi K , Mita Y , et al . DNA microarray genotyping of N-acetyltransferase 2 polymorphism using carbodiimide as the linker for assessment of isoniazid hepatotoxicity. Tuberculosis 2006;86:374–81.doi:10.1016/j.tube.2005.09.002 Yuliwulandari R , Susilowati RW , Wicaksono BD , et al . NAT2 variants are associated with drug-induced liver injury caused by anti-tuberculosis drugs in Indonesian patients with tuberculosis. J Hum Genet 2016;61:533–7.doi:10.1038/jhg.2016.10 Chatterjee S , Lyle N , Mandal A , et al . GSTT1 and GSTM1 gene deletions are not associated with hepatotoxicity caused by antitubercular drugs. J Clin Pharm Ther 2010;35:465–70.doi:10.1111/j.1365-2710.2009.01101.x Gupta VH , Singh M , Amarapurkar DN , et al . Association of GST null genotypes with anti-tuberculosis drug induced hepatotoxicity in Western Indian population. Ann Hepatol 2013;12:959–65.doi:10.1016/S1665-2681(19)31302-X Huang YS , Su WJ , Huang YH , et al . Genetic polymorphisms of manganese superoxide dismutase, NAD(P)H:quinone oxidoreductase, glutathione S-transferase M1 and T1, and the susceptibility to drug-induced liver injury. J Hepatol 2007;47:128–34.doi:10.1016/j.jhep.2007.02.009 Kim SH , Kim SH , Yoon HJ , et al . GSTT1 and GSTM1 null mutations and adverse reactions induced by antituberculosis drugs in Koreans. Tuberculosis 2010;90:39–43.doi:10.1016/j.tube.2009.12.001 Leiro V , Fernández-Villar A , Valverde D , et al . Influence of glutathione S-transferase M1 and T1 homozygous null mutations on the risk of antituberculosis drug-induced hepatotoxicity in a Caucasian population. Liver Int 2008;28:835–9.doi:10.1111/j.1478-3231.2008.01700.x Liu F , Jiao A-xia , Wu X-rong , et al . Impact of Glutathione S-Transferase M1 and T1 on Anti-Tuberculosis Drug–Induced Hepatotoxicity in Chinese Pediatric Patients. PLoS One 2014;9:9:e115410.doi:10.1371/journal.pone.0115410 Monteiro TP , El-Jaick KB , Jeovanio-Silva AL , et al . The roles of GSTM1 and GSTT1 null genotypes and other predictors in anti-tuberculosis drug-induced liver injury. J Clin Pharm Ther 2012;37:712–8.doi:10.1111/j.1365-2710.2012.01368.x Rana SV , Kamboj JK , Sharma SK , et al . Antioxidant status and GST gene polymorphisms in antitubercular treatment-induced hepatotoxicity patients. Hepatol Int 2013;7:876–82.doi:10.1007/s12072-013-9452-1 Roy B , Chowdhury A , Kundu S , et al . Increased risk of antituberculosis drug-induced hepatotoxicity in individuals with glutathione S-transferase M1 ’null' mutation. J Gastroenterol Hepatol 2001;16:1033–7.doi:10.1046/j.1440-1746.2001.02585.x An HR , Wu XQ , Wang ZY , et al . NAT2 and CYP2E1 polymorphisms associated with antituberculosis drug-induced hepatotoxicity in Chinese patients. Clin Exp Pharmacol Physiol 2012;39:535–43.doi:10.1111/j.1440-1681.2012.05713.x Bose PD , Sarma MP , Medhi S , et al . Role of polymorphic N-acetyl transferase2 and cytochrome P4502E1 gene in antituberculosis treatment-induced hepatitis. J Gastroenterol Hepatol 2011;26:312–8.doi:10.1111/j.1440-1746.2010.06355.x Chamorro JG , Castagnino JP , Musella RM , et al . Sex, ethnicity, and slow acetylator profile are the major causes of hepatotoxicity induced by antituberculosis drugs. J Gastroenterol Hepatol 2013;28:323–8.doi:10.1111/jgh.12069 Cho HJ , Koh WJ , Ryu YJ , et al . Genetic polymorphisms of NAT2 and CYP2E1 associated with antituberculosis drug-induced hepatotoxicity in Korean patients with pulmonary tuberculosis. Tuberculosis 2007;87:551–6.doi:10.1016/j.tube.2007.05.012 Huang YS , Chern HD , Su WJ , et al . Cytochrome P450 2E1 genotype and the susceptibility to antituberculosis drug-induced hepatitis. Hepatology 2003;37:924–30.doi:10.1053/jhep.2003.50144 Lee SW , Chung LS , Huang HH , et al . NAT2 and CYP2E1 polymorphisms and susceptibility to first-line anti-tuberculosis drug-induced hepatitis. Int J Tuberc Lung Dis 2010;14:622–6. Mishra S , Daschakraborty S , Shukla P , et al . N-acetyltransferase and cytochrome P450 2E1 gene polymorphisms and susceptibility to antituberculosis drug hepatotoxicity in an Indian population. Natl Med J India 2013;26:260–5. Santos NP , Callegari-Jacques SM , Ribeiro Dos Santos AK , et al . N-acetyl transferase 2 and cytochrome P450 2E1 genes and isoniazid-induced hepatotoxicity in Brazilian patients. Int J Tuberc Lung Dis 2013;17:499–504.doi:10.5588/ijtld.12.0645 Vuilleumier N , Rossier MF , Chiappe A , et al . CYP2E1 genotype and isoniazid-induced hepatotoxicity in patients treated for latent tuberculosis. Eur J Clin Pharmacol 2006;62:423–9.doi:10.1007/s00228-006-0111-5 Yamada S , Tang M , Richardson K , et al . Genetic variations of NAT2 and CYP2E1 and isoniazid hepatotoxicity in a diverse population. Pharmacogenomics 2009;10:1433–45.doi:10.2217/pgs.09.66 Zaverucha-do-Valle C , Monteiro SP , El-Jaick KB , et al . The role of cigarette smoking and liver enzymes polymorphisms in anti-tuberculosis drug-induced hepatotoxicity in Brazilian patients. Tuberculosis 2014;94:299–305.doi:10.1016/j.tube.2014.03.006 Sharma SK , Jha BK , Sharma A , et al . Genetic polymorphisms of CYP2E1 and GSTM1 loci and susceptibility to anti-tuberculosis drug-induced hepatotoxicity. Int J Tuberc Lung Dis 2014;18:588–93.doi:10.5588/ijtld.13.0344 Wang T , Yu HT , Wang W , et al . Genetic polymorphisms of cytochrome P450 and glutathione S-transferase associated with antituberculosis drug-induced hepatotoxicity in Chinese tuberculosis patients. J Int Med Res 2010;38:977–86.doi:10.1177/147323001003800324 Tang SW , Lv XZ , Zhang Y , et al . CYP2E1, GSTM1 and GSTT1 genetic polymorphisms and susceptibility to antituberculosis drug-induced hepatotoxicity: a nested case-control study. J Clin Pharm Ther 2012;37:588–93.doi:10.1111/j.1365-2710.2012.01334.x Brito TC , POSSUELO LIAG , Valim ARM , et al . Polymorphisms in CYP2E1, GSTM1 and GSTT1 and anti-tuberculosis drug-induced hepatotoxicity. Anais da Academia Brasileira de Ciências 2014;86:855–65.doi:10.1590/0001-3765201420130350 Forestiero FJ , Cecon L , Hirata MH , et al . Relationship of NAT2, CYP2E1 and GSTM1/GSTT1 polymorphisms with mild elevation of liver enzymes in Brazilian individuals under anti-tuberculosis drug therapy. Clin Chim Acta 2013;415:215–9.doi:10.1016/j.cca.2012.10.030 Rana SV , Sharma SK , Ola RP , et al . N-acetyltransferase 2, cytochrome P4502E1 and glutathione S-transferase genotypes in antitubercular treatment-induced hepatotoxicity in North Indians. J Clin Pharm Ther 2014;39:91–6.doi:10.1111/jcpt.12105 Sotsuka T , Sasaki Y , Hirai S , et al . Association of isoniazid-metabolizing enzyme genotypes and isoniazid-induced hepatotoxicity in tuberculosis patients. In Vivo 2011;25:803–12. Teixeira RL , Morato RG , Cabello PH , et al . Genetic polymorphisms of NAT2, CYP2E1 and GST enzymes and the occurrence of antituberculosis drug-induced hepatitis in Brazilian TB patients. Mem Inst Oswaldo Cruz 2011;106:716–24.doi:10.1590/s0074-02762011000600011 Xiang Y , Ma L , Wu W , et al . The Incidence of Liver Injury in Uyghur Patients Treated for TB in Xinjiang Uyghur Autonomous Region, China, and Its Association with Hepatic Enzyme Polymorphisms NAT2, CYP2E1, GSTM1 and GSTT1. PLoS One 2014;9:e85905:9.doi:10.1371/journal.pone.0085905