Skip to main content

HIV-1 subtype C predicted co-receptor tropism in Africa: an individual sequence level meta-analysis

Abstract

Background

Entry inhibitors, such as Maraviroc, hold promise as components of HIV treatment and/or pre-exposure prophylaxis in Africa. Maraviroc inhibits the interaction between HIV Envelope gp120 V3-loop and CCR5 coreceptor. HIV-1 subtype C (HIV-1-C) is predominant in Southern Africa and preferably uses CCR5 co-receptor. Therefore, a significant proportion of HIV-1-C CXCR4 utilizing viruses (X4) may compromise the effectiveness of Maraviroc. This analysis examined coreceptor preferences in early and chronic HIV-1-C infections across Africa.

Methods

African HIV-1-C Envelope gp120 V3-loop sequences sampled from 1988 to 2014 were retrieved from Los Alamos HIV Sequence Database. Sequences from early infections (< 186 days post infection) and chronic infections (> 186 days post infection) were analysed for predicted co-receptor preferences using Geno2Pheno [Coreceptor] 10% FPR, Phenoseq-C, and PSSMsinsi web tools. V3-loop diversity was determined, and viral subtype was confirmed by phylogenetic analysis. National treatment guidelines across Africa were reviewed for Maraviroc recommendation.

Results

Sequences from early (n = 6316) and chronic (n = 7338) HIV-1-C infected individuals from 10 and 15 African countries respectively were available for analyses. Overall, 518/6316 (8.2%; 95% CI 0.7–9.3) of early sequences were X4, with Ethiopia and Malawi having more than 10% each. For chronic infections, 8.3% (95% CI 2.4–16.2) sequences were X4 viruses, with Ethiopia, Tanzania, and Zimbabwe having more than 10% each. For sequences from early chronic infections (< 1 year post infection), the prevalence of X4 viruses was 8.5% (95% CI 2.6–11.2). In late chronic infections (≥ 5 years post infection), X4 viruses were observed in 36% (95% CI − 16.3 to 49.9), with two countries having relatively high X4 viruses: South Africa (43%) and Malawi (24%). The V3-loop amino acid sequence were more variable in X4 viruses in chronic infections compared to acute infections, with South Africa, Ethiopia and Zimbabwe showing the highest levels of V3-loop diversity. All sequences were phylogenetically confirmed as HIV-1-C and clustered according to their co-receptor tropism. In Africa, Maraviroc is registered only in South Africa and Uganda.

Conclusions

Our analyses illustrate that X4 viruses are present in significantly similar proportions in early and early chronic HIV-1 subtype C infected individuals across Africa. In contrast, in late chronic infections, X4 viruses increase 3–5 folds. We can draw two inferences from our observations: (1) to enhance the utility of Maraviroc in chronic HIV subtype C infections in Africa, prior virus co-receptor determination is needed; (2) on the flip side, research on the efficacy of CXCR4 antagonists for HIV-1-C infections is encouraged. Currently, the use of Maraviroc is very limited in Africa.

Background

Data from the Joint United Nations Programme on HIV/AIDS (UNIAIDS) shows that about 38 million people were living with HIV infection at the end of 2018, 68% of this number are in Africa [1]. The UNAIDS 90–90–90 target translates to 90% of all persons to be tested for HIV, 90% of those infected should be on treatment, and 90% of those on treatment should have their viral load suppressed to undetectable levels. In this scheme, treatment is expected to act as a prevention tool, since the chances of transmission is highly reduced when viral load is undetectable [2,3,4,5]. Combination antiretroviral therapy (cART) is the gold standard for the management of HIV infections. In most developing countries, nucleoside and non-nucleoside reverse transcriptase inhibitors are the backbone of all first and second line of antiretroviral regimens in adults, with a boosted protease component for children in first line treatment [6,7,8,9].

The goal of cART is to rapidly reduce HIV viral load to undetectable levels, thereby permitting the reconstitution of immune function as measured by rising levels of CD4+ cell counts. The increasing ease of access to cART across Africa has tremendously reduced morbidity and mortality due to HIV infection [10,11,12,13,14]. However, treatment is not curative, and a significant proportion of individuals will fail first line and second line regimens making them legible for salvage therapy. Maraviroc, an entry inhibitor, is gaining significance as part of treatment regimens in the United States and elsewhere [15,16,17], but there is little documentation of its use in Africa, even as salvage therapy.

Maraviroc inhibits viral entry by prohibiting the interaction between HIV Envelope gp120 V3-loop and the CCR5 co-receptor, following the interaction of Gp120 with the CD4 molecule [18,19,20,21]. High tolerability, safety, and efficacy in viral reduction in both treatment experienced and naïve patients have been demonstrated for Maraviroc, making it a valuable treatment option against HIV/AIDS [22]. It has been reported that amino acid substitutions, particularly by glycine, in the V3 loop crown motive may reduce the binding efficiency of gp120 to CCR5 [23]. Globally, about 47% of infections are due to HIV-1 subtype C (HIV-1-C) [24], and HIV-1-C overwhelmingly dominates infections in Southern Africa [25,26,27]. It is not conclusive how and why HIV-1-C appears to be more transmissible than other members of HIV Group M and how it became the dominant variant in Southern Africa. However, some Ex-vivo pathogenic fitness studies and long-term natural history cohorts in Uganda and Zimbabwe have suggested that subtype C is the least fit subtype among HIV-1 group M. Their lower virulence leads to longer asymptomatic periods which could explain their continuous dominance and expansion in the HIV global pandemic [28, 29]. Further, phylogeography reconstruction models of HIV polymerase sequences [30], have shown that HIV-1-C may have originated from Lubumbashi and Mbuji-Mayi (cities in the Democratic Republic of Congo) from where it was transmitted to Southern Africa, facilitated by the developed and busy road and rail networks in the 1960s. In addition [31], had demonstrated that conserved V1–V2 loops and V3-316T, which occur at higher frequencies in HIV-1-C, increase viral infectivity; and proposed that this could be responsible for the relatively high transmissibility of HIV-1-C heterosexually.

A significant majority of the initial infecting HIV-1-C viruses utilize CCR5. However, the presence of a significant proportion of CXCR4 utilizing viruses (X4) in chronic HIV-1-C infection [32,33,34] might compromise the effectiveness of CCR5 antagonists, such as Maraviroc, when included as components in salvage therapy. In the current analyses, co-receptor preference in early and chronic HIV-1-C infections across Africa was examined, using sequences from the Los Alamos HIV Sequence Database, with the view of understanding the co-receptor preference landscape of the epidemiologically important HIV-1-C across the continent. We also examine literature for indications of active use of Maraviroc in African countries. The association of viral tropism to stage of infection and mode of transmission were examined. We observed that X4 viruses are present in similar proportions in early (less than 6 months post infection) and early chronic (less than 1 year post infection) HIV-1 subtype C infected individuals across Africa. On the contrary, in late chronic infections, there is a significant 3–5 fold increase in X4 viruses. Although there is currently a limited use of Maraviroc across Africa, our findings could be useful in the development of treatment management guidelines in regions where HIV-1 subtype C drives the epidemic.

Methods

Sequence search and extraction

HIV-1 subtype C Gp120 V3-loop Sanger generated sequences were extracted from the Los Alamos HIV Sequence Database (https://www.hiv.lanl.gov/content/index) using the sequence search interface. Firstly, sequences were searched and extracted for each African country based on early and chronic infections. In this study, early infection was defined as the period comprising HIV infection, seroconversion, and recent infection [35] and chronic infection was defined as the period post recent infection. Most of the sequences retrieved were generated when tests to measure HIV RNA, and thus detect acute infections were not available in many African countries. Secondly, extracted sequences were categorized according to the following: route of infection (mother-to-child transmission (MTCT) and heterosexual transmission), and disease progression (slow and rapid progressors). Problematic sequences and those with no data on seroconversion dates were excluded from the analyses.

Co-receptor prediction and sequence analyses

Sequences were classified as early (< 186 days post infection) and chronic infections (> 186 days post infection). Sub-classifications included early chronic infection (186 days to 1 year post infection) and late chronic infections (> 5 years post infection). Sequences were separately analysed for predicted co-receptor preferences using Geno2Pheno [Coreceptor] 10% FPR (https://coreceptor.geno2pheno.org/), Phenoseq-C (http://tools.burnet.edu.au/phenoseq/), and PSSMsinsi (https://indra.mullins.microbiol.washington.edu/webpssm/) web tools. Phenoseq-C and PSSMsinsi are particularly HIV-1-C based co-receptor prediction tools [36,37,38,39]. An inferred concordance of all three tools was used to assign the co-receptor biotype. Gp120 V3-loop diversity was determined for both R5 and X4 viruses by entropy plotting, amino acid length, net charge, N-glycosylation and crown motif examination. The genetic subtype was confirmed by phylogenetic analysis. Figure 1 illustrates an outline of the procedures used in the study.

Fig. 1
figure 1

Flow chart illustrating the study procedure: sequence extraction, categorization, co-receptor prediction and diversity analysis

Active use of Maraviroc in Africa

To estimate the active use of Maraviroc in Africa, national ART guidelines and documents on ARV approvals were retrieved by Google search and reviewed for the recommendation and approval of Maraviroc. Search terms included: national ART guidelines AND [African country], and this was done for all African countries. We also reviewed the recent summary statement of HIV treatment regimens in Africa, 2017 [40].

Statistical analyses

To statistically assess the differences between frequencies, a Fisher’s exact test was conducted; confidence intervals at the 95% threshold were calculated to estimate the interval estimate of a population mean; using an online GraphPad QuickCalcs tool (https://www.graphpad.com/quickcalcs/contingency1.cfm). Significant differences were implied when p-values were < 0.05.

Results

HIV-1 subtype C Gp120 V3-loop co-receptor was predicted, analysed and categorized according to early or chronic infections in general, route of transmission (mother-to-child and heterosexual), and pathogenesis (slow and rapid progressors). In addition, the V3-loop amino acid diversity in terms of amino acid length, entropy, net charge and N-glycosylation sites were further analysed (Fig. 1). Sequences used for these analyses were obtained from the Los Alamos HIV Sequence database, spanned from 1998 through 2014 and originated from 15 out of 54 (27.8%) African countries. These countries were mostly from East and Southern Africa where subtype C is most prevalent. A total of 14,641 HIV-1-C gp120 V3-loop sequences were retrieved. Of these, 987 were excluded due to lack of seroconversion date, affording 13,654 sequences, of which 6316 were from early infections, and 7338 were from chronic infections (Table 1 and Fig. 1).

Table 1 Co-receptor prediction using sequences from African HIV-1 subtype C early (< 6 months of seroconversion) and chronically (> 6 months of seroconversion) infected individuals

Co-receptor biotype prediction

A total of 6316 HIV Gp120 V3-loop Sanger generated sequences from early HIV-1 subtype C infections were available from 10 countries, namely; Botswana, Ethiopia, Kenya, Malawi, Rwanda, Senegal, South Africa, Tanzania, Zambia, and Zimbabwe. For chronic infections, 7338 sequences were available from 15 countries, namely; Botswana, Burundi, DR Congo, Ethiopia, Gabon, Gambia, Guinea Bissau, Kenya, Malawi, Rwanda, South Africa, Tanzania, Uganda, Zambia, and Zimbabwe. Analyses of all early infection sequences showed that 518/6316 (8.2%; 95% CI 0.7–9.3) were of X4 variant while 5798 (91.8%; 95% CI 90.7–99.3) were R5 (Table 1). Ethiopia (3/20; 15.0%) and Malawi (339/2100; 16.1%) had more than 10% of X4 using viruses. For all chronic infections, 612/7338 (8.34%; 95% CI 2.4–16.2) of the sequences were X4-tropic, with four countries, Ethiopia, South Africa, Tanzania, and Zimbabwe each having about 10% or more (Table 1). Overall, there was no difference in the proportion of X4 viruses in early (8.2%) versus all chronic infections (8.3%) (p = 0.8; X2 = 0.064). When sequences, from early chronic infections (> 186 days to < 1 year post-infection) were considered, the prevalence of X4 viruses was 156/1832 (8.5%; 95% CI 2.6–11.2). For chronic infections ≥ 186 day, X4 viruses were present in 64/179 (36%; 95% CI − 16.3 to 49.9) of the study population (p = 0.0001; X2 = 65.257), with two countries having relatively high X4 prevalence: South Africa (43%; 58/136) and Malawi (24%, 6/25).

Co-receptor biotype prediction according to mode of transmission

Of the 6316 acute HIV infection sequences, 1052 were from cases of MTCT. Of these, only 46/1052 (4.4%; 95% CI − 1.5 to 6.5) were X4-tropic. Country-wise, Malawi and Tanzania marginally had more X4 viruses (6% and 6.1% respectively) in early infections (Table 2). In comparison, among chronic infections, 143/2115 (6.8%; 95% CI − 10.9 to 31.5) were X4-tropic viruses, with South Africa having a significantly higher number of X4 viruses (40.4%; p = 0.001; X2 = 30.288). The proportion of X4 viruses in early MTCT sequences was not significantly different from chronic MTCT cases; 4.4% (46/1052) versus 6.7% (143/2115) (p = 0.535; X2 = 0.385) respectively (Table 2).

Table 2 Co-receptor prediction using sequences from HIV-1 subtype C mother-to-child and heterosexual transmissions

A total of 5769 sequences from both early and chronic heterosexual transmissions were analysed for X4 tropism. Overall, 385/3289 (11.7%; 95% CI − 1.8 to 13.9) of early infections were X4-tropic, versus 231/2480 (9.3%; 95% CI − 0.9 to 13.5) for chronic infections (p = 0.6446; X2 = 0.213). Country-wise, Malawi had the most X4 viruses in early heterosexual transmissions (28.5%; p = 0.0047; X2 = 8.000), and Tanzania had the most X4 viruses in chronic heterosexual transmissions (32.4%; p = 0.0001; X2 = 30.042), both of which were significant (Table 2).

Co-receptor biotype prediction according to disease progression

Sequences from slow and rapid progressors from South Africa and Zambia were also available in significant numbers for analyses. There were no X4 tropic variants among the 145 early infection sequences from slow progressors. Among the 158 sequences from chronic slow progressors, 2/158 (1.3%; 95% CI 0–0) were X4 tropic, all from South Africa. The rapid progressors had a total of 169 early infections sequences none of which was X4; while 55/141 (39%) of the chronic sequences, were X4, from South Africa. A significantly higher proportion of sequences from chronic rapid progressors were of X4 tropic compared to those from chronic slow progressors (p = 0.001; X2 = 45.995) (Table 3).

Table 3 Co-receptor prediction using sequences from HIV-1 subtype C slow and rapid progressors during early and chronic infections

HIV-1 subtype C gp120 V3-loop diversity in Africa

In order to perform V3-loop diversity analyses, sequences were available for seven countries, namely; Botswana, Ethiopia, Malawi, Tanzania, South Africa, Zambia, and Zimbabwe (Table 4). The least amino acid substitutions (measured by entropy) was observed in X4 and R5 sequences for both early and chronic infections from Botswana ranging from 0–1 and 0–1.25 respectively; while Ethiopia and Malawi had sequences with the highest entropy (range 0–1.915) in both early and chronic sequences of X4 and R5 variants. High levels of entropy were seen in positions 12, 25, 29 and 32 for both the R5 and X4 viruses (Fig. 2).

Table 4 HIV-1 subtype C V3 loop N-glycosylation site and V3 loop crown motif variation in African countries
Fig. 2
figure 2figure 2

Inter-country variation within the V3 region of subtype HIV-1-C R5 and X4 viruses in early and chronic infection. a–d Are entropy plots for variations at each amino acid position for the different African countries

N-glycosylation sites play key roles in the interaction between the virus, CD4 molecule and the CCR5 and CXCR4 co-receptors; as well as aiding the virus to evade neutralization by host immune response. R5 and X4 viruses of early and chronic infections from Tanzania had a highly conserved N-glycosylation site while Botswana, Ethiopia, Malawi, South Africa, and Zambia had very few sequences that lost the N-glycosylation site (Table 4). On the contrary, sequences from Ethiopia showed the highest frequency of X4 viruses that had lost the N-glycosylation site for acute (67%) and chronic (46.9%) infections (Table 4, Fig. 3).

Fig. 3
figure 3

Amino acid alignment of 126 representative sequences of early and chronic R5 and X4 infections. Highlighted areas (blue) indicate the N-glycosylation site and the V3 loop crown motive

The crown motive (GPGQ) was very much conserved, except for Botswana for which about 77% of sequences from X4 viruses of early infections had Q → R (40.3%) and Q → K (35.8%) substitutions; while about 99% of sequences from chronic infections had Q → R (77%) and Q → K (22.6%) substitutions among the X4 viruses (Table 4). The amino acid lengths were shorter for early than chronic R5 or X4 viruses ranging between 34 and 35 for early and 34–38 for chronic infections. There was no particular trend in the V3-loop amino acid length across the countries (Table 5). The V3-loop net charge for both early and chronic X4 viruses were generally higher, ranging from 1 to 10, and lower in R5 viruses which ranged from 1 to 6 (Table 5).

Table 5 V3 loop amino acid diversity (amino acid length, deletions, insertions, net charge) of X4 and R5 viruses in early and chronic infections

As expected, a higher level of gp120 V3-loop amino acid variation was observed in X4 tropic viruses from chronic than early infections, with South Africa, Ethiopia and Zimbabwe showing the highest levels of V3-loop diversity (Fig. 2). All sequences were phylogenetically confirmed as HIV-1 subtype C. More than 92% of the sequences clustered according to their tropism. Only two R5 sequences (3%) and eight X4 sequences (13%) did not cluster according to their tropism based on the selection of sequences used for the phylogenetic analysis (Fig. 4).

Fig. 4
figure 4

Neigbour joining phylogenetic tree of 126 representative sequences of early and chronic R5 and X4 infections from 6 African countries. Majority of these sequences cluster based on predicted biotypes

Use of Maraviroc in Africa

A scoping review of literature on the active use of Maraviroc revealed Tanzania and Zambia as the only countries in Africa that include Maraviroc as a component of salvage therapy in their national ART guidelines, following an HIV tropism test. Although the Southern African HIV Clinicians Society recommends the use of Maraviroc in salvage therapy, this has not been incorporated in the treatment guidelines by the Health Ministries of Southern African Countries. In addition, Maraviroc is registered in South Africa and Uganda. Table 6 shows details for countries on the registration and or recommendation on the use of Maraviroc.

Table 6 Current registration and recommended use of Maraviroc in African countries

Discussion

Maraviroc is a CCR5 antagonist, that prevents HIV from utilizing the CCR5 co-receptor to enter target cells. During the acute phase of infection, HIV strains irrespective of genotype, utilize CCR5 as the main co-receptor. However, as disease progresses, CXCR4 utilizing viruses emerge in about 50% of infected individual [41,42,43]. By and large, in Africa Maraviroc is prescribed mostly for patients who have failed first and second line regimens, which are comprised of nucleoside reverse transcriptase inhibitors, non-nucleoside reverse transcriptase inhibitors, and protease inhibitors [44, 45]. Since, HIV-1-C drives the epidemic in Southern Africa and accounts for about 46% of infections worldwide [24], the current investigation was aimed at determining the distribution of R5 and X4 viruses in early versus chronic infections in HIV-1-C acquired through MTCT and heterosexual routes in Africa.

Using predicted biotypes of about 14,641 gp120 V3-loop sequences, filtered from the Los Alamos HIV database, our analyses showed that X4 variants are present in significantly similar proportions in early and early chronic (< 1 year post-infection) HIV-1-C infected individuals. However, in late chronic infections (5 years post infection), X4 variants increase 3–5 folds. A study by [46], studying paired RNA and proviral DNA from HIV-1-B antiretroviral naive patients with acute and chronic infections, showed that more than 90% of viruses in acute infections from plasma and peripheral blood mononuclear cells were R5, while patients with chronic infections had a significantly higher prevalence of X4 viruses than in patients with acute infection. Another study [47], evaluating 200 patients for co-receptor tropism, using an ultra-deep sequencing approach, also showed that both X4 and R5 viruses co-exists during acute infections, but with R5 viruses as the highly predominant variant. This shift from R5 to X4 viruses was also observed in the current study.

Several hypothesis have been advanced to explain the predominance of R5 viruses in early infection and X4 viruses in chronic infections: that X4 and R5 viruses are transmitted at the same time but X4 viruses are suppressed by the prevailing strong immune response at the time of infection and proliferate later in infection when the immune surveillance is weak [48,49,50]; that X4 viruses emerge from R5 viruses at the beginning of infection [51, 52]; and thirdly, X4 and R5 viruses have different target cell types, with more target cell types for X4 viruses (T-cells) increasing in abundance as infection becomes chronic [53, 54]. A closer look at the sequences in the current study showed that overall, the prevalence of X4 viruses in early infections and early chronic (between 6 months and 1 year post-infection) are similar. However, there is a dramatic rise in the frequency of X4 viruses when sequences from patients with more than 5 years of infection were considered. The inference is that within 4 years of infection with HIV-1-C, the proportion of X4 viruses becomes significantly higher ranging from 24 to 43%. In fact in a recent study, we reported a high frequency (43%) of X4 viruses identified by ultra-deep sequencing from a cohort of HIV-1-C chronically infected individuals from northern South Africa [34]. Overall in the current study, the prevalence of X4 viruses was not significantly different in early and chronic infections among individuals who were vertical infected, and was also similar among individuals who acquired infection heterosexually. It also appears the route of infection does not influence coreceptor tropism in early or chronic infections. We did not find X4 viruses in early infections from both slow and rapid progressors; but there was significantly more X4 viruses in chronic rapid progressors than in chronic slow progressors.

The findings from the current study should be examined in the context of several limitations. Firstly, there was a dearth of sequence data from longitudinal cohorts in Africa to assess evolution of co-receptor usage over time. Nevertheless, in seeking correlates of co-receptor switch between CCR5 and CXCR4 [55], reported no significant difference in co-receptor ‘switching’ over time among patients who were initially infected exclusively with R5 or X4 viruses, when age, viral load, and gender was considered. In the same vein [56], showed that viruses are either R5 or X4 but not dual tropic, and that dual tropism is due to mixture of both phenotypes. Secondly, data meeting the selection criteria for known acute infections and routes of transmission on HIV-1-C infections were available for only 15 African countries, with four countries (Botswana, Malawi, South Africa and Zambia) providing a highly disproportionate number of sequences. This limits the scope of applicability of the findings at least in Southern Africa where HIV-1-C predominantly drives the epidemic. Thirdly, due to the degree of false prediction of co-receptor usage by bioinformatic tools such as Geno2Pheno and position-specific scoring matrices, opinion on the clinical usefulness of predicted coreceptor usage varies [57,58,59]; so predictions will need to be seen in the context of other clinical parameters. Nevertheless, with the observation that X4 viruses exist in an appreciable proportion in chronic infections and that Maraviroc usage as salvage therapy might lead to resistance due to pre-selection of X4 strains [60], it is intriguing why Maraviroc should be reserved for management at the late stage of infection even in those few African countries in which it is recommended for salvage therapy. Nevertheless, new changes in treatment regimens are being introduced. In an effort to reduce drug resistance to non-nucleoside reverse transcriptase inhibitors (NNRTI), many low and middle income countries, and also high income countries are replacing NNRTI with dolutegravir, an integrase strand transfer inhibitors (INSTI), in first and second line treatment regimens [61, 62]. For example, recently in November 2019, South Africa switched from a fixed dose combination of standard tenofovir–lamivudine–nevirapine to a fixed dose combination of tenofovir–lamivudine–dolutegravir. This move may potentially delay, across the board, the use of Maraviroc in patient management.

Conclusion

Our data show that the use of Maraviroc is very limited in Africa, and confirms that for an improved utility of Maraviroc as salvage therapy among HIV-1-C patients in Africa, preliminary virus co-receptor determination is required. Alternatively, Maraviroc may be included as first line therapy in combination with nucleoside analogues; although this may not be beneficial if prevention of mother-to-child transmission is a desirable outcome, since there is no evidence of Maraviroc’s efficacy in the prevention of HIV mother-to-child transmission. Finally, research in CXCR4 antagonists is encouraged as universal access to treatment gains steam across Africa.

Availability of data and materials

Sequences used for this analysis were obtained from the Los Alamos HIV Sequence Database. The extracted dataset is available upon request.

Abbreviations

cART:

Combination antiretroviral therapy

CCR5:

C-C chemokine receptor type 5

CXCR4:

C-X-C chemokine receptor type 4

HIV-1-C:

HIV subtype C

R5:

CCR5 utilizing viruses

MTCT:

Mother to child transmission

X4:

CXC4 utilizing viruses

UNIAIDS:

Joint United Nations Programme on HIV/AIDS

References

  1. UNAIDS. Fact sheet—global AIDS update 2019. UNAIDS. 2019. https://www.unaids.org/en/resources/fact-sheet. Accessed Aug 2019.

  2. Attia S, Egger M, Müller M, Zwahlen M, Low N. Sexual transmission of HIV according to viral load and antiretroviral therapy: systematic review and meta-analysis. AIDS. 2009. https://0-doi-org.brum.beds.ac.uk/10.1097/QAD.0b013e32832b7dca.

    Article  PubMed  Google Scholar 

  3. Chagomerana MB, Miller WC, Tang JH, Hoffman IF, Harrington BJ, DiPrete B, et al. Prevalence of antiretroviral therapy treatment failure among HIV-infected pregnant women at first antenatal care: PMTCT option b+ in Malawi. PLoS ONE. 2018. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0209052.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Myer L, Phillips TK, McIntyre JA, Hsiao NY, Petro G, Zerbe A, et al. HIV viraemia and mother-to-child transmission risk after antiretroviral therapy initiation in pregnancy in Cape Town, South Africa. HIV Med. 2017. https://0-doi-org.brum.beds.ac.uk/10.1111/hiv.12397.

    Article  PubMed  Google Scholar 

  5. Yotebieng M, Mpody C, Ravelomanana NL, Tabala M, Malongo F, Kawende B, et al. HIV viral suppression among pregnant and breastfeeding women in routine care in the Kinshasa province: a baseline evaluation of participants in CQI-PMTCT study. J Int AIDS Soc. 2019. https://0-doi-org.brum.beds.ac.uk/10.1002/jia2.25376.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Murphy RA, Sunpath H, Castilla C, Ebrahim S, Court R, Nguyen H, et al. Second-line antiretroviral therapy: long-term outcomes in South Africa. J Acquir Immune Defic Syndr. 2012. https://0-doi-org.brum.beds.ac.uk/10.1097/QAI.0b013e3182615ad1.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Keiser O. Switching to second-line antiretroviral therapy in resource-limited settings: comparison of programmes with and without viral load monitoring. AIDS. 2009. https://0-doi-org.brum.beds.ac.uk/10.1097/QAD.0b013e32832e05b2.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Kanters S, Socias ME, Paton NI, Vitoria M, Doherty M, Ayers D, et al. Comparative efficacy and safety of second-line antiretroviral therapy for treatment of HIV/AIDS: a systematic review and network meta-analysis. Lancet HIV. 2017. https://0-doi-org.brum.beds.ac.uk/10.1016/S2352-3018(17)30109-1.

    Article  PubMed  Google Scholar 

  9. Ciaffi L, Koulla-Shiro S, Sawadogo AB, Ndour CT, Eymard-Duvernay S, Mbouyap PR, et al. Boosted protease inhibitor monotherapy versus boosted protease inhibitor plus lamivudine dual therapy as second-line maintenance treatment for HIV-1-infected patients in sub-Saharan Africa (ANRS12 286/MOBIDIP): a multicentre, randomised, parallel, open-la. Lancet HIV. 2017. https://0-doi-org.brum.beds.ac.uk/10.1016/S2352-3018(17)30069-3.

    Article  PubMed  Google Scholar 

  10. Adejumo OA, Malee KM, Ryscavage P, Hunter SJ, Taiwo BO. Contemporary issues on the epidemiology and antiretroviral adherence of HIV-infected adolescents in sub-Saharan Africa: a narrative review. J Int AIDS Soc. 2015. https://0-doi-org.brum.beds.ac.uk/10.7448/IAS.18.1.20049.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Rosen S, Fox MP, Gill CJ. Patient retention in antiretroviral therapy programs in sub-Saharan Africa: a systematic review. PLoS Med. 2007. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pmed.0040298.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Eaton JW, Johnson LF, Salomon JA, Bärnighausen T, Bendavid E, Bershteyn A, et al. HIV treatment as prevention: systematic comparison of mathematical models of the potential impact of antiretroviral therapy on HIV incidence in South Africa. PLoS Med. 2012. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pmed.1001245.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Granich R, Gupta S, Hersh B, Williams B, Montaner J, Young B, et al. Trends in AIDS deaths, new infections and ART coverage in the top 30 countries with the highest AIDS mortality burden; 1990–2013. PLoS ONE. 2015. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0131353.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Scanlon ML, Vreeman RC. Current strategies for improving access and adherence to antiretroviral therapies in resource-limited settings. HIV/AIDS Res Palliat Care. 2013. https://0-doi-org.brum.beds.ac.uk/10.2147/HIV.S28912.

    Article  Google Scholar 

  15. Gates TM, Cysique LA, Siefried KJ, Chaganti J, Moffat KJ, Brew BJ. Maraviroc-intensified combined antiretroviral therapy improves cognition in virally suppressed HIV-associated neurocognitive disorder. AIDS. 2016. https://0-doi-org.brum.beds.ac.uk/10.1097/QAD.0000000000000951.

    Article  PubMed  Google Scholar 

  16. Pozo-Balado MM, Rosado-Sánchez I, Méndez-Lagares G, Rodríguez-Méndez MM, Ruiz-Mateos E, Benhnia MR, et al. Maraviroc contributes to the restoration of the homeostasis of regulatory T-cell subsets in antiretroviral-naive HIV-infected subjects. Clin Microbiol Infect. 2016. https://0-doi-org.brum.beds.ac.uk/10.1016/j.cmi.2015.12.025.

    Article  PubMed  Google Scholar 

  17. Williams DW, Li Y, Dastgheyb R, Fitzgerald KC, Maki PM, Spence AB, et al. Associations between antiretroviral drugs on depressive symptomatology in homogenous subgroups of women with HIV. J Neuroimmune Pharmacol. 2020. https://0-doi-org.brum.beds.ac.uk/10.1007/s11481-019-09899-2.

    Article  PubMed  Google Scholar 

  18. Borm K, Jakobsen MR, Cashin K, Flynn JK, Ellenberg P, Ostergaard L, et al. Frequency and Env determinants of HIV-1 subtype C strains from antiretroviral therapy-naive subjects that display incomplete inhibition by maraviroc. Retrovirology. 2016. https://0-doi-org.brum.beds.ac.uk/10.1186/s12977-016-0309-2.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Flynn JK, Ellenberg P, Duncan R, Ellett A, Zhou J, Sterjovski J, et al. Analysis of clinical HIV-1 strains with resistance to maraviroc reveals strain-specific resistance mutations, variable degrees of resistance, and minimal cross-resistance to other CCR5 antagonists. AIDS Res Hum Retroviruses. 2017. https://0-doi-org.brum.beds.ac.uk/10.1089/aid.2017.0097.

    Article  PubMed  Google Scholar 

  20. Shaik MM, Peng H, Lu J, Rits-Volloch S, Xu C, Liao M, et al. Structural basis of coreceptor recognition by HIV-1 envelope spike. Nature. 2019. https://0-doi-org.brum.beds.ac.uk/10.1038/s41586-018-0804-9.

    Article  PubMed  Google Scholar 

  21. Garcia-Perez J, Rueda P, Alcami J, Rognan D, Arenzana-Seisdedos F, Lagane B, et al. Allosteric model of maraviroc binding to CC Chemokine Receptor 5 (CCR5). J Biol Chem. 2011;286(38):33409–21.

    Article  CAS  Google Scholar 

  22. Van Der Ryst E. Maraviroc—a CCR5 antagonist for the treatment of HIV-1 infection. Front Immunol. 2015;6:277.

    Google Scholar 

  23. Garcia-Perez J, Staropoli I, Azoulay S, Heinrich JT, Cascajero A, Colin P, et al. A single-residue change in the HIV-1 V3 loop associated with maraviroc resistance impairs CCR5 binding affinity while increasing replicative capacity. Retrovirology. 2015. https://0-doi-org.brum.beds.ac.uk/10.1186/s12977-015-0177-1.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Hemelaar J, Elangovan R, Yun J, Dickson-Tetteh L, Fleminger I, Kirtley S, et al. Global and regional molecular epidemiology of HIV-1, 1990–2015: a systematic review, global survey, and trend analysis. Lancet Infect Dis. 2019. https://0-doi-org.brum.beds.ac.uk/10.1016/S1473-3099(18)30647-9.

    Article  PubMed  Google Scholar 

  25. Lihana RW, Ssemwanga D, Abimiku A, Ndembi N. Update on HIV-1 diversity in Africa: a decade in review. AIDS Rev. 2012;14:83–100.

    PubMed  Google Scholar 

  26. Kahle E, Campbell M, Lingappa J, Donnell D, Celum C, Wondondo C, et al. HIV-1 subtype C is not associated with higher risk of heterosexual HIV-1 transmission: a multinational study among HIV-1 serodiscordant couples. AIDS. 2014. https://0-doi-org.brum.beds.ac.uk/10.1097/qad.0000000000000024.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Wilkinson E, Engelbrecht S, De Oliveira T. History and origin of the HIV-1 subtype C epidemic in South Africa and the greater southern African region. Sci Rep. 2015. https://0-doi-org.brum.beds.ac.uk/10.1038/srep16897.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Abraha A, Nankya IL, Gibson R, Demers K, Tebit DM, Johnston E, et al. CCR5- and CXCR4-tropic subtype C human immunodeficiency virus type 1 isolates have a lower level of pathogenic fitness than other dominant group M subtypes: implications for the epidemic. J Virol. 2009. https://0-doi-org.brum.beds.ac.uk/10.1128/jvi.02051-08.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Venner CM, Nankya I, Kyeyune F, Demers K, Kwok C, Chen P-L, et al. Infecting HIV-1 subtype predicts disease progression in women of sub-Saharan Africa. EBioMedicine. 2016;13:305–14.

    Article  Google Scholar 

  30. Faria RN, Vidal N, Lourenco J, Raghwani J, Sigaloff KCE, Tatem JA, et al. Distinct rates and patterns of spread of the major HIV-1 subtypes in Central and East Africa. PLoS Pathog. 2019;15(12):e1007976.

    Article  Google Scholar 

  31. Walter BL, Armitage AE, Graham SC, de Oliveira T, Skinhøj P, Jones EY, et al. Functional characteristics of HIV-1 subtype-C compatible with increased heterosexual transmissibility. AIDS. 2009;23(9):1047–57.

    Article  CAS  Google Scholar 

  32. Ketseoglou I, Lukhwareni A, Steegen K, Carmona S, Stevens WS, Papathanasopoulos MA. Viral tropism and antiretroviral drug resistance in HIV-1 subtype C-infected patients failing highly active antiretroviral therapy in Johannesburg, South Africa. AIDS Res Hum Retroviruses. 2014;30(3):289–93.

    Article  CAS  Google Scholar 

  33. Sollerkvist LP, Gaseitsiwe S, Mine M, Sebetso G, Mphoyakgosi T, Diphoko T, et al. Increased CXCR4 use of HIV-1 subtype C identified by population sequencing in patients failing antiretroviral treatment compared with treatment-naive patients in botswana. AIDS Res Hum Retroviruses. 2014. https://0-doi-org.brum.beds.ac.uk/10.1089/aid.2013.0203.

    Article  PubMed  Google Scholar 

  34. Matume ND, Tebit DM, Gray LR, Hammarskjold M, Rekosh D, Bessong PO. Next generation sequencing reveals a high frequency of CXCR4 utilizing viruses in HIV-1 chronically infected drug experienced individuals in South Africa. J Clin Virol. 2018. https://0-doi-org.brum.beds.ac.uk/10.1016/j.jcv.2018.02.008.

    Article  PubMed  Google Scholar 

  35. Fiebig EW, Wright DJ, Rawal BD, Garrett PE, Schumacher RT, Peddada L, et al. Dynamics of HIV viremia and antibody seroconversion in plasma donors: implications for diagnosis and staging of primary HIV infection. AIDS. 2003. https://0-doi-org.brum.beds.ac.uk/10.1097/00002030-200309050-00005.

    Article  PubMed  Google Scholar 

  36. Jensen MA, Coetzer M, van’t Wout AB, Morris L, Mullins JI. A reliable phenotype predictor for human immunodeficiency virus type 1 subtype C based on envelope V3 sequences. J Virol. 2006;80(10):4698–704.

    Article  CAS  Google Scholar 

  37. Jensen MA, Coetzer M, van’t Wout AB, Morris L, Mullins JI. A reliable phenotype predictor for human immunodeficiency virus type 1 subtype C based on envelope V3 sequences. J Virol. 2006. https://0-doi-org.brum.beds.ac.uk/10.1128/jvi.80.10.4698-4704.2006.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Riemenschneider M, Cashin KY, Budeus B, Sierra S, Shirvani-Dastgerdi E, Bayanolhagh S, et al. Genotypic prediction of co-receptor tropism of HIV-1 subtypes A and C. Sci Rep. 2016. https://0-doi-org.brum.beds.ac.uk/10.1038/srep24883.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Cashin K, Jakobsen MR, Sterjovski J, Roche M, Ellett A, Flynn JK, et al. Linkages between HIV-1 specificity for CCR5 or CXCR4 and in vitro usage of alternative coreceptors during progressive HIV-1 subtype C infection. Retrovirology. 2013;10(1):98.

    Article  Google Scholar 

  40. Strengthening High Impact Intervention for an AIDS-free Generation (AIDSFree) Project. Summary table of HIV treatment Regimens. 2017. https://aidsfree.usaid.gov/sites/default/files/2018.2.26_hiv-xwalk-update.pdf. Accessed 19 Jan 2020.

  41. Mild M, Kvist A, Esbjörnsson J, Karlsson I, Fenyö EM, Medstrand Patrik P. Differences in molecular evolution between switch (R5 to R5X4/X4-tropic) and non-switch (R5-tropic only) HIV-1 populations during infection. Infect Genet Evol. 2010. https://0-doi-org.brum.beds.ac.uk/10.1016/j.meegid.2009.05.003.

    Article  PubMed  Google Scholar 

  42. Verhofstede C, Nijhuis M, Vandekerckhove L. Correlation of coreceptor usage and disease progression. Curr Opin HIV AIDS. 2012. https://0-doi-org.brum.beds.ac.uk/10.1097/COH.0b013e328356f6f2.

    Article  PubMed  Google Scholar 

  43. Philpott S. HIV-1 coreceptor usage, transmission, and disease progression. Curr HIV Res. 2005. https://0-doi-org.brum.beds.ac.uk/10.2174/1570162033485357.

    Article  Google Scholar 

  44. Palladino C, Gómez MLN, Soler-Palacín P, González-Tomé MI, De Ory SJ, Espiau M, et al. Off-label use of maraviroc in HIV-1-infected paediatric patients in clinical practice. AIDS. 2015. https://0-doi-org.brum.beds.ac.uk/10.1097/QAD.0000000000000819.

    Article  PubMed  Google Scholar 

  45. Meini G, Rossetti B, Bianco C, Ceccherini-Silberstein F, Di Giambenedetto S, Sighinolfi L, et al. Longitudinal analysis of HIV-1 coreceptor tropism by single and triplicate HIV-1 RNA and DNA sequencing in patients undergoing successful first-line antiretroviral therapy. J Antimicrob Chemother. 2014. https://0-doi-org.brum.beds.ac.uk/10.1093/jac/dkt426.

    Article  PubMed  Google Scholar 

  46. Bon I, Turriziani O, Musumeci G, Clò A, Montagna C, Morini S, et al. HIV-1 coreceptor usage in paired plasma RNA and proviral DNA from patients with acute and chronic infection never treated with antiretroviral therapy. J Med Virol. 2015. https://0-doi-org.brum.beds.ac.uk/10.1002/jmv.24036.

    Article  PubMed  Google Scholar 

  47. Raymond S, Saliou A, Nicot F, Delobel P, Dubois M, Carcenac R, et al. Characterization of CXCR4-using HIV-1 during primary infection by ultra-deep pyrosequencing. J Antimicrob Chemother. 2013. https://0-doi-org.brum.beds.ac.uk/10.1093/jac/dkt290.

    Article  PubMed  Google Scholar 

  48. Wodarz D, Lloyd AL, Jansen VAA, Nowak MA. Dynamics of macrophage and T cell infection by HIV. J Theor Biol. 1999. https://0-doi-org.brum.beds.ac.uk/10.1006/jtbi.1998.0816.

    Article  PubMed  Google Scholar 

  49. Wodarz D, Nowak MA. The effect of different immune responses on the evolution of virulent CXCR4-tropic HIV. Proc R Soc B Biol Sci. 1998. https://0-doi-org.brum.beds.ac.uk/10.1098/rspb.1998.0552.

    Article  Google Scholar 

  50. Bou-Habib DC, Roderiquez G, Oravecz T, Norcross MA, Berman PW, Lusso P. Cryptic nature of envelope V3 region epitopes protects primary monocytotropic human immunodeficiency virus type 1 from antibody neutralization. J Virol. 1994;68(9):600–6013.

    Article  Google Scholar 

  51. Bunnik EM, Swenson LC, Edo-Matas D, Huang W, Dong W, Frantzell A, et al. Detection of inferred CCR5- and CXCR4-using HIV-1 variants and evolutionary intermediates using ultra-deep pyrosequencing. PLoS Pathog. 2011. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.ppat.1002106.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Poon AFY, Swenson LC, Bunnik EM, Edo-Matas D, Schuitemaker H, van’t Wout AB, et al. Reconstructing the dynamics of HIV evolution within hosts from serial deep sequence data. PLoS Comput Biol. 2012. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pcbi.1002753.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Bleul CC, Wu L, Hoxie JA, Springer TA, Mackay CR. The HIV coreceptors CXCR4 and CCR5 are differentially expressed and regulated on human T lymphocytes. Proc Natl Acad Sci USA. 1997. https://0-doi-org.brum.beds.ac.uk/10.1073/pnas.94.5.1925.

    Article  PubMed  Google Scholar 

  54. Blanpain C, Migeotte I, Lee B, Vakili J, Doranz BJ, Govaerts C, et al. CCR5 binds multiple CC-chemokines: MCP-3 acts as a natural antagonist. Blood. 1999;94(6):1899–905.

    Article  CAS  Google Scholar 

  55. Arif MS, Hunter J, Léda AR, Zukurov JPL, Samer S, Camargo M, et al. Pace of coreceptor tropism switch in HIV-1-infected individuals after recent infection. J Virol. 2017. https://0-doi-org.brum.beds.ac.uk/10.1128/jvi.00793-17.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Singh A, Page T, Moore PL, Allgaier RL, Hiramen K, Coovadia HM, et al. Functional and genetic analysis of coreceptor usage by dualtropic HIV-1 subtype C isolates. Virology. 2009;393(1):56–67.

    Article  CAS  Google Scholar 

  57. Low AJ, Dong W, Chan D, Sing T, Swanstrom R, Jensen M, et al. Current V3 genotyping algorithms are inadequate for predicting X4 co-receptor usage in clinical isolates. AIDS. 2007. https://0-doi-org.brum.beds.ac.uk/10.1097/QAD.0b013e3282ef81ea.

    Article  PubMed  Google Scholar 

  58. Delobel P, Nugeyre MT, Cazabat M, Pasquier C, Marchou B, Massip P, et al. Population-based sequencing of the V3 region of env for predicting the coreceptor usage of human immunodeficiency virus type 1 quasispecies. J Clin Microbiol. 2007. https://0-doi-org.brum.beds.ac.uk/10.1128/JCM.02090-06.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Poveda E, Briz V, Roulet V, Del Mar González M, Faudon JL, Skrabal K, et al. Correlation between a phenotypic assay and three bioinformatic tools for determining HIV co-receptor use. AIDS. 2007. https://0-doi-org.brum.beds.ac.uk/10.1097/QAD.0b013e32826fb741.

    Article  PubMed  Google Scholar 

  60. Westby M, Lewis M, Whitcomb J, Youle M, Pozniak AL, James IT, et al. Emergence of CXCR4-using human immunodeficiency virus type 1 (HIV-1) variants in a minority of HIV-1-infected patients following treatment with the CCR5 antagonist maraviroc is from a pretreatment CXCR4-using virus reservoir. J Virol. 2006. https://0-doi-org.brum.beds.ac.uk/10.1128/jvi.80.10.4909-4920.2006.

    Article  PubMed  PubMed Central  Google Scholar 

  61. 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. https://0-doi-org.brum.beds.ac.uk/10.1016/s1473-3099(18)30710-2.

    Article  PubMed  Google Scholar 

  62. Nance RM, Vannappagari V, Smith K, Johannes CB, Calingaert B, Saltus CW, et al. Virologic failure among people living with HIV initiating dolutegravir-based versus other recommended regimens in real-world clinical care settings. J Acquir Immune Defic Syndr. 2019;81(5):572–7.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We thank Joachim Buch of Max Planck Institute for Informatics, Germany, for assisting in coreceptor predication using the Pheno2Geno tool.

Funding

Research reported in this publication was support by the South African Medical Research Council (RCDI) through funding received from the South African National Treasury; additional support was obtained from the South African National Research Foundation (GUN109312). Nontokozo D. Matume was supported by the Research Capacity Development Initiative of the Medical Research Council (RCDI Project number: 57009). The contents are solely the responsibility of the authors and do not necessarily represent the official views of the University of Virginia, the South African Medical Research Council or the National Research Foundation.

Author information

Authors and Affiliations

Authors

Contributions

POB conceptualized the study; NDM and DMT abstracted sequences from the Los Alamos HIV Sequence Database and analysed the data; POB and NDM wrote the first draft of the manuscript; DMT critically read the draft manuscript for intellectual input. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Pascal O. Bessong.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

DMT is an Associate Editor of AIDS Research and Therapy. NDM and POB declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Matume, N.D., Tebit, D.M. & Bessong, P.O. HIV-1 subtype C predicted co-receptor tropism in Africa: an individual sequence level meta-analysis. AIDS Res Ther 17, 5 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12981-020-0263-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12981-020-0263-x

Keywords