Prevalence of HIV-1 pretreatment drug resistance among treatment naïve pregnant women in Bissau, Guinea Bissau

Abstract

Background

With increased access to antiretroviral treatment (ART) in sub-Saharan Africa emergence of HIV-1 pretreatment drug resistance constitutes a serious risk. This may lead to rapid virological failure in subjects initiating ART, and mother-to-child transmission despite prophylaxis.

Methods

Treatment-naïve pregnant women from four antenatal care clinics in Bissau, Guinea-Bissau, were enrolled from October 2016 to November 2017. Genotypic resistance testing and phylogenetic subtype analysis was performed on 48 specimens.

Results

Forty eight women met the survey inclusion criteria. All specimens were successfully amplified and genotyped. Specimens from five women were associated with HIV-1 drug resistance mutations. Four carried mutations exclusively linked to non-nucleoside reverse transcriptase inhibitors (NNRTIs) (K103N, K103N/S) and one carried mutations to both NNRTIs (G190S, K101E) and nucleoside reverse transcriptase inhibitors (NRTIs) (M184V). These results corresponded to 10.4% (95% CI: 4.5–22.2%), 2.1% (95% CI: 0.4–10.9%) and 0% (95% CI: 0.0–7.4%) drug resistance mutations to NNRTIs, NRTIs and protease inhibitors, respectively. HIV-1 circulating recombinant form 02AG was most commonly found, followed by HIV-1 sub-subtype A3. Subtype/CRF was not associated with drug resistance mutations.

Conclusion

Our study reports a 10.4% prevalence of pretreatment drug resistance to NNRTIs in HIV-1-infected pregnant women in the capital Bissau, Guinea Bissau. Since NNRTIs are part of first-line ART in the country, baseline resistance screenings or adjustment of national treatment guidelines should be considered as antiretroviral treatment programs are scaled up.

About the Authors:

Sten Wilhelmson

Roles Conceptualization, Data curation, Formal analysis, Methodology, Project administration, Writing – original draft, Writing – review & editing

Affiliation: The Department of Translational Medicine, Lund University, Malmö, Sweden

Fredrik Månsson

Roles Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing

Affiliation: The Department of Translational Medicine, Lund University, Malmö, Sweden

ORCID logo http://orcid.org/0000-0002-2802-0420

Jacob Lopatko Lindman

Roles Data curation, Project administration, Supervision, Writing – review & editing

Affiliation: The Department of Clinical Sciences Lund, Lund University, Lund, Sweden

ORCID logo http://orcid.org/0000-0002-2034-464X

Ansu Biai

Roles Data curation, Methodology, Project administration, Writing – review & editing

Affiliation: The National Public Health Laboratory, Bissau, Guinea-Bissau

Joakim Esbjörnsson

Roles Conceptualization, Methodology, Writing – review & editing

Affiliation: The Department of Laboratory Medicine, Lund University, Lund, Sweden

Hans Norrgren

Roles Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Writing – review & editing

Affiliation: The Department of Clinical Sciences Lund, Lund University, Lund, Sweden

Marianne Jansson

Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing – review & editing

Affiliation: The Department of Laboratory Medicine, Lund University, Lund, Sweden

Patrik Medstrand

Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: The Department of Translational Medicine, Lund University, Malmö, Sweden

ORCID logo http://orcid.org/0000-0003-4259-8744

the SWEGUB CORE group

¶Group members are listed in Acknowledgments.

Introduction

Guinea-Bissau is a small multi-ethnic West African country with 1.8 million inhabitants, of whom 300,000 reside in the capital, Bissau [1]. The country ranks amongst the poorest in the world [1]. The human immunodeficiency virus type 1 (HIV-1) prevalence among adults, 15–49 years, is estimated to 3.1%, but is higher among women than men (3.8% vs. 2.4%) [2]. The national antiretroviral treatment (ART) program in Guinea-Bissau was initiated in 2005. As in most low and middle income countries (LMIC), ART follows the WHO guidelines with use of two nucleoside reverse transcriptase inhibitors (NRTIs) and one non-nucleotide reverse transcriptase inhibitor (NNRTI) [3]. About 35% of eligible HIV positive adults and 15% of the children received ART in 2016, while the proportion of pregnant women who received prevention from mother to child transmission (PMTCT) prophylaxis was 85% [2]. Furthermore, high numbers of loss to follow-up, mortality and development of HIV drug resistance (DR) during treatment (acquired DR) have been reported in the country [4].

The emergence of HIV-1 DR constitutes a threat to ART programs in resource-limited settings. Various factors drive the emergence of HIV-1 DR, including poor adherence and logistical barriers [5], high replication and mutation rates of HIV-1, and low genetic barrier for development of resistance to standardized first-line regimes [6]. Acquired DR as well as DR in treatment-naïve people (pretreatment DR) have increased in LMIC over the last decade [7, 8]. From a public health perspective it is important to monitor levels of acquired and pretreatment DR in order to forecast long term success of ART programs [9]. On an individual level pretreatment DR is of particular importance in LMIC, where virological monitoring and drug resistance screening is normally not part of patient management. Thus, individuals infected with drug resistant HIV-1 strains are at risk of not benefitting from the treatment, which may lead to disease progression and onward transmission.

Various HIV subtypes and circulating recombinant forms (CRF) cause infection in humans. Recombination between HIV subtypes and CRFs have been estimated to account for approximately 20% of the global HIV infections [10]. CRF02_AG is the most prevalent variant in several West African countries and accounts for 39–83% of the infections in this area [11–13]. Three major HIV subtypes/CRFs have been described in Guinea-Bissau: CRF02_AG, sub-subtype A3, and recombinants of CRF02_AG and A3 referred to as A3/02 [12, 13]. This may be important in the HIV epidemic since different viral variants have been linked to differences in viral load [14–16], disease progression rate [12], vertical transmission rate [17] and propensity to develop resistance to ART [18].

Documenting the profile of DR in HIV-1-infected pregnant women is crucial for improving the efficacy of maternal ART and prophylaxis in infants, and is also a preferred approach to estimate pretreatment DR (19). This can also help policy makers in the process of designing future national HIV treatment guidelines. Thus, in the context of increasing prevalence of acquired DR, and to gain an understanding of the effectiveness of contemporary ART in Guinea-Bissau, the aim of the current study was to estimate the level of pretreatment DR among pregnant women in the country. Moreover, since resistance data linked to information regarding HIV-1 subtypes and recombinants circulating among pregnant women have not been reported in Guinea-Bissau previously, this was also studied.

Methods

Study design and participants

Pregnant women who tested positive for HIV-1 in antenatal screening at four antenatal care clinics in the capital Bissau: Bairro Militar Health Centre, Antula Health Centre, Quelele Health Centre and Plack-II Health Centre, were asked for participation in the study. All participants finalized a questionnaire regarding previous HIV testing, antiretroviral treatment/mother-to-child prophylaxis as well as questions of a socio-economical nature. The survey aimed to follow the World Health Organization (WHO) recommended threshold survey methodology [19]. However, the WHO threshold survey methodology was under revision at the time of this study, and hence, we used the pre-revised guidelines. Original inclusion criteria were laboratory confirmation of HIV infection, age <25 years and no previous pregnancies. Due to frequent stock outs of HIV tests and a slower inclusion rate as a result, and in order not to prolong the study period, we omitted the age limit of 25 years and also included women with previous pregnancies. A total of 52 reportedly antiretroviral-naïve HIV-infected pregnant women were enrolled from October 2016 to November 2017. All participants that tested HIV positive were counselled and informed about antiretroviral treatment (ART), and were offered ART through the local health centre or at centralised services within Bissau City.

Sample management

Determine (Abbott Diagnostic Division, Hoofddorp, Holland) was used for pretreatment HIV diagnosis at the antenatal care clinics. Samples were transported to the laboratory for national health (LNSP) and analyzed for CD4 absolute count, CD4% and haemoglobin count using FACSPresto—Near Patient CD4 counter (Becton Dickinson, NYSE:BDX, USA). A confirmatory HIV-1 discriminatory test was performed using Geenius HIV 1/2 confirmatory assay (Bio-RAD). Plasma was separated from whole blood by centrifugation and stored at -20 C until transported on dry ice for further storing in -80 C and genotyping at the Clinical Virology section at Lund University, Sweden.

Drug resistance genotyping

RNA was extracted from plasma using the QIAamp Viral RNA Mini Kit (Qiagen). Reverse transcription and PCR amplification of HIV-1 pol gene were done using One-Step SuperScript III RT/Platinum Taq High Fidelity Enzyme Mix (ThermoFisher Scientific), using JA269 and JA272 primers [20]. For nested PCR, High Fidelity Platinum Taq DNA Polymerase, (ThermoFisher Scientific) was used, with primers JA270 and JA271 [20], resulting in a PCR fragment of 1086 bases. The PCR products were sequenced in both directions with six primers described by Zhou et al. [21] using the BigDye terminator kit v 1.1 (Applied Biosystem) followed by sequence analysis on an ABI PRISM 3130 ×l genetic analyzer (Applied Biosystem). Sequence assembly and editing were performed using RECall V 2.0 HIV sequencing analysis tool [22]. The final length of all the sequences following removal of regions corresponding to the primers, editing and alignment was 1035 bases, corresponding to nucleotide positions 2268–3302 of HXB2 (GenBank accession number K03455). Sequence quality control was performed using the quality control program of the Los Alamos HIV sequence database [23]. The presence of drug resistance mutations (DRM) was assessed using the Stanford Genotypic Resistance Interpretation Algorithm [24]. DRM were examined according to the calibrated population resistance (CPR) tool v6.0 beta [25] (http://cpr.stanford.edu/cpr/servlet/CPR), based on the WHO surveillance transmitted drug resistance mutation list of 2009 [26]. Pretreatment DR was referred to as low (< 5%), moderate (5–15%), or high (>15%) [27].

HIV subtyping

All sequences were screened for recombination using RDP v.3.44 [28]. The sequences were subtyped through phylogenetic analysis with group M HIV reference sequences from Los Alamos HIV database [23] with the addition of sub-subtype A3 polymerase sequences [29]. Phylogenetic subtyping was verified using COMET and jpHMM [30, 31]. Sequences were aligned using ClustalX2 and then edited to a final length of 1035 bases for each sequence using BioEdit [32, 33]. A Maximum likelihood phylogenetic tree was constructed using the online version of PhyML with the GTR+I+Γ nucleotide substitution model (using estimated proportion of invariable sites and four gamma categories) and SPR to estimate the tree topology. Branch support was determined with aLRT-SH (approximate likelihood ratio test Shimodaira-Hasegawa like) implemented in PhyML. A branch in the phylogeny with an aLRT-SH value ≥0.9 was considered significant [34, 35]. Cluster analysis was performed using ClusterPicker, where phylogenetic linkage was defined as an aLRT-SH branchsupport ≥0.9 and a genetic distance threshold of 5% [36].

Statistical analysis

Statistical tests were performed using SPSS 24 (IBM Corp., Armonk, NY, USA). Resistance prevalence analyses were performed on 48 patients. Variables are expressed as medians with interquartile ranges (IQRs). Level of DR prevalence was determined with a confidence interval (CI) of 95% using the Wilson method. Univariate analysis for association between resistance mutations and sociodemographic/laboratorial variables was done using 2-tailed Fisher’s exact test.

Ethical approval

The study was approved by the ethical committees of the National Health Ethics Committee in Guinea-Bissau (Ref 038/CNES/INASA/2016) and the Regional Ethical Review Board, Lund University, Sweden (Dnr 2016/426). All participants received information about the study before inclusion, and provided oral and written informed consent. To ensure confidentiality, all study data was managed under code.

Availability of data

Nucleotide sequences reported in this study have been deposited in the Genbank repository (Accession Numbers: MH605452-MH605505).

Results

Study population

Fifty-two pregnant women met the survey inclusion criteria. Among them, 48 were confirmed HIV-1 infected and four were confirmed HIV-2 infected. The HIV-1 pol gene was successfully amplified and sequenced from plasma RNA of all 48 HIV-1 infected individuals. Among the 48 confirmed HIV-1 positive participants, 21 (43.8%) were enrolled at Bairro Militar Health Centre, nine (18.8%) at Antula Health Centre, ten (20.8%) at Quelele Health Centre and eight (16.7%) at Plack-II Health Centre, all located in the capital Bissau. Demographic characteristics of participants are shown in Table 1. Data in regard to absolute CD4 count, CD4%, age and previous pregnancies were missing for two, four, four and three of the women, respectively. The median CD4+ T-cell count was 415 (IQR: 292–562) cells/mL and the median CD4+ T-cell percentage (CD4%) was 19.4 (IQR: 16.1–25.9). The median age was 25 (IQR: 22–28) years and the number of women with previous pregnancies was 22/45 (48.9%).

[Figure omitted. See PDF.]

Table 1. Demographic characteristics of enrolled participants (n = 48).

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

Subtype distribution

Phylogenetic subtyping and verification as described in the methods section identified that the vast majority of study participants (88%; 42 of 48) carried a HIV-1 CRF02_AG polymerase sequence, while sub-subtype A3 (n = 2), subtype G (n = 1), B (n = 1) and CRF06_cpx (n = 1) were found less frequently. One sequence represented a putative recombinant CRF06_xpx/CRF02_AG pol sequence. There were no signs of phylogenetic clustering due to recent transmissions or laboratory contamination, since there were no clusters containing sequences with >95% nucleotide sequence similarity (Fig 1).

[Figure omitted. See PDF.]

Fig 1. Maximum likelihood phylogenetic analysis of HIV-1 pol sequences (specimen names are indicated at tips in black font) from 48 women residing in Bissau, together with HIV-1 reference sequences (in green font), showed that sequences of the study participants belonged to five different subtypes/CRFs (highlighted within a box shaded pink, where the subtype/CRF names are shown to the right).

One sequence (Plk_9) represented a putative recombinant between CRF06_cpx and CRF02_AG (indicated by an asterisk). Red internal branches indicate highly supported branches (aLRT-SH ≥0.90) while green internal branches indicate well-supported branches (0.87 ≤ aLRT-SH <0.90). The prefix of the taxa names indicates the collection site (BM, Bairro Militar Health Centre; A, Antula Health Centre; Q, Quelele Health Centre; and Plk, Plack-II Health Centre). The length of the horizontal bar below the phylogenetic tree corresponds to 0.02 nucleotide substitutions/site.

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

Resistance mutations

Five individuals (5 of 48; 10.4%) carried major DRM. Four carried mutations towards NNRTI only (K103N and K103N/S) and one carried mutations to both NNRTI (G190S, K101E) and NRTI (M184V) No DRM toward protease inhibitors (PIs) was found. This represents a prevalence of DRM to NNRTI, NRTI and PI of 10.4% (95% CI: 4.5–22.2%), 2.1% (95% CI: 0.4–10.9%) and 0% (95% CI: 0.0–7.4%) respectively, corresponding to low levels of DRM to NRTI and PI and moderate levels of DRM to NNRTI. The individuals who presented with pretreatment DR are shown in Table 2. Women ≥25 years were more likely to harbour DRM (p = 0.049, FET; Table 3). Number of previous pregnancies, marital status, education, CD4 absolute count, CD4% and subtype/CRF were not associated with DRMs (Table 3).

[Figure omitted. See PDF.]

Table 2. Study participants with pretreatment drug resistance.

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

[Figure omitted. See PDF.]

Table 3. Medical and socioeconomic data of study participants in relation to pretreatment drug resistance (PDR).

https://doi.org/10.1371/journal.pone.0206406.t003

Discussion

In this study we found 10.4% pretreatment DR in an ART-naïve group of pregnant women in Guinea Bissau. This is the first study of pretreatment HIV-1 DR in the country and the results raise some concerns. First line treatment in Guinea-Bissau, as in most LMIC relies on two NRTI and one NNRTI. In Guinea Bissau, the most commonly used NNRTI is Efavirenz [37] which is less effective in patients with the mutation K103N/S. The K103N/S mutation/s was found most frequently in this study and represents one of the most commonly detected pretreatment DR mutations in most LMIC since ART roll out [38].

We confirm our previous observations that CRF02_AG is the most common HIV form in Guinea-Bissau. In fact, our observation indicated an increased dominance of CRF02_AG in the country, from 57% in 1993–2008 to 88% in this study [13]. Among the specimens sequenced, we found no evidence of phylogenetic clustering, indicating that the study participants were not epidemiologically linked due to recent transmissions. However, further studies are necessary to confirm trends in prevalence of HIV subtypes/CRF in Guinea Bissau and association of clustering, if any. Moreover, the PCR protocol employed in our study was highly successful with a 100% amplification rate, indicating that the PCR protocol is useful for amplifying different HIV-1 subtypes/CRFs (21).

Previous studies in West Africa have found low levels of pretreatment DR in Guinea, moderate levels in Mali and levels varying from low to moderate in Cape Verde [39–43]. Levels of pretreatment DR in pregnant women were low in studies from Ghana and Nigeria, while moderate levels were observed in Burkina Faso [44–46]. Our study, which was performed twelve years after ART roll out in Guinea Bissau, found levels of pretreatment DR higher than the above mentioned studies, a difference that may be associated with the time of the studies after national ART scale up in the countries [8]. Similar observations have been made in other sub-Saharan countries after national ART scale-up and are mainly attributable to an increase in NNRTI resistance [38]. The specific mutations found in this study have also been found among patients failing treatment in Bissau, suggesting a link between acquired and transmitted drug resistance in the country [4, 47, 48].

WHO recommends performing threshold surveys in sites where a large proportion of individuals are likely to be young or recently infected and for these reasons the choice of antenatal clinics are representative investigation sites [19]. Several limitations are however evident in our study. Despite obtaining information of previous ART use and HIV diagnoses, almost half of the participants had previous pregnancies. We found no association between the number of previous pregnancies, marital status, education, CD4 absolute count, CD4% and subtype/CRF with DRMs. However, since the number of participants in the study was low, larger studies are needed to address these associations more accurately. It is possible that some participants concealed information regarding knowledge of HIV status and previous ART/PMTCT use which opens up for the possibility that a detected DR in an individual could be acquired and not transmitted by nature. Although we included participants from different clinics, they were all located in the capital of Bissau which may be less representative in the country as a whole due to regional differences in pretreatment DR [49]. There were also a number of challenges in recruiting participants at antenatal clinics in Bissau due to periods of stock-outs of HIV-antibody tests, limiting a full assessment of resistance in the region. In addition, acute and recent infections were not evaluated which may have led to an underestimation of the pretreatment DR level. Moreover, we included participants who were older than 25 year, while a 25 year age limit is recommended by the WHO. The purpose of this criterion is to limit the number of years of potential exposure to HIV transmission. However, and which has been discussed by Kasang et al. previously, it is questionable whether the age criterion is representative for a “sentinel group”, while the clinical reality is related to pre-existing HIVDR in all patients eligible for ART, irrespective of age and HIVDR origin (i.e. transmitted or acquired HIVDR) [50]. Regardless of whether the level of HIVDR reported here is predominately caused by transmission or selection due to undisclosed ART, it is likely to affect the first-line ART efficacy for a substantial number of individuals in the country.

In conclusion, even moderate levels of pretreatment DR in a presumably naïve population should raise concerns since it may represent a threat to the ART program in the country. As more individuals will start ART, more durable first line regimens are required. To respond to the threat of NNRTI drug resistance, the WHO recommends countries with a pre-treatment NNRTI resistance prevalence of 10% or higher to introduce a non–NNRTI-containing first-line regimen, for example the integrase inhibitor Dolutegravir [51]. Whether Guinea Bissau has reached this critical level should be further investigated. The results of the current study indicate that a substantial proportion of HIV-1 infected individuals are at risk for a suboptimal treatment compromising health and increasing the risk of transmission of drug resistant virus. Our findings provide novel data that may inform the development of more effective public health strategies in HIV prevention in Guinea-Bissau.

Acknowledgments

The listed authors and the members of the Sweden Guinea-Bissau Cohort Research (SWEGUB CORE) group, including Babetida N’Buna, Antonio Biague, Ansu Biai, Cidia Camara, Joakim Esbjörnsson, Marianne Jansson, Sara Karlson, Jacob Lopatko Lindman, Patrik Medstrand, Fredrik Månsson, Hans Norrgren, Gülsen Özkaya Sahin, Zacarias Jose da Silva and Sten Wilhelmson, are indebted to the staff of the antenatal care clinics and the National Public Health Laboratory (LNSP) in Bissau, Guinea-Bissau. We also thank Sara Karlson and Dawit Assefa for technical support.

Citation: Wilhelmson S, Månsson F, Lopatko Lindman J, Biai A, Esbjörnsson J, Norrgren H, et al. (2018) Prevalence of HIV-1 pretreatment drug resistance among treatment naïve pregnant women in Bissau, Guinea Bissau. PLoS ONE 13(10): e0206406. https://doi.org/10.1371/journal.pone.0206406

References

1. World Bank. Guinea-Bissau 2018 [2018, January 16]. https://data.worldbank.org/country/guinea-bissau.

2. UNAIDS. Country factsheets Guinea-Bissau. 2016.

3. World Health Organization. What ART regimen to start with (first-line ART). 2018.

4. Jespersen S, Tolstrup M, Honge BL, Medina C, Te Dda S, Ellermann-Eriksen S, et al. High level of HIV-1 drug resistance among patients with HIV-1 and HIV-1/2 dual infections in Guinea-Bissau. Virol J. 2015;12:41. pmid:25889017

5. World Health Organization. Global report on early warning indicators of HIV drug resistance 2016 [2018, June 28]. http://apps.who.int/iris/bitstream/handle/10665/246219/9789241511179-eng.pdf?sequence=1.

6. Wallis CL, Godfrey C, Fitzgibbon JE, Mellors JW. Key Factors Influencing the Emergence of Human Immunodeficiency Virus Drug Resistance in Low- and Middle-Income Countries. J Infect Dis. 2017;216(suppl_9):S851–S6. pmid:29207000

7. Frentz D, Boucher CA, van de Vijver DA. Temporal changes in the epidemiology of transmission of drug-resistant HIV-1 across the world. AIDS Rev. 2012;14(1):17–27. pmid:22297501.

8. Gupta RK, Jordan MR, Sultan BJ, Hill A, Davis DH, Gregson J, et al. Global trends in antiretroviral resistance in treatment-naive individuals with HIV after rollout of antiretroviral treatment in resource-limited settings: a global collaborative study and meta-regression analysis. Lancet. 2012;380(9849):1250–8. pmid:22828485

9. world health Organization. Global action plan on HIV drug resistance 2017–2021 2017 [2018, June 28]. http://apps.who.int/iris/bitstream/handle/10665/255883/9789241512848-eng.pdf?sequence=1.

10. Hemelaar J, Gouws E, Ghys PD, Osmanov S, Who-Unaids Network for HIV Isolation Characterisation. Global trends in molecular epidemiology of HIV-1 during 2000–2007. AIDS. 2011;25(5):679–89. pmid:21297424

11. Hemelaar J, Gouws E, Ghys PD, Osmanov S. Global and regional distribution of HIV-1 genetic subtypes and recombinants in 2004. AIDS. 2006;20(16):W13–23. pmid:17053344.

12. Palm AA, Esbjornsson J, Mansson F, Kvist A, Isberg PE, Biague A, et al. Faster progression to AIDS and AIDS-related death among seroincident individuals infected with recombinant HIV-1 A3/CRF02_AG compared with sub-subtype A3. J Infect Dis. 2014;209(5):721–8. pmid:23935204.

13. Esbjornsson J, Mild M, Mansson F, Norrgren H, Medstrand P. HIV-1 molecular epidemiology in Guinea-Bissau, West Africa: origin, demography and migrations. PLoS One. 2011;6(2):e17025. pmid:21365013

14. Hu DJ, Vanichseni S, Mastro TD, Raktham S, Young NL, Mock PA, et al. Viral load differences in early infection with two HIV-1 subtypes. AIDS. 2001;15(6):683–91. pmid:11371682.

15. Neilson JR, John GC, Carr JK, Lewis P, Kreiss JK, Jackson S, et al. Subtypes of human immunodeficiency virus type 1 and disease stage among women in Nairobi, Kenya. J Virol. 1999;73(5):4393–403. pmid:10196337

16. Morrison CS, Demers K, Kwok C, Bulime S, Rinaldi A, Munjoma M, et al. Plasma and cervical viral loads among Ugandan and Zimbabwean women during acute and early HIV-1 infection. AIDS. 2010;24(4):573–82. pmid:20154581

17. Renjifo B, Fawzi W, Mwakagile D, Hunter D, Msamanga G, Spiegelman D, et al. Differences in perinatal transmission among human immunodeficiency virus type 1 genotypes. J Hum Virol. 2001;4(1):16–25. pmid:11213929.

18. Wainberg MA, Brenner BG. Role of HIV Subtype Diversity in the Development of Resistance to Antiviral Drugs. Viruses. 2010;2(11):2493–508. pmid:21994627

19. Bennett DE, Myatt M, Bertagnolio S, Sutherland D, Gilks CF. Recommendations for surveillance of transmitted HIV drug resistance in countries scaling up antiretroviral treatment. Antivir Ther. 2008;13 Suppl 2:25–36. pmid:18575189.

20. Murillo W, de Rivera IL, Parham L, Jovel E, Palou E, Karlsson AC, et al. Prevalence of drug resistance and importance of viral load measurements in Honduran HIV-infected patients failing antiretroviral treatment. HIV Med. 2010;11(2):95–103. pmid:19686436.

21. Zhou Z, Wagar N, DeVos JR, Rottinghaus E, Diallo K, Nguyen DB, et al. Optimization of a low cost and broadly sensitive genotyping assay for HIV-1 drug resistance surveillance and monitoring in resource-limited settings. PLoS One. 2011;6(11):e28184. pmid:22132237

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

23. Los Alamos National Database. HIV sequence database. 2017.

24. Stanford University HIV Drug resistance Database. HIVdb Program Genotypic Resistance Interpretation Algorithm. 2018.

25. Gifford RJ, Liu TF, Rhee SY, Kiuchi M, Hue S, Pillay D, et al. The calibrated population resistance tool: standardized genotypic estimation of transmitted HIV-1 drug resistance. Bioinformatics. 2009;25(9):1197–8. pmid:19304876

26. 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(3):e4724. pmid:19266092

27. Myatt M, Bennett DE. A novel sequential sampling technique for the surveillance of transmitted HIV drug resistance by cross-sectional survey for use in low resource settings. Antivir Ther. 2008;13 Suppl 2:37–48. pmid:18575190.

28. Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P. RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics. 2010;26(19):2462–3. pmid:20798170

29. Meloni ST, Kim B, Sankale JL, Hamel DJ, Tovanabutra S, Mboup S, et al. Distinct human immunodeficiency virus type 1 subtype A virus circulating in West Africa: sub-subtype A3. J Virol. 2004;78(22):12438–45. pmid:15507630

30. Schultz AK, Zhang M, Bulla I, Leitner T, Korber B, Morgenstern B, et al. jpHMM: improving the reliability of recombination prediction in HIV-1. Nucleic acids research. 2009;37(Web Server issue):W647–51. pmid:19443440

31. Struck D, Lawyer G, Ternes AM, Schmit JC, Bercoff DP. COMET: adaptive context-based modeling for ultrafast HIV-1 subtype identification. Nucleic acids research. 2014;42(18):e144. pmid:25120265

32. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series. 1999:95–8.

33. Jeanmougin F, Thompson JD, Gouy M, Higgins DG, Gibson TJ. Multiple sequence alignment with Clustal X. Trends in biochemical sciences. 1998;23(10):403–5. pmid:9810230.

34. Esbjornsson J, Mild M, Audelin A, Fonager J, Skar H, Bruun Jorgensen L, et al. HIV-1 transmission between MSM and heterosexuals, and increasing proportions of circulating recombinant forms in the Nordic Countries. Virus Evol. 2016;2(1):vew010. pmid:27774303

35. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59(3):307–21. pmid:20525638.

36. Ragonnet-Cronin M, Hodcroft E, Hue S, Fearnhill E, Delpech V, Brown AJ, et al. Automated analysis of phylogenetic clusters. BMC Bioinformatics. 2013;14:317. pmid:24191891

37. Health GBMo. General directions for prevention and promotion of health—National AIDS Program. 2016.

38. Rhee SY, Blanco JL, Jordan MR, Taylor J, Lemey P, Varghese V, et al. Geographic and temporal trends in the molecular epidemiology and genetic mechanisms of transmitted HIV-1 drug resistance: an individual-patient- and sequence-level meta-analysis. PLoS medicine. 2015;12(4):e1001810. pmid:25849352

39. Alvarez P, Fernandez McPhee C, Prieto L, Martin L, Obiang J, Avedillo P, et al. HIV-1 Variants and Drug Resistance in Pregnant Women from Bata (Equatorial Guinea): 2012–2013. PLoS One. 2016;11(10):e0165333. pmid:27798676

40. de Pina-Araujo II, Guimaraes ML, Bello G, Vicente AC, Morgado MG. Profile of the HIV epidemic in Cape Verde: molecular epidemiology and drug resistance mutations among HIV-1 and HIV-2 infected patients from distinct islands of the archipelago. PLoS One. 2014;9(4):e96201. pmid:24763617

41. Haidara A, Chamberland A, Sylla M, Aboubacrine SA, Cisse M, Traore HA, et al. High level of primary drug resistance in Mali. HIV Med. 2010;11(6):404–11. pmid:20146734.

42. Maiga AI, Fofana DB, Maiga AC, Diallo F, Ait-Arkoub Z, Daou F, et al. Transmitted antiretroviral drug resistance in newly HIV-infected and untreated patients in Segou and Bamako, Mali. AIDS Res Hum Retroviruses. 2013;29(1):182–6. pmid:22823755

43. Oliveira V, Bartolo I, Borrego P, Rocha C, Valadas E, Barreto J, et al. Genetic diversity and drug resistance profiles in HIV type 1- and HIV type 2-infected patients from Cape Verde Islands. AIDS Res Hum Retroviruses. 2012;28(5):510–22. pmid:21902592.

44. Bonney EY, Addo NA, Ntim NA, Addo-Yobo F, Bondzie P, Aryee KE, et al. Low level of transmitted HIV Drug resistance at two HIV care centres in Ghana: a threshold survey. Ghana Med J. 2013;47(2):82–6. pmid:23966745

45. Imade GE, Sagay AS, Chaplin B, Chebu P, Musa J, Okpokwu J, et al. Short communication: Transmitted HIV drug resistance in antiretroviral-naive pregnant women in north central Nigeria. AIDS Res Hum Retroviruses. 2014;30(2):127–33. pmid:24164431

46. Somda A, Sangare L, Soro M, Yameogo S, Bazie B, Bigirimana F, et al. Surveillance of transmitted drug-resistant HIV among young pregnant women in Ouagadougou, Burkina Faso. Clin Infect Dis. 2012;54 Suppl 4:S317–9. pmid:22544195

47. Tilghman MW, Perez-Santiago J, Osorio G, Little SJ, Richman DD, Mathews WC, et al. Community HIV-1 drug resistance is associated with transmitted drug resistance. HIV Med. 2014;15(6):339–46. pmid:24417811

48. Winand R, Theys K, Eusebio M, Aerts J, Camacho RJ, Gomes P, et al. Assessing transmissibility of HIV-1 drug resistance mutations from treated and from drug-naive individuals. AIDS. 2015;29(15):2045–52. pmid:26355575

49. TenoRes Study G. Global epidemiology of drug resistance after failure of WHO recommended first-line regimens for adult HIV-1 infection: a multicentre retrospective cohort study. The Lancet Infectious diseases. 2016;16(5):565–75. pmid:26831472

50. Kasang C, Kalluvya S, Majinge C, Stich A, Bodem J, Kongola G, et al. HIV drug resistance (HIVDR) in antiretroviral therapy-naive patients in Tanzania not eligible for WHO threshold HIVDR survey is dramatically high. PLoS One. 2011;6(8):e23091. pmid:21886779

51. World Health Organization. Guidelines on the Public Health Response to Pretreatment HIV Drug Resistance. 2017.