Figures
Abstract
Adverse reactions are the main cause of treatment discontinuation among HIV+ individuals. Genes related to drug absorption, distribution, metabolism and excretion (ADME) influence drug bioavailability and treatment response. We have investigated the association between single nucleotide polymorphisms (SNPs) in 29 ADME genes and intolerance to therapy in a case-control study including 764 individuals. Results showed that 15 SNPs were associated with intolerance to nucleoside and 11 to non-nucleoside reverse transcriptase inhibitors (NRTIs and NNRTIs), and 8 to protease inhibitors (PIs) containing regimens under alpha = 0.05. After Bonferroni adjustment, two associations remained statistically significant. SNP rs2712816, at SLCO2B1 was associated to intolerance to NRTIs (ORGA/AA = 2.37; p = 0.0001), while rs4148396, at ABCC2, conferred risk of intolerance to PIs containing regimens (ORCT/TT = 2.64; p = 0.00009). Accordingly, haplotypes carrying rs2712816A and rs4148396T alleles were also associated to risk of intolerance to NRTIs and PIs, respectively. Our data reinforce the role of drug transporters in response to HIV therapy and may contribute to a future development of personalized therapies.
Citation: Arruda MB, Campagnari F, de Almeida TB, Couto-Fernandez JC, Tanuri A, Cardoso CC (2016) Single Nucleotide Polymorphisms in Cellular Drug Transporters Are Associated with Intolerance to Antiretroviral Therapy in Brazilian HIV-1 Positive Individuals. PLoS ONE 11(9): e0163170. https://doi.org/10.1371/journal.pone.0163170
Editor: Sanjay B. Maggirwar, University of Rochester Medical Center, UNITED STATES
Received: July 11, 2016; Accepted: September 2, 2016; Published: September 20, 2016
Copyright: © 2016 Arruda et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was supported by the “Brazilian Ministry of Science and Technology” [FINEP 01/2010 – contract number 03.12.0057.00]. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Deoxi Biotecnologia provided support in the form of salaries for authors [FC], but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The specific roles of this author are articulated in the ‘author contributions’ section.
Competing interests: Francine Campagnari is employed by Deoxi Biotecnologia. There are no patents, products in development or marketed products to declare. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials. FC has worked in the present study as part of her PhD project.
Introduction
The infection by HIV-1 is a pandemic condition that affects nearly 37 million people worldwide [1]. The treatment of HIV positive individuals is based on the combined antiretroviral therapy (cART), which consists on the administration of three antiretrovirals (ARVs) from at least two different classes. These classes include the nucleoside reverse transcriptase inhibitors (NRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTIs), protease inhibitors (PIs), fusion and entry inhibitors, and integrase inhibitors. Introduction of HAART in 1996 has led to clear health improvements to the HIV+ individuals, with reduced morbidity and increase in life expectancy [2]. However, the emergence of drug resistance and the development of adverse reactions to the ARVs are still of major concern since they may impair therapy effectiveness.
Adverse reactions (ARs) are one of the most important factors associated with reduced life quality among individuals undergoing antiretroviral therapy. In fact, ARs have been reported as the most common cause of treatment discontinuation and changes in cART regimen [3–6]. Gastrointestinal intolerance, anemia and hypersensitivity reactions are commonly observed since the early beginning of the treatment. The long-term use of the ARVs may also lead to metabolic disorders, including lipodystrophy and dyslipidemia, and peripheral neuropathies [5–7]. Moreover, ARs can also limit treatment effectiveness by reducing adherence [8–11] and, consequently, favoring the emergence of drug resistant viruses.
Variations in genes coding for molecules involved in drug absorption, distribution, metabolism and excretion (ADME) have been consistently associated to both short-term and long-term adverse reactions to the ARVs. The impact of such variations is clearly reflected by the interindividual variations in antiretroviral plasma concentrations after taking the same dosage. Therefore, tolerance to treatment can also be substantially variable. Although much remains unknown, there is extensive data regarding the influence of ABCB1 variations and plasma or intracellular levels of protease inhibitors [12]. Similar results have been observed for a second ABC transporter, coded by ABCC2. In addition, single nucleotide polymorphisms (SNPs) and haplotypes of CYP2B6 have been consistently associated to efavirenz clearance and also to central nervous system adverse reactions to this drug [13–15]. Moreover, atazanavir-associated hyperbilirubinemia was associated to a microsatellite at UGT1A1 promoter [16]. Mitochondrial DNA haplogroups and variations in genes coding for inflammatory mediators and apolipoproteins have also been suggested as predictors for lipid disorders [17–19].
According to recent estimates, there are about 781,000 people living with HIV in Brazil, and 405,000 are currently undergoing treatment [20]. However, just a few studies have been conducted to describe the role of host genetics in response to cART in our population [16, 21–25]. Therefore, the prevalence of SNPs in ADME genes among Brazilian HIV positive individuals remains largely unknown.
The present study was designed to investigate the association between 346 SNPs in 29 ADME genes and intolerance to antiretroviral therapy among Brazilians. For this purpose, 764 individuals undergoing antiretroviral therapy were enrolled in a case-control study. Our data showed a clear association between ABCC2 SNPs and treatment modification due to intolerance to protease inhibitors, while SLCO2B1 genetic variations increased the risk of intolerance to NRTIs.
Methods
Ethics Statement
The present study was approved by Fiocruz Institutional Review Board (IRB) and a written informed consent was obtained from all subjects. All analyses were conducted according to the principles expressed in the Declaration of Helsinki.
Subjects and study design.
A total of 764 HIV-1+ individuals were enrolled in this study. All samples were obtained from the biorepository of the AIDS and Molecular Immunology Laboratory at Fiocruz, Rio de Janeiro, Brazil. This laboratory is a part of the Brazilian Network for HIV-1 Genotyping (RENAGENO), which was implemented in 2001 with the aim of providing free HIV genotyping tests for patients in HAART failure [26]. Patients from both genders, born in Rio de Janeiro and with a minimum age of 18 years old were considered eligible. A patient was included in the study if he/she had: (1) a complete history of cART schemes since first line therapy available in the file, (2) buffy coat samples available for DNA extraction. Patients under treatment for hepatitis or tuberculosis and women who exhibited adverse reactions during pregnancy were excluded from the study. The patients who developed adverse reactions to any antiretroviral were classified as cases (N = 359), while the control group (N = 405) included patients to whom the treatment was safe. Considering a minor allele frequency of 0.05 and alpha = 0.0001 (Bonferroni adjustment for 346 candidate SNPs), the minimum OR value to reach 80% of power with this sample size was 2 under an additive model.
A subject was not included in the control group if he/she has not been treated for at least 6 months. A total of 195 controls (48%) were still undergoing first-line therapy (first scheme) at the moment of data collection, with a medium time of 3.9 ± 2 years using the same drug combination. For those who had already changed treatment, the medium time of permanence in each scheme was 2.87 ± 1.87 years.
Clinical description of specific adverse reactions was available for a limited number of subjects (N = 65 for cases of intolerance to nucleoside reverse transcriptase inhibitors, 26 for non-nucleoside reverse transcriptase inhibitors, and 32 for protease inhibitors). Therefore, the outcome “intolerance” was used for case definition. Among the cases of intolerance to NRTIs, 31 developed anemia attributed to zidovudine and 10 exhibited lipodistrophy due to stavudine, which was also associated to peripheral neuropathy in 9 individuals. Fifteen individuals developed dyslipidemia and the remaining 8 cases exhibited rash, gastrointestinal adverse effects and lactic acidosis (N = 3, 3 and 2, respectively). Among the subjects who developed intolerance to NNRTIs, the most prevalent adverse reaction was cutaneous rash after nevirapine use (N = 16), followed by neuropsychiatric disorders in response to efavirenz (N = 7) and gastrointestinal disorders (N = 3). In the group of intolerance to protease inhibitors, gastrointestinal adverse reactions including diarrhea and nausea due to lopinavir/ritonavir use were the most prevalent (N = 11), followed by nephrolithiasis due to indinavir (N = 9). Lipodistrophy and dyslipidaemia were observed in 5 and 3 cases, respectively. The remaining 4 cases exhibited jaundice due to atazanavir and cutaneous rash (N = 2, each).
First of all, an overall analysis was developed including all cases and controls (“All ARVs”), regardless the drug class. Afterwards, three different case-control studies were conducted to characterize association between pharmacogenetic markers and intolerance to: (1) nucleoside reverse transcriptase inhibitors (NRTIs), (2) non- nucleoside reverse transcriptase inhibitors (NNRTIs) or (3) protease inhibitors (PIs). For this analysis, all individuals exhibiting adverse reactions to drugs belonging to the same class were grouped in a single case-control.
Genetic ancestry was determined using a panel of 28 ancestry informative SNPs previously validated for Brazilian population [27]. Distribution of each gender and proportions of European and African ancestries are represented in Table 1.
Single Nucleotide Polymorphisms (SNPs) Selection
Selection of candidate genes and SNPs was based mainly on literature search, as a way to replicate previous associations. We have also searched for additional SNPs in coding and regulatory regions of all candidate genes using the SNPper tool. In order to increase gene coverage, tag SNPs were selected from HapMap data bank using the following parameters: minor allele frequency of 0.05 in CEU (Utah residents with Northern and Western European ancestry) or YRI (Yoruba in Ibadan, Nigeria) populations and r2 cutoff of 0.8. Following this strategy, 346 SNPs in 29 candidate genes were selected (S1 Table).
DNA Extraction and SNP Genotyping
Genomic DNA was obtained from buffy coat samples using the QIAamp® DNA Blood Mini Kit (QIAGEN) according to manufacturer’s instructions. A total of 345 SNPs was genotyped using Illumina Golden Gate® assays (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. Briefly, the assay includes two allele-specific oligos (ASO) and a locus-specific oligo (LSO) that hybridizes several bases downstream from the SNP site. The three oligonucleotide sequences contain regions of complementarity with universal PCR primer sites. DNA samples (250 ng) were activated and hybridized to the assay oligonucleotides designed for each SNP locus. Following hybridization, several wash steps were performed before base extension and ligation of the ASO and LSO primers. Amplification of the full length products was performed using Cy3- and Cy5-labeled universal PCR primers and Titanium® DNA polymerase (CloneTech, Mountain View, CA, USA). Finally, Golden Gate® assay products were hybridized to the BeadChip and the fluorescence was analyzed using an iScan instrument (Illumina, San Diego, CA, USA). Genotype calling was performed using the software Genome Studio (Illumina, San Diego, CA, USA). The SNP rs3745274, located at CYP2B6 gene, did not pass Illumina criteria for assay design and, therefore, was genotyped using SNaPshot® analysis (Thermo Scientific, Massachusetts, USA). Briefly, PCR amplification was performed using QIAGEN Multiplex PCR kit (QIAGEN, Hilden, Germany) and specific forward (5’-GATTGAACACCTACTCTGCCCAGCC-3’) and reverse (5’-AGACGATGGAGCAGATGATGTTGGC-3) primers at 0.2 μM. PCR products were cleaned using 3U of FastAP (Thermo Scientific, Massachusetts, USA) and 1U of Exonucelase I (Affimetrix, California, USA). SNaPshot reactions were performed using the SNaPshot™ Multiplex kit (Thermo Scientific, Massachusetts, USA) according to manufacturer’s instructions and the specific primer for rs3745274 (5’-CCCTCATGGACCCCACCTTCCTCTTCCA-3’) at 2 μM. SNaPshot products were treated with FastAP 1U and submitted to capillary electrophoresis on ABI3130 Genetic Analyzer (Thermo Scientific, Massachusetts, USA) using the standard fragment analysis protocol. Genotyping was performed using GeneMapper software, version 4.0 (Thermo Scientific, Massachusetts, USA)”.
Statistical Analyses
Deviations from Hardy-Weiberg Equilibrium (HWE) were assessed by χ2 tests. Comparisons between cases and controls were performed using unconditional logistic regression models controlling for potential confounders such as gender, age at diagnosis and proportions of European and African ancestries. A stepwise logistic regression analysis was performed to select covariates. All SNPs were analyzed under dominant and codominant models. A Bonferroni adjustment was applied to avoid type I error. A Cochran–Armitage trend test was conducted to describe possible allele dose effects. Pairwise linkage disequilibrium (LD) patterns were determined using the r2 statistics and a cutoff of r2 ≥ 0.8 to define “tags”. SNPs with nominal p-values < 0.05 in both models (dominant and codominant) were included in haplotype analyses. Haplotype frequencies were estimated by maximum likelihood and compared using the same logistic regression models applied for isolated SNPs. All analyses were performed using the R for windows, version 2.14.1, with the packages ‘‘genetics”, “gap”, “SNPassoc”, ‘‘haplo.stats” and ‘‘coin”.
Results
SNPs in ADME Genes Are Associated with Intolerance to Different cART Regimens
A total of 346 SNPs were genotyped in our study. After quality control, 76 SNPs were excluded from analysis due to HWE deviations (N = 38) or because they were monomorfic (N = 38) in our control sample. The minor allele frequencies are shown for each SNP in S1 Table. The remaining 270 candidate SNPs were analyzed separately using logistic regression models adjusted for sex and genetic ancestry. First of all, patients were compared according to the development of adverse reactions to any antiretroviral (All ARVs), and a total of 12 SNPs were significantly associated to this outcome in both multivariate models (dominant and codominant adjusted for gender and genetic ancestry) under a 0.05 significance level (S2 Table). When statistical analyses were stratified according to the ARVs classes, 15 SNPs were associated to adverse reactions to NRTIs, 11 to NNRTIs and 8 to PIs (p < 0.05 in dominant and codominant models; S2 Table). After Bonferroni adjustment for multiple comparisons (270 SNPs; alpha = 0.0002), two associations remained statistically significant (S2 Table, Tables 2 and 3). The SNP rs2712816 (G>A) at SLCO2B1 gene was associated to increased risk of adverse reactions to NRTIs (OR = 2.37; p = 0.0001 for GA/AA genotypes after adjustment for gender and genetic ancestry; S2 Table), while rs4148396 (C>T), located at ABCC2, conferred risk of intolerance to PIs (OR = 2.64; p = 0.00009 for CT/TT genotypes; S2 Table). OR values remained virtually the same when patients who developed adverse reactions to both NRTIs and PIs (N = 25) were removed from analysis (data not shown).
Association between SLCO2B1 and Intolerance to Nucleoside Reverse Transcriptase Inhibitors
In addition to rs2712816, three SNPs at SLCO2B1 (rs12422149, rs1676885 and rs949069) were associated to adverse reactions to NRTIs at a 0.05 significance level (S2 Table). SNP rs1676885 was associated to increased risk of intolerance (adjusted OR = 1.55 for AG/GG), while results obtained for rs949069 (adjusted OR = 0.57 for GA/AA) and rs12422149 (adjusted OR = 0.5 for GA/AA) suggested a protective effect. In all cases, the results obtained suggested an allele-dose effect, with more prominent effects observed among homozygotes for the minor allele (Table 2). This trend was confirmed by a Cochran-Armitage test, which resulted in adjusted p-values of 0.0005 for rs2712816, 0.003 for rs949069 and 0.01 for rs12422149 and rs1676885.
Results of linkage disequilibrium (LD) analyses have showed low association between the four markers tested (r2 < 0.8). Therefore, all SNPs were included in haplotype analysis. Comparisons of haplotype frequencies in cases and controls corroborated the data obtained when SNPs were analyzed separately (Table 2). Indeed, both haplotypes carrying rs2712816A allele increased the susceptibility to adverse reactions to NRTIs when compared to the baseline rs2712816G/rs949069G /rs12422149G /rs1676885A (adjusted OR = 1.39 for A/G/G/A and 1.49 for A/G/G/G; p = 0.04 for both haplotypes).
Association between ABCC2 and Intolerance to Protease Inhibitors
Besides rs4148396, four other SNPs at ABCC2 (rs2073337, rs2804398, rs2804400, rs3740066) were associated with intolerance to PIs (p < 0.05; S2 Table, Table 3). From these, rs2804398 and rs2804400 were in almost perfect linkage disequilibrium with rs4148396 (r2 > 0.9) while rs2073337 was in moderate/strong LD (r2 = 0.66). Therefore, rs4148396 was defined as a tag for the entire bin. SNP rs3740066 was also included in the following analyses, since the low LD patterns with rs4148396 (r2 = 0.37) indicate that this variation belongs to an independent bin. The results obtained from both SNPs indicated an association with increased risk of intolerance to PIs containing regimens (adjusted OR = 2.64 for rs4148396 T carriers and 1.82 for rs3740066 A carriers; S2 Table). Moreover, the OR values suggested that the homozygotes were at higher risk of developing intolerance to PIs as compared to the heterozygotes (3.68 vs 2.47 for rs4148396 and 2.27 vs 1.74 for rs3740066; Table 3), suggesting an allele dose effect. This trend was confirmed by a Cochran-Armitage trend test (p = 0.0002 for rs4148396 and 0.03 for rs3740066).
Next, rs4148396 and rs3740066 alleles were combined in haplotypes. After analysis, the combination of rs4148396T/rs3740066A was associated to increased risk of adverse reactions to PIs (adjusted OR = 2.09; p = 0.0006), as compared to the baseline rs4148396C/rs3740066G, reinforcing the risk effect observed for rs4148396T and rs3740066A in single SNP analysis (Table 3).
Discussion
In this study, we have investigated the association between SNPs in ADME genes and intolerance to cART regimens. Our results showed a consistent association between SLCO2B1 gene variations and intolerance to NRTIs, while ABCC2 polymorphisms conferred risk of intolerance to protease inhibitors containing regimens.
To date, there is no data regarding NRTIs transport through OATP transporters. In fact, OATP1B1, coded by SLCO1B1, plays a role in protease inhibitors uptake [28] and polymorphisms in this gene are associated to variations in plasma concentrations of these drugs [23]. Indeed, the polymorphism rs10444413 was associated with intolerance to PIs in our cohort considering a 0.05 significance level (S2 Table).
The organic anion transporter polypeptide OATP2B1, coded by SLCO2B1, is expressed in hepatocytes [29] and apical membrane of intestinal epithelial cells [30], among other tissues. PIs were also shown to inhibit OATP2B1 activity, suggesting that they might be a target for this transporter [31]. Several reports have shown a clear role for OATP2B1 in absorption of statins and steroid hormones, especially androgens, as well as its precursors. Accordingly, SNPs at SLCO2B1 may predict response to these compounds [32,33]. The missense variation rs12422149 has been consistently associated to worse prognosis for prostate cancer after androgen-deprivation therapy [33]. Differential patterns of hormones absorption might indirectly influence susceptibility to lipid disorders, which are commonly observed in response to ARVs such as stavudine. In addition, since statins are commonly used in the management of lipid disorders associated to HIV infection and ARVs use [34], SNPs at SLCO2B1 may also limit the effectiveness of this treatment in individuals undergoing antiretroviral therapy. Further analyses are required to test these hypotheses and define the role of OATP2B1 in response to antiretroviral therapy.
Beyond the classical candidate ABCB1 (P-glycoprotein), the ABCC transporters, or multidrug resistance-associated proteins (MRP), also promote the efflux of many drugs from cells. The MRP2 transporter, encoded by ABCC2, has been consistently associated to the efflux of protease inhibitors from hepatocytes and peripheral blood mononuclear cells [35–38]. In the present study, the most prominent effect was observed for the SNP rs4148396, which was associated with intolerance to PIs even after Bonferroni adjustment (ORCT/TT = 2.64; p = 0.00009; Table 3). This SNP, as well as rs3740066 were previously associated to neurotoxicity after treatment with the adjuvant FOLFOX4 in patients with colorectal cancer [39]. Besides ABCC2, SNPs in other two ABC transporters were also associated to intolerance to PIs (ABCC4 and ABCG2, S2 Table), reinforcing the role of efflux transporters in adverse reactions to this class of drugs [12].
SNPs at ABCC2 have been consistently associated to tenofovir-induced tubulopathy [40,41]. Moreover, this transporter is also inhibited by NRTIs and NNRTIs [42], which may increase the risk of drug-drug interactions in patients undergoing cART. Notably, SNPs in this gene were also associated to intolerance in comparisons including all antiretrovirals at a 0.05 significance level (S2 Table).
We were unable to detect classical associations such as CYP2B6 and response to efavirenz and ABCB1 SNPs and response to protease inhibitors. A limitation of our study was the inability to perform case-control analyses specific for each drug. In this case, for example, the inclusion of individuals with intolerance to nevirapine may have diluted the effect of CYP2B6 SNPs, which is more prominent for efavirenz response. Another possibility is the outcome under investigation, which was intolerance instead of a clearly defined adverse reaction. As a consequence, we would have limited power to detect subtle or specific effects. More specifically, a reduced number of individuals exhibiting central nervous system reactions to efavirenz among the cases of intolerance to NNRTIs might also explain the lack of association with CYP2B6 gene. Nevertheless, we cannot rule out the possibility that these effects are not valid in our population. The effect of CYP2B6 and ABCB1 SNPs are dependent on genetic background [43]. In this study, the genetic background was determined using a panel of ancestry informative SNPs previously validated for the Brazilian population [27], and no associations were found even after adjustment for genetic ancestry.
In Brazil, the treatment is provided to all HIV+ individuals free of charge and a National Network for HIV Genotyping was implemented to detect HIV resistance and avoid virologic failure [26]. However, just a few studies were conducted to determine the role of host genetics in treatment response. The number of people undergoing cART is expected to increase in the next few years due to new policies that recommend starting treatment at the time of diagnosis to prevent transmission. In addition, the increase in life expectancy of HIV+ individuals will also increase the possibility of drug-drug interactions due to polypharmacy. In this scenario, it is crucial to determine the distribution of ADME markers in our population.
To our knowledge, this is the largest study conducted to describe the role of ADME genes in response to antiretroviral therapy among Brazilians until today, considering the number of candidate genes and polymorphisms. In the future, this information may help to define a panel of markers that may be used along with clinical follow up and help clinicians to manage the different profiles of response to cART using personalized regimens.
Supporting Information
S1 Table. Candidate genes and SNPs selected for case-control analysis.
https://doi.org/10.1371/journal.pone.0163170.s001
(XLSX)
S2 Table. Results of statistically significant associations between candidate SNPs and intolerance to the different antiretroviral classes.
https://doi.org/10.1371/journal.pone.0163170.s002
(DOCX)
Author Contributions
- Conceptualization: MBA FC AT CCC.
- Data curation: MBA CCC.
- Formal analysis: MBA TBA CCC.
- Funding acquisition: FC AT CCC.
- Methodology: MBA FC TBA JCCF CCC.
- Project administration: MBA FC CCC.
- Resources: FC JCCF AT CCC.
- Supervision: MBA AT CCC.
- Visualization: CCC.
- Writing – original draft: MBA TBA CCC.
- Writing – review & editing: FC AT JCCF.
References
- 1.
Joint United Nations Programme on HIV/AIDS (UNAIDS). AIDS by the numbers. Jc2571/1/E [Internet]. 2015;1–11. Available: http://search.unaids.org.
- 2. Mocroft A, Ledergerber B, Katlama C, Kirk O, Reiss P, D’Arminio Monforte A, et al. Decline in the AIDS and death rates in the EuroSIDA study: An observational study. Lancet. 2003;362(9377):22–9. pmid:12853195
- 3. d’Arminio Monforte A, Lepri AC, Rezza G, Pezzotti P, Antinori A, Phillips AN, et al. Insights into the reasons for discontinuation of the first highly active antiretroviral therapy (HAART) regimen in a cohort of antiretroviral naïve patients. I.CO.N.A. Study Group. Italian Cohort of Antiretroviral-Naïve Patients. Aids. 2000;14(5):499–507. pmid:10780712
- 4. Mocroft A, Youle M, Moore A, Sabin CA, Madge S, Lepri AC, et al. Reasons for modification and discontinuation of antiretrovirals: results from a single treatment centre. Aids. 2001;15(2):185–94. pmid:11216926
- 5. Cardoso SW, Grinsztejn B, Velasque L, Veloso VG, Luz PM, Friedman RK, et al. Incidence of modifying or discontinuing first HAART regimen and its determinants in a cohort of HIV-infected patients from Rio de Janeiro, Brazil. AIDS Res Hum Retroviruses. 2010;26(8):865–74. pmid:20672997
- 6. Ribeiro FA, Tupinambás U, Fonseca MO, Greco DB. Durability of the first combined antiretroviral regimen in patients with AIDS at a reference center in Belo Horizonte, Brazil, from 1996 to 2005. Brazilian J Infect Dis. 2012;16(1):27–33.
- 7. Haubrich RH, Riddler SA, Dirienzo AG, Komarow L, Haas DW, Mellors JW, et al. Metabolic Outcomes in a Randomized Trial of Nucleoside, Nonnucleoside and Protease Inhibitor-Sparing Regimens for Initial HIV Treatment. Aids. 2009;23(9):1109–18. pmid:19417580
- 8. Ammassari A, Murri R, Pezzotti P, Trotta MP, Ravasio L, De Longis P, et al. Self-reported symptoms and medication side-effects influence adherence to highly active antiretroviral therary in persons with HIV infection. Journal of Acquired Immune Deficiency Syndrome. 2001. p. 445–9.
- 9. Guimarães MDC, Rocha GM, Campos LN, de Freitas FMT, Gualberto FAS, Teixeira R d’ AR, et al. Difficulties reported by HIV-infected patients using antiretroviral therapy in Brazil. Clinics (Sao Paulo). 2008;63(2):165–72.
- 10. Pádua C a M De, César CC, Bonolo PF, Acurcio F a, Guimarães MDC. Self-reported adverse reactions among patients initiating antiretroviral therapy in Brazil. Brazilian J Infect Dis. 2007;11(1):20–6.
- 11. Bonolo P de F, Cesar CC, Acurcio FA, Ceccato M das GB, de Padua CAM, Alvares J, et al. Non-adherence among patients initiating antiretroviral therapy: a challenge for health professionals in Brazil. AIDS. 2005;19 Suppl 4:S5–13. pmid:16249654
- 12. Michaud V, Bar-magen T, Turgeon J, Flockhart D, Desta Z, Wainberg MA. The Dual Role of Pharmacogenetics in HIV Treatment : Mutations and Polymorphisms Regulating Antiretroviral Drug Resistance and Disposition. 2012;64(3):803–33.
- 13. Haas DW, Ribaudo HJ, Kim RB, Tierney C, Wilkinson GR, Gulick RM, et al. Pharmacogenetics of efavirenz and central nervous system side effects: an Adult AIDS Clinical Trials Group study. Aids. 2004;18(18):2391–400. pmid:15622315
- 14. Haas DW, Smeaton LM, Shafer RW, Robbins GK, Morse GD, Labbe L, et al. Pharmacogenetics of long-term responses to antiretroviral regimens containing Efavirenz and/or Nelfinavir: an Adult Aids Clinical Trials Group Study. J Infect Dis [Internet]. 2005;192(11):1931–42. pmid:16267764
- 15. Ribaudo HJ, Haas DW, Tierney C, Kim RB, Wilkinson GR, Gulick RM, et al. Pharmacogenetics of plasma efavirenz exposure after treatment discontinuation: an Adult AIDS Clinical Trials Group Study. Clin Infect Dis. 2006;42(December 2005):401–7.
- 16. Rodríguez-Nóvoa S, Martín-Carbonero L, Barreiro P, González-Pardo G, Jiménez-Nácher I, González-Lahoz J, et al. Genetic factors influencing atazanavir plasma concentrations and the risk of severe hyperbilirubinemia. AIDS. 2007;21(1):41–6. pmid:17148966
- 17. Tarr PE, Taffé P, Bleiber G, Furrer H, Rotger M, Martinez R, et al. Modeling the influence of APOC3, APOE, and TNF polymorphisms on the risk of antiretroviral therapy-associated lipid disorders. J Infect Dis. 2005;191(9):1419–26. pmid:15809899
- 18. Hulgan T, Tebas P, Canter J a, Mulligan K, Haas DW, Dubé M, et al. Hemochromatosis gene polymorphisms, mitochondrial haplogroups, and peripheral lipoatrophy during antiretroviral therapy. J Infect Dis. 2008;197(6):858–66. pmid:18419350
- 19. Rotger M, Bayard C, Taffé P, Martinez R, Cavassini M, Bernasconi E, et al. Contribution of genome-wide significant single-nucleotide polymorphisms and antiretroviral therapy to dyslipidemia in HIV-infected individuals: A longitudinal study. Circ Cardiovasc Genet. 2009;2(6):621–8. pmid:20031643
- 20.
Ministério da Saúde. Boletim Epidemilógico HIV/AIDS. Bol Epidemiológico HIV/AIDS. 2015;1–100.
- 21. Turatti L, Sprinz E, Lazzaretti RK, Kuhmmer R, Agnes G, Silveira JM, et al. Short communication: UGT1A1*28 variant allele is a predictor of severe hyperbilirubinemia in HIV-infected patients on HAART in southern Brazil. AIDS Res Hum Retroviruses. 2012;28(9):1015–8. pmid:22050734
- 22. Estrela RDC, Ribeiro FS, Barroso PF, Tuyama M, Gregório SP, Dias-Neto E, et al. ABCB1 polymorphisms and the concentrations of lopinavir and ritonavir in blood, semen and saliva of HIV-infected men under antiretroviral therapy. Pharmacogenomics. 2009;10(2):311–8. pmid:19207033
- 23. Kohlrausch FB, De Cássia Estrela R, Barroso PF, Suarez-Kurtz G. The impact of SLCO1B1 polymorphisms on the plasma concentration of lopinavir and ritonavir in HIV-infected men. Br J Clin Pharmacol. 2010;69(1):95–8. pmid:20078617
- 24. Lazzaretti RK, Gasparotto AS, Sassi MGDM, Polanczyk C a., Kuhmmer R, Silveira JM, et al. Genetic markers associated to dyslipidemia in HIV-infected individuals on HAART. Sci World J. 2013;2013.
- 25. Coelho AVC, Silva SPS, de Alencar LC a, Stocco G, Crovella S, Brandão L a C, et al. ABCB1 and ABCC1 variants associated with virological failure of first-line protease inhibitors antiretroviral regimens in Northeast Brazil patients. J Clin Pharmacol. 2013;53(12):1286–93. pmid:23996099
- 26. Couto-Fernandez JC, Silva-de-Jesus C, Veloso VG, Rachid M, Gracie RSG, Chequer-Fernandez SL, et al. Human immunodeficiency virus type 1 (HIV-1) genotyping in Rio de Janeiro, Brazil: assessing subtype and drug-resistance associated mutations in HIV-1 infected individuals failing highly active antiretroviral therapy. Mem Inst Oswaldo Cruz. 2005;100(1):73–8. pmid:15867968
- 27. Lins TC, Vieira RG, Abreu BS, Grattapaglia D, Pereira RW. Genetic composition of Brazilian population samples based on a set of twenty eight ancestry informative SNPs. Am J Hum Biol. 2010;22(2):187–92. pmid:19639555
- 28. Hartkoorn RC, Kwan WS, Shallcross V, Chaikan A, Liptrott N, Egan D, et al. HIV protease inhibitors are substrates for OATP1A2, OATP1B1 and OATP1B3 and lopinavir plasma concentrations are influenced by SLCO1B1 polymorphisms. Pharmacogenet Genomics. 2010;20(2):112–20. pmid:20051929
- 29. Kullak-Ublick G a, Ismair MG, Stieger B, Landmann L, Huber R, Pizzagalli F, et al. Organic anion-transporting polypeptide B (OATP-B) and its functional comparison with three other OATPs of human liver. Gastroenterology. 2001;120(2):525–33. pmid:11159893
- 30. Kobayashi D, Nozawa T, Imai K, Nezu J, Tsuji A, Tamai I. Involvement of human organic anion transporting polypeptide OATP-B (SLC21A9) in pH-dependent transport across intestinal apical membrane. J Pharmacol Exp Ther. 2003;306(2):703–8. pmid:12724351
- 31. Kis O, Zastre J a, Ramaswamy M, Bendayan R. pH Dependence of Organic Anion-Transporting Polypeptide 2B1 in Caco-2 Cells : Potential Role in Antiretroviral Drug Oral Bioavailability and Drug–Drug Interactions. J Pharmacol Exp Ther. 2010;334(3):1009–22. pmid:20507927
- 32. Tsamandouras N, Dickinson G, Guo Y, Hall S, Rostami-Hodjegan A, Galetin A, et al. Identification of the effect of multiple polymorphisms on the pharmacokinetics of simvastatin and simvastatin acid using a population-modeling approach. Clin Pharmacol Ther. 2014;96(1):90–100. pmid:24598718
- 33. Yang M, Xie W, Mostaghel E, Nakabayashi M, Werner L, Sun T, et al. SLCO2B1 and SLCO1B3 may determine time to progression for patients receiving androgen deprivation therapy for prostate cancer. J Clin Oncol. 2011;29(18):2565–73. pmid:21606417
- 34. Saves M, Raffi F, Capeau J, Rozenbaum W, Ragnaud J-M, Perronne C, et al. Factors related to lipodystrophy and metabolic alterations in patients with human immunodeficiency virus infection receiving highly active antiretroviral therapy. Clin Infect Dis. 2002;34(10):1396–405. pmid:11981737
- 35. Williams GC, Liu A, Knipp G, Sinko PJ. Direct evidence that saquinavir is transported by multidrug resistance-associated protein (MRP1) and canalicular multispecific organic anion transporter (MRP2). Antimicrob Agents Chemother. 2002;46(11):3456–62. pmid:12384350
- 36. Janneh O, Owen A, Chandler B, Hartkoorn RC, Hart CA, Bray PG, et al. Modulation of the intracellular accumulation of saquinavir in peripheral blood mononuclear cells by inhibitors of MRP1, MRP2, P-gp and BCRP. Aids. 2005;19(18):2097–102. pmid:16284458
- 37. Agarwal S, Pal D, Mitra AK. Both P-gp and MRP2 mediate transport of Lopinavir, a protease inhibitor. Int J Pharm. 2007;339(1–2):139–47. pmid:17451894
- 38. Lubomirov R, di Iulio J, Fayet A, Colombo S, Martinez R, Marzolini C, et al. ADME pharmacogenetics: investigation of the pharmacokinetics of the antiretroviral agent lopinavir coformulated with ritonavir. Pharmacogenet Genomics. 2010;20(4):217–30. pmid:20139798
- 39. Cecchin E, D’Andrea M, Lonardi S, Zanusso C, Pella N, Errante D, et al. A prospective validation pharmacogenomic study in the adjuvant setting of colorectal cancer patients treated with the 5-fluorouracil/leucovorin/oxaliplatin (FOLFOX4) regimen. Pharmacogenomics J. 2013;13(5):403–9. pmid:22868256
- 40. Izzedine H, Hulot J-S, Villard E, Goyenvalle C, Dominguez S, Ghosn J, et al. Association between ABCC2 gene haplotypes and tenofovir-induced proximal tubulopathy. J Infect Dis. 2006;194(11):1481–91. pmid:17083032
- 41. Rodríguez-Nóvoa S, Labarga P, Soriano V, Egan D, Albalater M, Morello J, et al. Predictors of kidney tubular dysfunction in HIV-infected patients treated with tenofovir: a pharmacogenetic study. Clin Infect Dis. 2009;48(11):e108–16. pmid:19400747
- 42. Weiss J, Herzog M, König S, Storch CH, Ketabi-Kiyanvash N, Haefeli WE. Induction of multiple drug transporters by efavirenz. J Pharmacol Sci. 2009;109(2):242–50. pmid:19234366
- 43. Ribaudo HJ, Liu H, Schwab M, Schaeffeler E, Motsinger-reif AA, Ritchie MD, et al. Impact of CYP2B6, ABCB1 and CYP3A5 Polymorphisms on Efavirenz Pharmacokinetics and Treatment Response: An AIDS Clinical Trials Group Study. 2010;202(5):717–22.