Genetic diversity and drug resistance of HIV-1 among infected pregnant women newly diagnosed in Luanda, Angola

Abstract

Monitoring genetic diversity and drug resistance mutations (DRMs) is critical for understanding HIV epidemiology. Here, we report HIV-1 genetic diversity and DRMs in blood samples from 42 HIV-positive pregnant women naive to antiretroviral therapy (ART), in Luanda. The samples were subjected to nested-PCR, followed by sequencing of HIV-1 pol gene, targeting the protease and reverse transcriptase fragments. HIV-1 diversity was analyzed using the REGA HIV-1 subtyping tool and DRMs were identified using the Calibrated Population Resistance tool. A total of 34 sequences were obtained. The data revealed wide HIV-1 subtypes heterogeneity, with subtype C (38%, 13/34) the most frequent, followed by the subtypes F1 (18%, 6/34), A1 (9%, 3/34), G (9%, 3/34), D (6%, 2/34) and H (3%, 1/34). In addition, recombinants strains were detected, with CRF02_AG (6%, 2/34) the most frequent, followed by CRF37_cpx, F1/C, A1/G and H/G, all with 3% (1/34). A total of 6/34 (18%) of the sequences presented DRMs. The non-nucleoside reverse transcriptase inhibitors presented 15% (5/34) of resistance. Moreover, 1/34 (3%) sequence presented resistance against both non-nucleoside reverse transcriptase inhibitors and nucleoside reverse transcriptase inhibitors, simultaneously. Despite the small sample size, our results suggest the need to update currently used ART regimens. Surveillance of HIV-1 subtypes and DRMs are necessary to understand HIV epidemiology and to guide modification of ART guidelines in Angola.

About the Authors:

Cruz S. Sebastião

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Writing – original draft, Writing – review & editing

Affiliations NOVA Medical School, Faculdade de Ciências Médicas, Universidade NOVA de Lisboa, Lisboa, Portugal, Laboratório de Biologia Molecular, Instituto Nacional de Investigação em Saúde, Luanda, Angola, Instituto Superior de Ciências da Saúde, Universidade Agostinho Neto, Luanda, Angola, Centro de Investigação em Saúde de Angola, Luanda, Angola

Zoraima Neto

Roles Investigation, Methodology, Writing – review & editing

Affiliation: Laboratório de Biologia Molecular, Instituto Nacional de Investigação em Saúde, Luanda, Angola

Carlos S. de Jesus

Roles Methodology, Supervision, Writing – review & editing

Affiliation: Laboratorio de AIDS & Imunologia Molecular, Instituto Oswaldo Cruz, Rio de Janeiro, Brazil

Marinela Mirandela

Roles Data curation, Investigation, Methodology, Writing – review & editing

Affiliation: Laboratório de Biologia Molecular, Instituto Nacional de Investigação em Saúde, Luanda, Angola

Domingos Jandondo

Roles Investigation, Methodology, Writing – review & editing

Affiliation: Laboratório de Biologia Molecular, Instituto Nacional de Investigação em Saúde, Luanda, Angola

José C. Couto-Fernandez

Roles Data curation, Methodology, Writing – review & editing

Affiliation: Laboratorio de AIDS & Imunologia Molecular, Instituto Oswaldo Cruz, Rio de Janeiro, Brazil

Amilcar Tanuri

Roles Formal analysis, Methodology, Writing – review & editing

Affiliation: Laboratório de Virologia Molecular, Departamento de Genética, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brasil

Joana Morais

Roles Conceptualization, Data curation, Formal analysis, Investigation, Writing – review & editing

Affiliations Laboratório de Biologia Molecular, Instituto Nacional de Investigação em Saúde, Luanda, Angola, Centro de Investigação em Saúde de Angola, Luanda, Angola, Faculdade de Medicina, Universidade Agostinho Neto, Luanda, Angola

Miguel Brito

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Supervision, Validation, Writing – review & editing

* E-mail: [email protected]

Affiliations Centro de Investigação em Saúde de Angola, Luanda, Angola, Health and Technology Research Center, Escola Superior de Tecnologia da Saúde de Lisboa, Instituto Politécnico de Lisboa, Lisboa, Portugal

ORCID logo http://orcid.org/0000-0001-6394-658X

Introduction

The human immunodeficiency virus (HIV) has become a major global public health problem [1], affecting about 36.9 million people in the world [2]. In Angola, a total of 310,000 cases were reported in 2018 [2]. HIV is classified into types (HIV-1 and HIV-2), groups (M, N, O and P), subtypes (A-D, F-H, J and K), sub-subtypes (A1, A2, F1 and F2), circulating recombinant forms (CRFs) and unique recombinant forms (URFs) [3]. HIV-1 is responsible for the vast majority of HIV infections [4]. All subtypes of HIV-1 group M (except B), several CRFs and URFs have been described in Angola [5–11].

Universal access to antiretroviral therapy (ART) has successfully decreased mortality and morbidity associated with HIV [2,12]. The first-line of the ART drugs used in Angola includes the nucleoside reverse transcriptase inhibitors (NRTIs), tenofovir (TDF) and lamivudine (3TC), and a non-nucleoside reverse transcriptase inhibitor (NNRTI), either efavirenz (EFV) or nevirapine (NVP) [13,14]. In addition, zidovudine (AZT) has been used to prevent vertical transmission [13,14].

The emergence of HIV-1 subtypes with drug resistance mutations (DRMs) during pregnancy represents a challenge for the efficacy of ART, especially in low- and middle-income countries [15]. There is a lack of recent data on HIV-1 genetic diversity and prevalence of DRMs in Angola [15,16]. In this study, we investigated the genetic diversity and DRM prevalence in blood samples from HIV-positive pregnant women naive to ART in Luanda, to better understand HIV epidemiology and to allow a timely modification of ART guidelines in Angola.

Materials and methods

Study design and sample collection

A cross-sectional study was carried out at the Lucrecia Paim Maternity clinic, located in Luanda, capital city of Angola, during the months of April to June of 2018. The study involved 1612 pregnant women who were screened for HIV infection using the rapid antibody detection test Determine HIV1/2™ (Alere, Japan) and the Unigold™ HIV (Trinity Biotech, Ireland) during prenatal care. Sociodemographic characteristics and blood samples were collected from HIV-positive pregnant women. The main criterion for inclusion of HIV-positive pregnant women was that they had not been previously exposed to any ART. The blood samples were collected in a tube with EDTA, centrifuged and the plasma was aliquoted and stored at -80°C. The blood samples preparation was performed at the Molecular Biology Laboratory, of the National Institute for Health Research of Angola (INIS). Following the recommendations of the National Institute of Fighting against AIDS (INLS), the HIV-positive women, were prescribed ART with TDF, 3TC and EFV, and were medicated with AZT until child birth [13,14].

RNA extraction, cDNA synthesis, PCR and sequencing

Total viral RNA was extracted from 140μL of plasma using QIAamp Viral RNA kit (QIAGEN, Germany) following the manufacturer instructions. The cDNA synthesis was carried out using 10μL of the RNA in a final reaction volume of 20μL. The mix contained 25mM DNTP mix, 5X M-MLV buffer, 10mM of dithiothreitol (DTT), 40U of RNase OUT™ (Life Technologies, USA), 0.1mM of MMRTR6 primer (5’-TTTTACATCATTAGTGTGGG-3’), and 200U of M-MLV enzyme (Life Technologies, USA) [17].

The obtained cDNA was subjected to a nested-PCR, targeting the protease (PR) and reverse transcriptase (RT) fragments of the HIV-1 pol gene, with an expected size of 1302 bp, using the protocol previously described [17]. Successful amplification was checked using a 1% agarose gel. The amplicons were purified using the NZYGelpure Kit (Nzytech, Portugal), and sequenced using the ABI BigDye Terminator v3.1 reaction kit (Applied Biosystems, USA). For each sample, eight primers were used for the complete sequencing of the PR (nucleotide range: 2253–2549) and the first 335 codons of RT (nucleotide range: 2550–3554), considering the genome of the HXB2 strain (nucleotide range: 2252–3554) [17]. Sequencing was performed on an ABI 3500 sequencer (Applied Biosystems, USA) at the Molecular Biology Laboratory of the INIS, in Luanda.

HIV-1 subtyping, phylogenetic and resistance mutation analysis

The electropherograms were analyzed using the software RECALL v2.25 [18]. Classification of HIV subtypes was conducted using the REGA HIV-1 subtyping tool v3.0 (http://dbpartners.stanford.edu:8080/RegaSubtyping/stanford-hiv/typingtool/) [19]. The nucleotide sequences obtained were aligned with HIV-1 M-group nucleotide sequences downloaded from NCBI (https://www.ncbi.nlm.nih.gov/nuccore/) and Los Alamos (https://www.hiv.lanl.gov/content/index) databases. The sequences obtained in the study were deposited in GenBank (NCBI) and were assigned the accession numbers MK543512 to MK543545. A phylogenetic tree was inferred using the Neighbor-Joining (NJ) method [20], with Tamura-Nei genetic distances [21]. Clade support was assessed using 1000 bootstrapped replicates [22]. The phylogeny was estimated using MEGA software v7.0 [23], and recombinant viruses were characterized by boot scanning using SimPlot [24].

Drug resistance mutations to HIV-1 were identified using the Calibrated Population Resistance tool (CPR) v8.0 (https://hivdb.stanford.edu/cpr/form/PRRT/) [25]. Additionally, analyses of drug resistance profile were performed using the Stanford genotypic resistance interpretation algorithm (https://hivdb.stanford.edu/hivseq/by-sequences/) [26].

Statistical analysis

Chi-square (X2) tests were performed to evaluate the association between sociodemographic characteristics and HIV prevalence at 5% statistical significance in the SPSS v25 statistical program (IBM SPSS Statistics, USA).

Ethical considerations

The pregnant women were informed of the study and consented (oral and written) to participation and follow-up until delivery. Consent from parents or guardians of the minors (under 15 years) was also obtained. All pregnant women underwent pre- and post-test counseling individually. The HIV study results were provided to the clinical staff to ensure appropriate patient clinical management. The study protocol was reviewed and approved by the National Ethics Committee of Angola (nr.13/2018), general directorate from Lucrecia Paim Maternity clinic (nr.083/GDG/MLP/2018) and the Ethics Research Committee of the NOVA Medical School/Faculdade de Ciências Médicas—Universidade NOVA de Lisboa (nr.51/2019/CEFCM).

Results

Sociodemographic characteristics

From the 1612 pregnant women tested for HIV, 42 (2.6%) tested positive and not been exposed to any ART previously. All pregnant women were from Luanda province. A total of 25/42 HIV-positive pregnant women were in the age group of 25–34 years, 18/42 had basic and secondary education, respectively, and 19/42 were unemployed. Moreover, a total of 12/42 pregnant women were in the second trimester of gestation, followed by 11/42 in the third trimester and 10/42 pregnant women diagnosed and the sample was obtained just before labor (Table 1).

[Figure omitted. See PDF.]

Table 1. Sociodemographic characteristics of pregnant women tested for HIV in Luanda, Angola, 2018.

https://doi.org/10.1371/journal.pone.0225251.t001

Genetic diversity analysis

From the 42 plasma samples subjected to nested-PCR, a total of 34 amplicons and respective sequences were obtained. It was not possible to obtain amplicons from the remaining 8 samples, even after repeated PCR attempts using different primers. The genotyping analysis revealed that 28/34 sequences were HIV-1 pure subtypes and 6/34 recombinants strains. From the detected pure subtypes, subtype C (38%, 13/34) was the most frequent HIV-1 subtype, followed by the subtypes F1 (18%, 6/34), A1 (9%, 3/34), G (9%, 3/34), D (6%, 2/34) and H (3%, 1/34). From the recombinants strains detected, CRF02_AG (6%, 2/34), was marginally more frequent, followed by CRF37_cpx, F1/C, A1/G and H/G, all with 3% (1/34) (Fig 1).

[Figure omitted. See PDF.]

Fig 1. HIV-1 diversity among infected pregnant women in Luanda, Angola, 2018.

The analysis involved 34 nucleotide sequences. (A) Frequency of HIV-1 subtypes and recombinants. (B) Recombinants distribution. The subtypes were identified by the REGA v3.0, and recombinants strains by SimPlot.

https://doi.org/10.1371/journal.pone.0225251.g001

Genetic distance analysis showed genetic similarity (more than 70%) of HIV-1 subtype C isolates with isolated from Botswana, Mozambique, Tanzania, South Africa and India. The subtype F1 was more similar to isolates from Brazil and Spain. The subtype A1 was most similar to isolates from Cameroon, South Africa, and Pakistan. The subtype G was most similar to isolates from Portugal, and subtype D, subtype H, and recombinant strain H/G to isolates from Democratic Republic of Congo. The recombinants strains A1/G, CRF02_AG, and CRF37_cpx were similar to isolates from Cameroon and Ivory Coast (Fig 2).

[Figure omitted. See PDF.]

Fig 2. Evolutionary relationships of HIV-1 subtypes among infected pregnant women in Luanda, Angola, 2018.

The analysis involved 56 nucleotide sequences. All HIV-1 strains identified in this study were from M-group. The phylogenetic tree was inferred using the NJ method, and Tamura-Nei genetic distances with 1000 bootstrap replicates. Analyses were conducted in MEGA v7.0. Bootstrap values more than 70% and recombinants strains are indicated.

https://doi.org/10.1371/journal.pone.0225251.g002

Antiretroviral resistance mutations and resistance profile

The sequence analysis showed that 6/34 (18%) of the sequences presented DRMs. Of these, 5/34 (15%) sequences presented resistance to specific NNRTIs. Moreover, 1/34 (3%) sequence presented resistance against both NNRTIs and NRTIs, simultaneously. Mutations at positions K103N (2/34), G190A (2/34), Y181I (1/34), and P225H (1/34) were observed against the class of NNRTIs, whereas mutations at positions M41L, D67N, T69D and T215S were observed against NRTIs. A non-polymorphic protease inhibitor (PI) selected mutation was observed at position I85V. High levels of drug resistance were observed for NVP (5/34) and EFV (2/34), intermediate resistance for EFV (3/34) and AZT (1/34) and low resistance for Abacavir (ABC) (1/34) and TDF (1/34). In addition, the susceptibility profile associated to the second generation of NNRTIs was also identified and 1/34 and 3/34 sequence presented high- and low-resistance to Etravirine (ETR) and Rilpivirine (RPV), respectively (Table 2).

[Figure omitted. See PDF.]

Table 2. Drug resistance mutations to NRTI, NNRTI and PI according to the HIV-1 subtypes and resistance profile.

https://doi.org/10.1371/journal.pone.0225251.t002

Discussion

This study presents an important update on molecular epidemiology of circulating HIV-1 strains in Luanda. The HIV prevalence was 2.6% (42/1612). A significant difference in HIV positivity was observed between pregnant women of different age groups and pregnancy stage (P<0.05). On the other hand, residence, level of education and occupation did not show significant differences (P>0.05) (Table 1).

Genetic analysis of 34 isolates revealed a wide diversity of HIV-1 strains (Fig 1), similar to that observed in previous studies performed in Angola [5–11]. Though it is hard to prove statistical significance with our small sample size, our study indicates an increase in HIV-1 subtype C, but a slight decrease in subtype F1 in Luanda [5,8]. The genetic similarity of HIV-1 subtype C with isolates from Botswana, Mozambique, Tanzania, South Africa, India and the subtype F1 with isolates from Brazil and Spain (Fig 2), may be attributed to high mobility between countries or to the fact that, after colonial war, and the end of civil war in 2002, thousands of refugees returned to Angola [27].

The RT inhibitors are important components of ART regimen [12–14]. The identification of pre-treatment drug resistance in pregnant women naïve to ART (Table 2), may threaten ART based strategy to HIV control in Luanda [12–14]. The K103N and G190A mutations are associated with EFV and NVP resistance [28]. The Y181I and P225H mutations are often associated with second generation RPV and ETR resistance [29]. The thymidine analogue associated mutations (TAMs) at positions M41L and D67N have the greatest impact on susceptibility of AZT and Stavudine (d4T) [28]. The T69D mutation when present with T215S mutation, is associated with broad resistance to NRTIs [28]. The non-polymorphic PI-selected mutation at position I85V has minimal effects on PIs susceptibility [30,31].

The identification of K103N and G190A mutations, may be attributed to the long use of NNRTIs as part of the first-line ART regimens in Angola [13,14]. Displacement of people to countries where ART is available for a longer period also may help to explain the origin of the HIV-1 subtypes with DRMs in Luanda [27,32].

The reasons for PCR failure of 8/42 samples were not identified. It may be that the high HIV-1 genetic diversity compromised the binding of primers, even though they were targeted at highly conserved regions of the HIV-1. Other studies with more representative sampling monitoring of the HIV-1 subtypes and DRMs are necessary to guide a timely modification of ART guidelines in Angola. Despite the small sample size, our findings suggest that the Angolan Ministry of Health should prompt consideration of moving to generic integrase strand transfer inhibitors (INSTIs) in the ART regimen in Angola [33].

Conclusions

Our results show a wide HIV-1 subtypes heterogeneity, with subtype C the most frequent. A total of 6/34 (18%) of the sequences presented DRMs. Of these, 15% (5/34) were associated with resistance against NNRTIs. Moreover, 1/34 (3%) sequence presented resistance against both NNRTIs and NRTIs, simultaneously. Our findings suggest the need to update currently used ART regimens. Better understanding is needed of emergence of HIV-1 subtypes and DRMs, to allow a timely modification of ART guidelines in Angola.

Acknowledgments

Thanks to pregnant women who participated in this study. To thank the Postgraduate Science for Development Program (PGCD) and the Foundation for Science and Technology (FCT), by the PhD scholarship awarded to Cruz Sebastião (SFRH/BD/135296/2017), and the project Pró-África (CNPq: 440145/2015-5) for supporting some analysis. Moreover, thanks to INIS, CISA, INLS, FIOCRUZ, NOVA Medical School/Faculdade de Ciências Médicas, ISCISA/Agostinho Neto University, and Lucrecia Paim Maternity clinic, for institutional support. Thanks to Valter Nuaila, Marcio Sitoe, Fredilson Melo and Sarah Hill, for scientific support.

Citation: Sebastião CS, Neto Z, de Jesus CS, Mirandela M, Jandondo D, Couto-Fernandez JC, et al. (2019) Genetic diversity and drug resistance of HIV-1 among infected pregnant women newly diagnosed in Luanda, Angola. PLoS ONE 14(11): e0225251. https://doi.org/10.1371/journal.pone.0225251

References

1. United Nations. Political Declaration on HIV and AIDS: On the Fast Track to Accelerating the Fight against HIV and to Ending the AIDS Epidemic by 2030. 2016;17020: 1–26.

2. United Nations Joint Programme on HIV/AIDS (UNAIDS). Unaids Data 2018. 2018.

3. Leitner T, Hahn B, Mullins J, Wolinsky S, Foley B, Apetrei C, et al. HIV Sequence Compendium 2015 Editors. Theor Biol Biophys Los Alamos Natl Lab. 2015.

4. Sharp PM, Hahn BH. Origins of HIV and the AIDS pandemic. Cold Spring Harb Perspect Med. 2011;1: a006841. pmid:22229120

5. Afonso JM, Bello G, Monick GL, Sojka M, Morgado MG. HIV-1 Genetic Diversity and Transmitted Drug Resistance Mutations among Patients from the North, Central and South Regions of Angola. PLoS One. 2012;7: 11. pmid:22952625

6. Bártolo I, Rocha C, Bartolomeu J, Gama A, Marcelino R, Fonseca M, et al. Highly divergent subtypes and new recombinant forms prevail in the HIV/AIDS epidemic in Angola: New insights into the origins of the AIDS pandemic. Infect Genet Evol. 2009;9: 672–682. pmid:18562253

7. Bártolo I, Zakovic S, Martin F, Palladino C, Carvalho P, Camacho R, et al. HIV-1 diversity, transmission dynamics and primary drug resistance in Angola. PLoS One. 2014;9: 1–17. pmid:25479241

8. Castelbranco EPAF, da Silva Souza E, Cavalcanti AMS, Martins AN, de Alencar LCA, Tanuri A. Frequency of Primary Resistance to Antiretroviral Drugs and Genetic Variability of HIV-1 Among Infected Pregnant Women Recently Diagnosed in Luanda-Angola. AIDS Res Hum Retroviruses. 2010;26: 1313–1316. pmid:20929349

9. Bártolo I, Epalanga M, Bartolomeu J, Fonseca M, Mendes A, Gama A, et al. High Genetic Diversity of Human Immunodeficiency Virus Type 1 in Angola. AIDS Res Hum Retroviruses. 2005;21: 306–310. pmid:15943573

10. Abecasis A, Paraskevis D, Epalanga M, Fonseca M, Burity F. HIV-1 genetic variants circulation in the North of Angola. 2005;5: 231–237. pmid:15737914

11. Garrido C, Zahonero N, Fernándes D, Serrano D, Silva AR, Ferraria N, et al. Subtype variability, virological response and drug resistance assessed on dried blood spots collected from HIV patients on antiretroviral therapy in Angola. J Antimicrob Chemother. 2008;61: 694–698. pmid:18218644

12. World Health Organization. The use of antiretroviral drugs for treating and preventing hiv infection. World Heal Organ. 2016.

13. National Institute of Fighting against AIDS. Plano Estratégico Nacional Para o Controlo das Infecções de Transmissão Sexual, VIH e SIDA Instituto Nacional de Luta Contra a Sida. 2006.

14. National Institute of Fighting against AIDS. Normas De Tratamento Antirretroviral. 2015; 159.

15. Ssemwanga D, Lihana RW, Ugoji C, Abimiku A, Nkengasong J, Dakum P, et al. Update on HIV-1 acquired and transmitted drug resistance in Africa. AIDS Rev. 2015;17: 3–20. pmid:25427100

16. Bbosa N, Kaleebu P, Ssemwanga D. HIV subtype diversity worldwide. Curr Opin HIV AIDS. 2019;14: 153–160. pmid:30882484

17. Delatorre E, Silva-de-Jesus C, Couto-Fernandez JC, Pilotto JH, Morgado MG. High HIV-1 Diversity and Prevalence of Transmitted Drug Resistance Among Antiretroviral-Naive HIV-Infected Pregnant Women from Rio de Janeiro, Brazil. AIDS Res Hum Retroviruses. 2017;33: 68–73. pmid:27392995

18. Woods CK, Brumme CJ, Liu TF, Chui CKS, Chu AL, Wynhoven B, et al. Automating HIV drug resistance genotyping with RECall, a freely accessible sequence analysis tool. J Clin Microbiol. 2012;50: 1936–1942. pmid:22403431

19. Pineda-peña A, Rodrigues N, Imbrechts S, Libin P, Barroso A, Deforche K, et al. Infection, Genetics and Evolution Automated subtyping of HIV-1 genetic sequences for clinical and surveillance purposes : Performance evaluation of the new REGA version 3 and seven other tools. 2013.

20. Saitou N, Nei M. The Neighbor-joining Method: A New Method for Reconstructing Phylogenetic Trees. 1987;4: 406–425.

21. Tamura K, Nei M. Estimation of the Number of Nucleotide Substitutions in the Control Region of Mitochondrial DNA in Humans and. 1993;10.

22. Felsenstein J. Confidence Limits on Phylogenies: an Approach Using the Bootstrap. Int J Org Evol. 1985;39: 783–791. pmid:28561359

23. Kumar S, Stecher G, Tamura K, Medicine E. MEGA7 : Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. 2016; 1–11.

24. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, et al. Full-Length Human Immunodeficiency Virus Type 1 Genomes from Subtype C-Infected Seroconverters in India, with Evidence of Intersubtype Recombination. Am Soc Microbiol. 1999;73: 152–160.

25. Bennett DE, Camacho RJ, Otelea D, Kuritzkes DR, Fleury H, Kiuchi M, et al. Drug resistance mutations for surveillance of transmitted HIV-1 drug-resistance: 2009 update. PLoS One. 2009;4. pmid:19266092

26. Liu TF, Shafer RW. Web Resources for HIV Type 1 Genotypic-Resistance Test Interpretation. Clin Infect Dis. 2006;42: 1608–1618. pmid:16652319

27. Perrin L, Kaiser L, Yerly S. Travel and the spread of HIV-1 genetic variants. Lancet Infect Dis. 2003;3: 22–27. pmid:12505029

28. Johnson VA, Calvez V, Günthard HF, Paredes R, Pillay D, Shafer R, et al. 2011 Update of the Drug Resistance Mutations in HIV-1. Top Antivir Med. 2011; 156–164. pmid:22156218

29. Alcaro S, Alteri C, Artese A, Ceccherini-silberstein F, Costa G, Ortuso F, et al. Docking Analysis and Resistance Evaluation of Clinically Relevant Mutations Associated with the HIV-1 Non- nucleoside Reverse Transcriptase Inhibitors Nevirapine, Efavirenz and Etravirine. 2011; 2203–2213. pmid:21953939

30. Shahriar R, Rhee S, Liu TF, Fessel WJ, Scarsella A, Towner W, et al. Nonpolymorphic Human Immunodeficiency Virus Type 1 Protease and Reverse Transcriptase Treatment-Selected Mutations. 2009;53: 4869–4878. pmid:19721070

31. Rhee S, Taylor J, Fessel WJ, Kaufman D, Towner W, Troia P, et al. HIV-1 Protease Mutations and Protease Inhibitor Cross-Resistance. 2010;54: 4253–4261. pmid:20660676

32. Van de Vijver D, Wensing AM, Boucher C, others. The epidemiology of transmission of drug resistant HIV-1. Reviews. 2006;2007: 17–36.

33. Inzaule SC, Hamers RL, Doherty M, Shafer RW, Bertagnolio S, Rinke de Wit TF. Curbing the rise of HIV drug resistance in low-income and middle-income countries: the role of dolutegravir-containing regimens. Lancet Infect Dis. 2019;19: e246–e252. pmid:30902440