Antiretroviral-Drug Resistant Mutations and Genetic Diversity in HIV-1 Infected Individuals in Nigeria ()
1. Introduction
HIV-1 strains can be divided into three distinct and highly divergent groups: M (major), O (outlier), and N (nonM/O). Recently, a new HIV-1 group closely related to gorilla simian immunodeficiency virus (SIVgor) was identified in a Cameroonian which showed no evidence of recombination with other HIV-1 lineages and was designated as HIV-1 group P [1] . Several genetic variants had been recognized within HIV-1 group M, including nine subtypes (A, B, C, D, F, G, H, J, and K) and many CRFs and URFs [2] with differing geographic distribution. Subtype B is the predominant subtype in the western world (Western Europe, the Americas, Japan and Australia) and is responsible for only 12% of global infections, whilst non-B subtypes predominate in the rest of the world [3] . Antiretroviral drugs were developed based on subtype-B isolates and had successfully reduced the incidence of opportunistic infections and HIV related morbidity and mortality [4] . Despite its remarkable achievement, drug resistance is the major obstacle in the management of HIV patient [5] . Drug resistant HIV causes therapeutic problems in treated and also therapy-naïve patients [6] . Emerging data has shown that viral subtype may likewise influence the effectiveness of antiretroviral treatment and that pre-existing mutations could reduce the effectiveness of ARV [7] [8] . Reduced susceptibility to ARVs has been reported in non-B subtypes than in subtype B, and that genetic differences among subtypes could result in differential patterns of resistant mutations in response to ARV pressure [9] . Meanwhile, most data on HIV-1 drug resistance mechanisms are from subtype-B viruses due to availability of genotypic and phenotypic antiretroviral drug resistance testing in the western world.
Nigeria is the most populous country in Africa with a population of 174 million in 2013 report [10] and HIV prevalence of 4.1% in 2010 antenatal survey [11] indicating 3.1 million people are living with HIV/AIDS in Nigeria. HIV 1 and 2 circulate in Nigeria [12] . In the light of expansion to increasing the number on antiretroviral therapy (ART) and creation of universal access, a total of 1,512,720 persons were estimated to require ART, of which 359,181 have been successfully placed on treatment using WHO public health approach to ART [13] . Studies on genetic diversity in Nigeria have revealed subtypes A, G, CRF02_AG, C, J, F2 and Group O among the population [14] -[16] . Apart from common subtypes, other rare complex recombinant forms like CRF06_cpx, CRF01_AE, CRF11_cpx, and A3 have also been identified [16] . Few studies have been done on HIV drug resistance and variants characterization in Nigeria because no easily applicable tool is available to carry out this study in a large population. Earlier data showed a 10% to 17% baseline resistance to major classes of ARVs of viruses in ART naïve Nigerians [16] [17] . One of the studies revealed a significant degree of primary and secondary drug resistance mutations in the RT and PR; as well as polymorphisms at positions of previously characterized drug resistance mutations in ART-naïve HIV-infected subjects from Ibadan and Saki in Oyo State [17] .
Although various studies have reported existence and types of resistance mutations in Nigeria, there is still need for continuous evaluation of both the magnitude and patterns of drug resistance especially with ARV scaleup. This will aid in the development of appropriate prevention strategies to limit wide-scale treatment failure as well as plan for appropriate treatment strategies. This will also provide information regarding the current state of HIV diversity in this region. In this study, we evaluated antiretroviral drug resistant among HIV-infected ARTtreated and naïve patients and further assessed the circulating HIV-1 variants in Nigeria.
2. Methodology
2.1. Study Site
A cross sectional study was conducted on 202 HIV infected patients attending ART clinics for treatment from 2008 to 2010 from all the six geopolitical zones of Nigeria (Table 1).
Table 1. Patient status and baseline characteristics.
ART: Antiretroviral therapy; non-ART: No Antiretroviral therapy. (No treatment or drugnaïve); IQR: Interquartile range; *missing—Data on the patient origin were not available on the questionnaire.
2.2. Study Population and Sample Collection
Study population included both ART naive and experienced patients (≥3 months). Patients who were sick at the time of sample collection and those with incomplete questionnaire were excluded from the study.
Baseline characteristic collected includes age, sex, place of birth, state of origin, ethnicity and mode of infection if known, ARV drug, date of commencement. Blood was drawn for baseline immunological and virological testing at sites. Selection for genotyping assay was based on samples that met the criteria of a viral load ≥ 1000 copies per ml.
2.3. Viral Load and Viral RNA Isolation
Viral load (plasma HIV-1 RNA) was quantified using b-DNA method, Versant HIV I-RNA 3.0 (Bayer Diagnostics, Emeryville, CA, USA). Viral RNA particles was extracted from 48 viral load plasma samples that met the criteria of detectable viral load copies greater or equal to 1000 copies per ml using QIA amp viral RNA mini kit, following the manufacturer’s protocol (Qaigen Inc. Chatworth, CA).
2.4. HIV-1 Genotyping
HIV-1 genotyping for the pol-gene region was done using the CDC in-house genotyping system to generate 1050 base pair fragment consisting of 1 - 99 amino acids of the protease region and codons 1 - 251 of the RT region. Detailed procedure for this assay has previously been described [18] [19] . Briefly, 15 µl of nucleic acid extracts from plasma were subjected to a one-step RT-PCR using RT forward and reverse primers and SuperScriptTM III one step RT-PCR systems with Platinum® Taq high fidelity following the manufacturer’s protocol (Invitrogen, Carlsbad, CA). For nested PCR, 4 µl of the RT-PCR product was used in the 2nd round nested PCR using inner primers to yield ~1062 base pair amplicon. Sequencing was then done using six overlapping primers and analyzed using the ABI 3100 genetic analyzer. Confirmation of base-calling and sequence editing was done using the Sequencher V 4.6 (Gene codes) software while drug resistance interpretation was performed using the Stanford genotyping drug resistance interpretation algorithm (V 4.2.6) (http://hivdb.stanford.edu/) whilst using the International AIDS Society 2009 mutation list for confirmation [20] .
2.5. Phylogenetic Analyses
Sequence alignment was done using sequence alignment editor (Se-Al V 2.0) to generate base pairs sequences devoid of gaps [21] . Phylogenetic analysis was then performed with Phylogenetic Analysis Using Parsimony (PAUP) Version 4.0 to generate neighbor-joining phylogenetic trees [22] . National Center for Biotechnology Information (NCBI) and Rega HIV subtyping tool version 2.0 was used to further evaluate accurately the classification of recombinant and non-recombinant genotypes using boot scanning methods. The reliability of the tree was tested with 1000 bootstrap analyses.
2.6. Nucleotide Sequence Accession Numbers
Sequences generated from this study have been deposited in the GenBank database with accession numbers HQ 845921 to HQ 845960. The GenBank numbers are used in this report.
2.7. Statistical Analysis
Statistical analyses were performed using STATA version 10.0 (Stata Corporation, College Station, TX). Differences between groups were compared using Pearson’s Chi-squared test. All statistical tests were two-sided and P value < 0.05 was considered statistically significant.
3. Results
3.1. Patient Characteristics
Of the 202 HIV infected patients, female patients were 130 (64%). The median age of the patients was 33 yrs old (interquartile range (IQR): 28 - 40). The median (IQR) CD4 cell count and plasma viral load level were 348 (217 - 535) cells/µl and 66 (50 - 766) copies/ml respectively (Table 1). ART patients sample analyzed were 160 (79.2%) were while 42 (20.8%) were drug-naïve. Among the ART patients, 153 (94.6%) were on first line regimen of highly active antiretroviral therapy (HAART) based on the drugs available in the National guideline (combination of 3 ARVs of two nucleoside/nucleotide reverse transcriptase inhibitors and one non-nucleoside reverse transcriptase inhibitor). Those on second line were 7 (5.4 %). Median month on treatment was 14 months (IQR: 8 - 25) (Table 1).
3.2. Evaluation of Drug Resistance Mutations
Of the 48 samples that met genotyping criteria, 40 were successfully sequenced (15 ART and 25 drug-naïve). Drug resistance mutation was identified in 5 ART patients (5/15, 33%) and 2 drug-naïve (2/25, 8%) at the RT region of the pol gene (Table 2). No mutation conferring resistance to Protease inhibitors (PIs) were seen in this population. Of the resistance mutations observed 6 (6/7) confers resistance to nucleotide reverse transcriptase inhibitors (NRTIs) and 6 (6/7) non-nucleotide reverse transcriptase inhibitors NNRTIs respectively of the RT region. Most commonly occurring mutation was M184V (5/40, 13%) and K103N (5/40, 13%) which confers resistance to NRTIs and NNRTIs respectively. Complex multi-drug resistance mutations Q151M occurred with other secondary mutations V75M and F77L in one of the ART patient; while thymidine analogue mutations (TAMs): M41L, L210W and T215Y occurred in two of the ART patient. Other NNRTI mutations observed were Y181C (2), G190G (1) (Table 2 & Table 3).
Among the ARV naive patients, singleton M41L and K103N resistant mutations conferring resistance to
Table 2 . Summary of resistance profile in ART experienced and naïve patients.
AZT; Zidovudine, 3TC; Lamivudine, NVP; Nevirapine, TDF; Tenofovir; FCT; Emtricitabine, LPV/TRV; Lopinavir-Ritonavir, d4T; Stavudine, ABC; Abacavir. Italic indicates key mutations.
Table 3. Resistance mutations selected by drug.
RT; Reverse Transcriptase, NRTI; Nucleotide Reverse Transcriptase Inhibitors, Non-Nucleotide Reverse Transcriptase Inhibitors.
NRTI and NNRTI respectively were identified.
Polymorphic substitutions were observed in both the PR: I13V, K20I, M36I, R41K, H69K, V82I and L89M and the RT: R211K, A98G/S/AS regions of the pol gene. Subtype specific polymorphisms A98AS/S/G (4/20; 20%;) and V82I (18/20; 90%) were observed among subtype G strains.
3.3. Genetic Variability
Phylogenetic analysis revealed Group M isolates by different clusters (Figure 1). The predominant subtype was subtype G (50%, 20/40), followed by circulating recombinant forms CRF02_AG (35%, 14/40,) CRF06_cpx (7.5%, 3/40), CRF-18-cpx (2.5%, 1/40), subtype J (2.5%, 1/40), and a unique recombinant AD (2.5%, 1/40). Subtype J and CRF06_cpx were the two strains found in the drug-naïve individuals.
4. Discussion
This study revealed key resistance mutations in HIV infected patients and also established that diverse recombinant forms of HIV are on the increase in Nigeria.
Analysis of patients ARVs regimen status revealed that majority of the patients (94.6%) were on first line regimen of highly active antiretroviral therapy (HAART) while the rest were on second line treatment based on the drug available as in the National guideline. The nucleoside backbone used for treatment in this setting are, AZT/ 3TC, TDF/FTC, d4T/3TC and ABC/3TC. Lopinavir/ritonavir serves as the protease inhibitors for the
Neighbor Joining tree of 40 Nigerian samples sequenced are reported and preceded by the prefix “NG,” underline, indicating the country (NG—Nigeria) of specimen collection and other code that follows represented the code for site where samples were collected and the hospital assigned number (e.g., NG.DOC.13741). The scale bar indicates an evolutionary distance of 10 nucleotides per position in the sequence.
Figure 1. Phylogenetic analysis of Nigerian HIV-1 sequences of the pol region.
nucleoside backbone combination of the second line. Virtually all the regimen except Tenofovir based has Lamivudine (3TC) in their combinations.
Analysis of the sequence revealed mutations conferring resistance at the RT region of the pol gene. This agrees with extensive use of NRTI in the first line regimen for the treatment of infected individuals in Nigeria. Mutations conferring resistance observed were M41L, Q151M, M184V, L210W, T215Y, K103N, and Y181C. Absence of resistance mutation at the PR region may be due to limited use of PI in Nigeria and only found among the second line regimens. M184V resistance mutation was commonly observed in the RT region among the ART experienced patients. This mutation is selected by use of lamivudine (3TC), Emtricitabine (FTC), Abacavir (ABC) but confers high level resistance to lamivudine and is usually the first to develop in isolates from patients receiving incompletely suppressive lamivudine-containing regimens [23] . Thymidine analogue mutations (TAMs) M41L, L210W and T215Y were observed in two of the patients on first line regimen of CRF_02AG strains. Presence of TAMs calls for public health concern as it promotes resistance to almost all NRTI analogues and is most often selected by zidovudine and stavudine-containing regimens. The two patients NG.SCBH.8960 and NG.DOC.13816 already failing therapy as the current mutation is selected by zidovudine based regimen.
Q151M a complex resistance mutation was detected in a subtype G treated patient (NG.DOC.10208) on second line regimen and occurred together with secondary mutation at positions 75 and 77. Q151M confer high level resistance to most NRTIs except Lamivudine and Tenofovir [24] [25] thus its presence in the study is also a major public health concern. Presence of Q151M mutation has not been highlighted by past studies in Nigeria.
Mutation conferring resistance to NNRTIs were seen in 6 of the samples occurring as K103N (12.5%) and Y181C (5%) as well as G190AG (2.5%) and M230L (2.5%). K103N was the second common mutation after M184V and is known to be the predominant mutation in patients receiving NNRTIs [26] and may cause cross resistance to all other NNRTI drug [27] .
Y181C is known to cause resistance to NNRTIs mostly nevirapine and delavirdine [28] . G190AG found in sample NG.DOC.10208 is an uncommon mutation, but G190A which is a close mutation is known to have high level resistance to nevirapine and an intermediate levels of resistance to efavirenz [28] [29] .
Drug resistance mutations were detected in the RT gene of 2 (2/25; 8%) drug-naïve persons. Circulation of HIV isolates containing resistance mutations in drug-naïve individuals was previously described in Nigeria [30] [31] and is quite worrisome as transmitted drug resistance will affect the choice of the most effective treatment regime. Early data suggested a 10% to 17% baseline resistance to major classes of ARVs in viruses from Nigerians not undergoing ART [16] [17] . Other studies have observed that about 10% of newly infected patients’ harbors virus that is at least partially resistant to one or more antiretroviral drugs [32] . One of the drug-naïve patient (NG.DOC.13742) having CRF-06-cpx displayed a singleton M41L mutation at the NRTI and the second one (NG.DOC.13805) having subtype J revealed K103N mutation at NNRTI region. M41L is known to be associated with development of resistance to AZT [33] , while V60I, an unusual mutation was found at the RT region. Some studies have mentioned that viruses with a single resistance associated amino acid substitution (singleton mutation) may constitute a particular interpretation problem as some of them may appear both as natural sequence polymorphism and as drug resistance mutation, and that M41L may represent a rare “natural” sequence variant when it occurs on its own [34] . It is estimated that every possible single point mutation occurs between 104 and 105 times per day in an untreated HIV-1-infected individual and that double mutants also occur commonly [35] . Therefore, most drug-resistance mutations are probably present prior to the start of therapy. The mutational patterns noticed in the drug-naïve individuals were not complex. This was also reported similarly in other studies both from Nigeria [30] and other African countries [36] [37] .Transmission of drug resistance from treatment-experienced patients to newly infected persons has been observed repeatedly in countries with long term use of antiretroviral drug [38] thus, raising major public health concerns [39] . Therefore clinical management of patients who have never been exposed to therapy may be unfavourably affected [40] . The occurrence of this intrinsic drug-resistant virus in HIV drug-naïve patients could imply that many drug experienced patients continue to engage in risk related behaviour. The public health implications of this include the need to intensify programs directed at such patients to address prevention of transmission.
Other mutations/polymorphic substitutions observed in both the PR (I13V, K20I, M36I, R41K, H69K, V82I and L89M) and the RT (R211K) regions of the pol gene. Mutations at these positions do not cause high level drug resistance, but they contribute to drug resistance when present together with certain primary protease mutations [41] . The higher frequency of protease polymorphism in non-B isolates, including protease positions 20, 36, 63, 82 and 93, has raised concern that PI treatment of non-subtype B infected persons could be less effective than that of subtype B [7] . A similar mutations/polymorphic substitutions in the protease region of some Nigerian isolates had been reported [30] .
The Phylogenetic analysis agrees with results of previous studies in Nigeria that the predominant variants of HIV-1 are subtype G and the circulating recombinant form CRF02_AG [17] [42] . However subtype A which had been reported as the second most prevalent in the country was not found in the study, this may be due to the limited sample size used for the study. Subtype J, a known rare subtype was observed in one of the drug-naïve. Meanwhile, the continuous presence of CRF02_AG as compared to previous study is suggesting the dominance of these recombinants in the country and in West Africa [43] and the epidemiologic importance in future research. Increasing number of recombinants was noticed in the study, the CRF06_cpx and a new form CRF18_ cpx. CRF06_cpx were found in two drug-naïve individuals. This is a complex mosaic virus that had been found in continents like Europe and Australia and likewise in some West Africa Countries such as Ivory Coast, Niger, Mali and Senegal [42] [44] .
CRF18_cpx, a novel HIV-1 circulating recombinant form was found in the North Eastern region of Nigeria, basically on the border boundary to Cameroun. This virus had been known to originate from Central Africa and had been found in Cameroun [45] . Cameroun share boundary with Nigeria with high diversity of circulating HIV-1 genetic forms and this closeness may likely influence the growing increase in HIV diversity in Nigeria. The identification of novel CRF may have practical implications in the management of patients and also in development of vaccine.
A unique recombinant form (URFs) of subtype A and D was observed in one of the sample; URF of subtype AD is known to circulate in East Africa [46] and Ghana [47] . Its presence may reflect an incidence of coinfections or superinfections with diverse viral strains thus increasing the genetic complexity of the HIV-1 epidemic in West Africa.
A major limitations of this study is that it was cross sectional, thus patients were not followed up to assess the outcome of their response to treatment and if the resistance mutations seen in the drug-naïve will be observed later in the course of treatment. The criteria for genotyping selection of a viral load ≥ 1000 copies per ml may have missed some mutation and variant in circulation.
The findings provide current information with regard to drug resistance profile among treatment experienced and naive population as well as the genetic diversity which is helpful in guiding the use of the different ARV drugs in use in Nigeria.
5. Conclusion
In summary, our results indicate that HIV-1 genetic diversity in Nigeria is increasing and presence of drug resistance mutations agrees with components of the first-line drug taken by majority. Thus, continuous assessment is needed to monitor the evolving diverse strains and potential impact of drug resistance on ART programs in Nigeria in order to ensure continuance efficacy of the limited available ART.
Acknowledgements
The authors would like to acknowledge the technical support provided by Dr. Clement Zeh and his HIV Research laboratory team at Centre for Disease Control and Prevention, HIV Research Laboratory, Kisumu, Kenya for the considerable help and support in this work. I also wish to thank the laboratory staff of the following facilities who assisted in the collection and initial processing of the samples: St. Charles Borromeo, Onitsha, Nigeria; Muritala Mohammed Specialist Hospital, Kano; Mambilla Baptist Hospital, Gembu, Taraba; Daughters of Charity, Kubwa, Abuja (DOC); Federal Medical Center, Yola; Holy Family Catholic Hospital, Ikom; Mainland Hospital, Yaba and Maitama District Hospital, Abuja. The views expressed herein is that of the authors and do not necessarily reflect those of any institution.
Authors’ Contributions
Conceived and designed of experiment: ORN-M, SOO, LEO, HII, SMA; Performed the experiment: ORN-M, ANK, SCI, analyzed the data: ORN-M, ANK, SCI, TB, HAM, SMA; wrote the paper: ORN-M; reviewed the manuscript: ORNM, S.C. Inzaule, EAO, HAM, SMA.
Ethical Approval
Ethical approval for the main study was obtained from ethical review committee for Ahmadu Bello University. Written informed consent was obtained from each participant prior to the study initiation.
Competing Interest
The authors declare that there is no financial or personal relationship(s) that may have inappropriately influenced the writing of this article.
NOTES
*Corresponding author.