- Research article
- Open Access
- Open Peer Review
Clinical and genetic factors associated with increased risk of severe liver toxicity in a monocentric cohort of HIV positive patients receiving nevirapine-based antiretroviral therapy
BMC Infectious Diseasesvolume 18, Article number: 556 (2018)
Nevirapine has been used as antiretroviral agent since early ‘90. Although nevirapine is not currently recommended in initial anti-HIV regimens, its use remains consistent in a certain number of HIV-1-positive subjects. Thus, our aim was to determine clinical and genetic factors involved in the development of severe nevirapine induced liver toxicity.
We retrospectively analyzed all HIV positive patients who were followed at the Infectious Diseases Unit, DIBIC Luigi Sacco, University of Milan from May 2011 to December 2015. All patients treated with nevirapine who underwent a genotyping for the functional variants mapping into ABCB1, CYP2B6, CYP3A4 and CYP3A5 genes were included in the analysis. Severe hepatotoxicity was defined as ACTG grade 3–4 AST/ALT increase during the first three months of nevirapine treatment. The causality assessment between NVP exposure and drug-induced liver injury was performed by using the updated Roussel Uclaf Causality Assessment Methods. Hardy Weinberg equilibrium was tested by χ2 test. A multivariable logistic regression model was constructed using a backward elimination method.
Three hundred and sixty-two patients were included in the analysis, of which 8 (2.2%) experienced a severe liver toxicity. We observed no differences between patients with and without liver toxicity as regards gender, ethnicity, age and immune-virological status. A higher prevalence of HCV coinfection (75.0% vs 30.2%; p = .0013) and higher baseline AST (58 IU/L vs 26 IU/L; p = 0.041) and ALT (82 IU/L vs 27 IU/L; p = 0.047) median levels were observed in patients with liver toxicity vs those without toxicity. The genotypes CT/TT at ABCB1 rs1045642 single nucleotide polymorphism (SNP), showed a protective effect for liver toxicity when compared with genotype CC (OR = 0.18, 95%CI 0.04–0.76; p = 0.020) in univariate analysis. In the multivariate model, HCV coinfection was independently associated with higher risk of developing liver toxicity (aOR = 8.00, 95%CI 1.27–50.29; p = 0.027), whereas ABCB1 rs1045642 CT/TT genotypes (aOR = 0.10, 95%CI 0.02–0.47; p = 0.004) was associated with a lower risk.
According to our findings HCV coinfection and ABCB1 rs1045642 SNP represent independent determinants of severe liver toxicity related to nevirapine. This genetic evaluation could be included as toxicity assessment in HIV-1-positive subjects treated with nevirapine.
Nevirapine (NVP) is a potent non-nucleoside reverse transcriptase inhibitor widely prescribed in low-income countries for HIV treatment and prevention of mother-to-child transmission of HIV . In high-resource countries, NVP is no longer included among antiretrovirals recommended for initial antiretroviral treatment (ART), although it still remains a valid component of regimens used as ART simplification strategy due to its high efficacy, good metabolic profile, convenience, and low cost [2,3,4].
Although generally well tolerated and effective, some individuals exposed to NVP develop in the short-term hypersensitivity reactions which can manifest as hepatotoxicity and/or severe cutaneous adverse reactions . Hepatotoxicity, in particular, has been reported more commonly with NVP than with other antiretroviral drugs [6, 7]. Higher baseline and nadir CD4 cell count have been found to independently influence the risk of NVP-related adverse reactions and the use of NVP in naive women (with CD4 > 250 cells/uL) and males (with CD4 > 400 cells/uL) is disallowed . The role of the immune system, i.e. higher CD4 cells count, in the development of NVP induced skin and liver toxicity is corroborated by the higher incidence of these events in HIV negative patients receiving NVP as a component of post exposure prophylaxis . Nevertheless, when NVP is used in ART-experienced patients with a controlled viremia the risk of development hepatic or cutaneous adverse events significantly decrease and there is no evident association with CD4 cell count [9,10,11].
A series of demographic and clinical factors have been found to correlate with an increased risk of NVP-toxicity. In particular, risk factors for ALT and AST elevations during NVP therapy included alteration of liver enzymes before NVP start, co-infection with hepatitis B or hepatitis C viruses, female gender and low body weight [6, 12].
A number of studies have also investigated the possible predictive role of genetic polymorphisms of CYP enzymes or drug-transporters involved in NVP metabolism in predisposing to NVP-related adverse effects. NVP is metabolized by cytochrome P450 enzymes CYP2B6 and CYP3A4 with a minor contribution from CYP3A5 . Single nucleotide polymorphisms (SNPs) have been found to impact NVP pharmacokinetics in ethnic mix populations [14,15,16]. The genotype TT (c.516/rs3745274) in the CYP2B6 gene, in particular, has been associated with higher plasma concentrations of NVP and its possible role in increasing the risk of hepatotoxicity has been hypothesized . However, there is contrasting evidence of interactions between the presence of variant alleles of CYP2B6 and the development of NVP-induced hepatotoxicity [17,18,19]. The role of CYP3A4 and CYP3A5 variants in determining NVP plasma concentration and the development of liver toxicity is more controversial, with only one report of association between CYP3A5 variants and transaminase values in African patients exposed to NVP . Moreover, although effects of the efflux transporter P-glycoprotein encoded by the ATP Binding Cassette Subfamily B Member 1 (ABCB1) gene on NVP pharmacokinetics remains controversial, two studies have found a protective effect of ABCB1 c.3435 T allele against NVP-related hepatotoxicity [19, 20].
The majority of the studies that evaluated the correlation between pharmacogenetic profiles and NVP-related hepatotoxicity have been performed mainly on African population. However, genetic variant frequencies can differ markedly between different populations and only few data are available on the mentioned pharmacogenetic profiles in non-African populations.
With this in mind, we assessed clinical and pharmacogenetic factors associated with the risk of severe NVP induced liver toxicity in a population of HIV-positive patients attending a clinical center in Italy.
This study was conducted on a cohort of adult HIV-positive patients attending the Infectious Diseases Unit, DIBIC Luigi Sacco, University of Milan between May 1 2011 and December 31 2015. Patients who have ever received or were receiving a NVP-containing cART at our clinical center were eligible for the analysis.
Patient’s demographic (age, gender, and ethnicity), epidemiological (HIV acquisition risk) and clinical (CDC stage, body mass index, coinfections, previous and current antiretroviral regimens, immune-virological and hemato-biochemical parameters) data registered during medical visits (on average every three months) are routinely collected in a structured database, allowing the use of the database for clinical, epidemiological or therapeutic studies.
Severe hepatotoxicity was defined as ACTG grade 3–4 AST and/or ALT increase (AST or ALT elevation above 5 time the upper reference limit) during the first three months of nevirapine treatment. The causality assessment between NVP exposure and drug-induced liver injury was performed by using the updated Roussel Uclaf Causality Assessment Methods (RUCAM) . According to the RUCAM, patients were firstly assessed for hepatocellular, cholestatic or mixed liver injury. Subsequently, the score was applied and single cases of NVP-induced liver injury were classified accordingly to the RUCAM total score interpretation and causality grading: ≤0, excluded; 1–2, unlikely; 3–5, possible; 6–8, probable and ≥ 9 highly probable .
All patients who met the inclusion criteria underwent a genotyping for the functional variants mapping in ABCB1 (c.3435/rs1045642), CYP2B6 (c.516/rs3745274), CYP3A5 (*3/rs776746) and CYP3A4 (*22/rs35599367) genes. Genomic DNA was isolated from peripheral blood cells using an automatic DNA extraction system (Maxwell® 16 System, Promega) according to the manufacturer’s instructions. All genotypes were determined by Real-Time PCR, using a panel of LightSNiP from TIB-MolBiol (assays based on SimpleProbe®). At the end of the amplification a melting curve analysis was performed (LightCycler 480, Roche).
Baseline clinical characteristics and genotypes of the two groups of interest, i.e. patients who developed severe hepatotoxicity and those who did not, were compared using the χ2 or Fisher’s exact test for categorical variables and the Mann-Whitney test for continuous variables.
Hardy Weinberg equilibrium was tested by χ2 test.
The association of clinical and genotypic variables with the development of sever liver toxicity was tested by means of a univariate logistic regression model, and all variables were incorporated into a multivariate logistic regression model with a backward elimination method. Statistical significance was defined at 2-sided P value < 0.05. The risks were expressed as adjusted odds ratios (aOR) with relative confidence intervals (95% CI). To perform statistical analysis we used the SAS software version 9.3.
The study was reviewed and approved by our ethics committee (Comitato Etico Interaziendale, Milano area 1); all subjects signed a dedicated informed consent.
A total of 362 patients were included in the analysis. Most of them were ART-experienced at the time of initiation of NVP-containing regimen, whereas a minority (16.1%) were ART-naïve. Overall 8 (2.2%) patients experienced a severe liver toxicity during the first three months from NVP initiation. Table 1 shows the comparison of patients who experienced a severe hepatotoxicity and those who did not. There was no significant difference as regards gender, ethnicity, age and baseline immune-virological status between the two groups although females showed a trend towards a higher frequency of NVP-induced hepatotoxicity (75.0% vs 35.9%; p = 0.055). Conversely, patients who developed severe hepatotoxicity were more frequently HCV-coinfected (75.0% vs 30.2%: p = 0.013), showed higher baseline AST and ALT median levels (58 IU/L vs 26 IU/L p = 0.041, 82 IU/L vs 27 IU/L p = 0.047, respectively), and a lower median baseline body mass index value (19.9 kg/m2 vs 22.6 kg/m2; p = 0.017).
According to RUCAM, the 8 cases of NVP-induced liver toxicity were classified as hepatocellular injury and the likelihood of NVP-induced liver toxicity resulted possible for 3 patients and probable for 5 patients as shown in Table 2. A brief narration for each of the 8 cases is reported in Table 3.
Distribution of different genotypes of ABCB1 rs1045642, CYP2B6 rs3745274 and CYP3A4/A5 combined are shown in Table 4. A statistical significant difference was observed in patient with and without NVP-induced hepatotoxicity according to ABCB1 rs1045642 genotypes (p = 0.019).
In univariate analysis (Table 5), male gender (OR = 0.19 95%CI 0.04–0.94; p = 0.042), HCV coinfection (OR = 6.93 95%CI 1.38–34.87; p = 0.019), AST (OR = 1.02, 95%CI 1.01–1.03; p = 0.008) and ALT (OR = 1.01, 95%CI 1.00–1.02; p = 0.015) median level at the beginning of NVP have been associated with an increased risk of severe NVP-induced liver toxicity. On the other hand, genotypes CT/TT of ABCB1 rs1045642 (OR = 0.18, 95%CI 0.04–0.76; p = 0.020) and greater value of body mass index (OR = 0.7 95%CI 0.51–0.94, p = 0.02) showed a protective effect.
In the multivariate logistic regression model (Table 5), HCV coinfection was confirmed to be independently associated with a higher risk of developing liver toxicity (aOR = 8.00, 95%CI 1.27–50.29; p = 0.027), whereas ABCB1 CT/TT genotypes (aOR = 0.10, 95%CI 0.02–0.47; p = 0.004) has been associated with a lower risk; a higher body mass index (aOR = 0.72, 95%CI 0.519–1.000; p = 0.050) has been barely related to a lower risk (Table 5). On the contrary, the association between gender, baseline AST levels and NVP-induced liver toxicity wasn’t confirmed in the final multivariate model.
In this study conducted in a monocentric cohort of HIV-1 positive patients exposed to NVP we observed, during the first three months of treatment, an incidence of severe liver toxicity of 2.2%. This finding is similar to that reported by other cohorts when NVP was used in experienced patients [10, 22]. We did not observe a significant association between the development of hepatotoxicity and high CD4 cells count at NVP start, supporting the observation of low frequency of NVP induced liver toxicity in experienced patients . The mechanisms involved in the development of severe hepatotoxicity are not well explained and it could be that in patients never exposed to antiretroviral therapy immune-mediated process leading to immune-reconstitution could elicit the development of liver toxicity . On the contrary, in experienced patients with a stable immune-virological situation, hepatotoxicity could be driven by a direct effect of the drug in susceptible patients . We confirm previous findings supporting the importance of HCV coinfection as independent factor associated to the development of NVP-related liver toxicity; HCV infection could play a direct role causing liver injury and also could interfere with the metabolism of the drug [8, 11, 25]. The enhanced risk of development hepatotoxicity in HCV coinfected patients treated with NVP seems to be independent from NVP plasmatic concentrations, since comparable NVP concentration are observed in patients with and without HCV coinfection .
A significant correlation between low value of body mass index (< 18.5) such as for increased NVP plasma concentration and increased risk for hepatotoxicity has been previously reported [12, 18], in our study we observed a trend of body weight in predisposing to NVP-hepatotoxicity albeit not confirmed in the multivariate model.
In accordance with previous studies, we did not found a statistically significant association between polymorphisms in CYPs genes (CYP2B6, CYP3A4, CYP3A5) and the development of severe hepatotoxicity. Although these genes are involved in NVP metabolism  and their functional variants may significantly affect NVP plasma concentrations [28,29,30], their role in predisposing NVP induced hepatotoxicity remains unclear [19, 31].
The ABCB1 gene encodes for P-glycoprotein, one of the most important efflux pomp involved in the transport of both NVP and efavirenz and a modification of this protein could determine an alteration in the intracellular concentration of these drugs .
Interestingly, the functional variant c.3435 C > T of ABCB1 gene was associated with an increases risk for severe liver toxicity in a previous study conducted by Hass et al. in South-Africa . Moreover, the variant ABCB1 c.3435 C > T resulted protective for NVP-associated hepatic adverse events in another study conducted in Mozambique . No statistically significant association between ABCB1 c.3435 C > T and NVP adverse events was found in another study conducted by Yuan J et al. in an ethnic mixed population. However, these authors found a significant association between ABCB1 c.3435 C > T variant and hepatic adverse events among Africans but not Asians or Caucasians, despite these latter groups showed increased T-allele frequencies .
Overall our study, which was conducted on a prevalently European Caucasian population, disagrees with the study by Yan et al. On the contrary it supports the protective role of T allele of ABCB1 rs1045642 for NVP-hepatotoxicity evidenced by Hass et al. and Ciccacci et al.
The present study has some limitations. In particular, due to the retrospective design we cannot exclude the presence of possible bias related to loss of data. A sub-optimal performance of RUCAM, which is more fitted for a prospective patient evaluation, could be hypothesized because of the retrospective nature of the study. Moreover, if on one hand the limited number of NVP induced hepatic adverse events supports the good safety profile of this drug in ART-experienced patients, on the other hand our study could not exclude that the prevalence of the investigated polymorphisms could be driven by chance. Moreover, the prevalence of Caucasian ethnicity limited the comparison of our findings between different ethnic groups.
Beyond to clinical conditions well known to drive the development of hepatotoxicity during NVP treatment, i.e. HCV coinfection and body mass index, pharmacogenomic profiles could also play a role in this phenomenon. Our results suggest the independent role of ABCB1 rs1045642 as a predictive marker of severe liver toxicity related to NVP. Further validation studies, to assess possible clinical application of this marker in countries in which NVP is still widely used and/or in patients with other risk factors for nevirapine related toxicity, are warranted.
ATP Binding Cassette Subfamily B
Adjusted Odds Ratio
Cytochrome P450 enzyme
Human immunodeficiency virus type 1
Inter quartile range
Roussel Uclaf Causality Assessment Methods
Single nucleotide polymorphisms
Milinkovic A, Martínez E. Nevirapine in the treatment of HIV. Expert Rev Anti-Infect Ther. 2004;2:367–73.
De Boissieu P, Dramé M, Raffi F, et al. Dat’AIDS study group. Long-term efficacy and toxicity of abacavir/lamivudine/nevirapine compared to the most prescribed ARV regimens before 2013 in a French Nationwide cohort study. Medicine (Baltimore). 2016;95(37):e4890.
Reliquet V, Allavena C, Morineau-Le Houssine P, Mounoury O, Raffi F. Twelve-year experience of nevirapine use: benefits and convenience for long-term management in a French cohort of HIV-1-infected patients. HIV Clin Trials. 2010;11(2):110–7.
Llibre JM, Bravo I, Ornelas A, et al. Effectiveness of a treatment switch to Nevirapine plus Tenofovir and Emtricitabine (or lamivudine) in adults with HIV-1 suppressed viremia. PLoS One. 2015;10(6):e0128131.
Patel SM, Johnson S, Belknap SM, et al. Serious adverse cutaneous and hepatic toxicities associated with nevirapine use by non-HIV-infected individuals. J Acquir Immune Defic Syndr. 2004 Feb 1;35(2):120–5.
Stern JO, Robinson PA, Love J, et al. A comprehensive hepatic safety analysis of nevirapine in different populations of HIV infected patients. J Acquir Immune Defic Syndr. 2003;34(Suppl. 1):S21–33.
Dieterich DT, Robinson PA, Love J, Stern JO. Drug-induced liver injury associated with the use of nonnucleoside reverse-transcriptase inhibitors. Clin Infect Dis. 2004;38(Suppl. 2):S80–9.
Wu PY, Cheng CY, Liu CE, et al. Multicenter study of skin rashes and hepatotoxicity in antiretroviral-naïve HIV-positive patients receiving non-nucleoside reverse-transcriptase inhibitor plus nucleoside reverse-transcriptase inhibitors in Taiwan. PLoS One. 2017 Feb 21;12(2):e0171596.
Kesselring AM, Wit FW, Sabin CA, Lundgren JD, Gill MJ, Gatell JM, et al. Risk factors for treatment-limiting toxicities in patients starting nevirapine-containing antiretroviral therapy. AIDS. 2009;23:1689–99.
Mocroft A, Staszewski S, Weber R, et al. EuroSIDA study group. Risk of discontinuation of nevirapine due to toxicities in antiretroviral-naïve and -experienced HIV-infected patients with high and low CD4+ T-cell counts. Antivir Ther. 2007;12(3):325–33.
Van Welzen B, Mudrikova T, Arends J, Hoepelman A. No increased risk of hepatotoxicity in long-term use of nonnucleoside reverse transcriptase inhibitors in HIV-infected patients. HIV Med. 2012;13(7):448–52.
Sanne I, Mommeja-Marin H, Hinkle J, et al. Severe hepatotoxicity associated with nevirapine use in HIV-infected subjects. J Infect Dis. 2005;191(6):825–9.
Riska P, Lamson M, MacGregor T, Sabo J, Hattox S, Pav J, Keirns J. Disposition and biotransformation of the antiretroviral drug nevirapine in humans. Drug Metab Dispos. 1999;27:895–901.
Bertrand J, Chou M, Richardson DM, et al. ANRS 12154 study group. Multiple genetic variants predict steady-state nevirapine clearance in HIV-infected Cambodians. Pharmacogenet Genomics. 2012;22(12):868–76.
Penzak SR, Kabuye G, Mugyenyi P, et al. Cytochrome P450 2B6 (CYP2B6) G516T influences nevirapine plasma concentrations in HIV-infected patients in Uganda. HIV Med. 2007;8(2):86–91.
Wyen C, Hendra H, Vogel M, et al. Impact of CYP2B6 983T>C polymorphism on non-nucleoside reverse transcriptase inhibitor plasma concentrations in HIV-infected patients. J Antimicrob Chemother. 2008;61:914–8.
Rotger M, Colombo S, Furrer H, et al. Influence of CYP2B6 polymorphism on plasma and intracellular concentrations and toxicity of efavirenz and nevirapine in HIV-infected patients. Pharmacogenet Genomics. 2005;15:1–5.
Gozalo C, Grard L, Loiseau P, et al. For the ANRS 081 study group pharmacogenetics of toxicity, plasma trough concentration and treatment outcome with Nevirapine-containing regimen in anti-retroviral-nave HIV-infected adults: an exploratory study of the TRIANON ANRS 081 trial. Basic & Clinical Pharmacology & Toxicology. 2011;109:513–20.
Haas DW, Bartlett JA, Andersen JW, et al. Pharmacogenetics of nevirapine-associated hepatotoxicity: an adult AIDS Clinical Trials Group collaboration. Clin Infect Dis. 2006;43:783–6.
Ciccacci C, Borgiani P, Ceffa S, et al. Nevirapine-induced hepatotoxicity and pharmacogenetics: a retrospective study in a population from Mozambique. Pharmacogenomics. 2010;11:23–31.
Danan G, Teschke R. RUCAM in drug and herb induced liver injury: the update. Int J Mol Sci. 2016;17(1):14.
Brück S, Witte S, Brust J, et al. Hepatotoxicity in patients prescribed efavirenz or nevirapine. Eur J Med Res. 2008;13:343–8.
Bekker Z, Walubo A, du Plessis JB. The role of the immune system in Nevirapine-induced subclinical liver injury of a rat model. ISRN Pharmaceutics. 2012;2012:932542.
Núñez M. Hepatotoxicity of antiretrovirals: incidence, mechanisms and management. J Hepatol. 2006;44(1 Suppl):S132–9 Epub 2005 Nov 28.
Sulkowski MS, Thomas DL, Chaisson RE, Moore RD. Hepatotoxicity associated with antiretroviral therapy in adults infected with human immunodeficiency virus and the role of hepatitis C or B virus infection. JAMA. 2000;283:74–80.
Vogel M, Bertram N, Wasmuth JC, Wyen C, Voigt E, Schwarze-Zander C, Sudhop T, Fätkenheuer G, Rockstroh JK, Reichel C. Nevirapine pharmacokinetics in HIV-infected and HIV/HCV-coinfected individuals. J Antimicrob Chemother. 2009 May;63(5):988–91.
Wen B, Chen Y, Fitch WL. Metabolic activation of nevirapine in human liver microsomes: dehydrogenation and inactivation of cytochrome P450 3A4. Drug Metab Dispos. 2009 Jul;37(7):1557–62.
Schipani A, Wyen C, Mahungu T, et al. German competence network for HIV/AIDS. Integration of population pharmacokinetics and pharmacogenetics: an aid to optimal nevirapine dose selection in HIV-infected individuals. J Antimicrob Chemother. 2011;66(6):1332–9.
Mahungu T, Smith C, Turner F, et al. Cytochrome P450 2B6 516G-->T is associated with plasma concentrations of nevirapine at both 200 mg twice daily and 400 mg once daily in an ethnically diverse population. HIV Med. 2009 May;10(5):310–7.
Giacomelli A, Rusconi S, Falvella FS, et al. Clinical and genetic determinants of nevirapine plasma trough concentration. SAGE Open Med. 2018;6:2050312118780861.
Yuan J, Guo S, Hall D, et al. Nevirapine Toxicogenomics study team. Toxicogenomics of nevirapine-associated cutaneous and hepatic adverse events among populations of African, Asian, and European descent. AIDS. 2011;25(10):1271–80.
Owen A, Pirmohamed M, Khoo HS, Back DJ. Pharmacogenetics of HIV therapy. Pharmacogenet Genomics. 2006;16:693–703.
Zanger UM, Schwab M. Cytochrome P450 enzymes in drug metabolism: regulation of gene expression, enzyme activities, and impact of genetic variation. Pharmacol Ther. 2013 Apr;138(1):103–41.
Kitzmiller JP, Sullivan DM, Phelps MA, Wang D, Sadee W. CYP3A4/5 combined genotype analysis for predicting statin dose requirement for optimal lipid control. Drug Metabol Drug Interact. 2013;28(1):59–63.
We thank Mrs. Tiziana Formenti for the excellent technical help.
No financial support to this study.
Availability of data and materials
The data sets used and/or analyzed during the current study available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The study was reviewed and approved by our ethics committee (Comitato Etico Interaziendale, Milano area 1); all subjects signed a dedicated informed consent.
Consent for publication
S.R. has received consultancy payments and speaking fee from Bristol-Myers Squibb, Gilead, ViiV Healthcare, Merck Sharp Dohme, ABBvie and Janssen. M.G. has received consultancy payments and speaking fee from Bristol-Myers Squibb, Gilead, ViiV Healthcare, Merck Sharp Dohme, ABBvie, Janssen and Roche.
Preliminary data of this study were presented as poster presentation (PE10/15) at the 16th European AIDS conference, October 25–27, 2017 Milan, Italy.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.