- Open Access
Subtype-specific differences in Gag-protease replication capacity of HIV-1 isolates from East and West Africa
Retrovirology volume 18, Article number: 11 (2021)
The HIV-1 epidemic in sub-Saharan Africa is heterogeneous with diverse unevenly distributed subtypes and regional differences in prevalence. Subtype-specific differences in disease progression rate and transmission efficiency have been reported, but the underlying biological mechanisms have not been fully characterized. Here, we tested the hypothesis that the subtypes prevalent in the East Africa, where adult prevalence rate is higher, have lower viral replication capacity (VRC) than their West African counterparts where adult prevalence rates are lower.
Gag-protease sequencing was performed on 213 and 160 antiretroviral-naïve chronically infected participants from West and East Africa respectively and bioinformatic tools were used to infer subtypes and recombination patterns. VRC of patient-derived gag-protease chimeric viruses from West (n = 178) and East (n = 114) Africa were determined using a green fluorescent protein reporter-based cell assay. Subtype and regional differences in VRC and amino acid variants impacting VRC were identified by statistical methods. CRF02_AG (65%, n = 139), other recombinants (14%, n = 30) and pure subtypes (21%, n = 44) were identified in West Africa. Subtypes A1 (64%, n = 103), D (22%, n = 35), or recombinants (14%, n = 22) were identified in East Africa. Viruses from West Africa had significantly higher VRC compared to those from East Africa (p < 0.0001), with subtype-specific differences found among strains within West and East Africa (p < 0.0001). Recombination patterns showed a preference for subtypes D, G or J rather than subtype A in the p6 region of gag, with evidence that subtype-specific differences in this region impact VRC. Furthermore, the Gag A83V polymorphism was associated with reduced VRC in CRF02_AG. HLA-A*23:01 (p = 0.0014) and HLA-C*07:01 (p = 0.002) were associated with lower VRC in subtype A infected individuals from East Africa.
Although prevalent viruses from West Africa displayed higher VRC than those from East Africa consistent with the hypothesis that lower VRC is associated with higher population prevalence, the predominant CRF02_AG strain in West Africa displayed higher VRC than other prevalent strains suggesting that VRC alone does not explain population prevalence. The study identified viral and host genetic determinants of virus replication capacity for HIV-1 CRF02_AG and subtype A respectively, which may have relevance for vaccine strategies.
Sub-Saharan Africa bears the brunt of the HIV-1 epidemic with approximately 70% of the estimated people living with the virus globally . The epidemic in sub-Saharan Africa is heterogeneous, with significant regional differences in prevalence and multiple HIV-1 subtypes that are unevenly distributed across the continent [2,3,4]. Specifically, West Africa is reported to have an adult HIV-1 prevalence rate of approximately 2%, although slight variations have been reported between countries . CRF02_AG is a circulating recombinant form (CRF) that constitutes the most common HIV-1 strain circulating in West Africa, however, almost all known pure subtypes have been identified within the region, as well as several CRFs and unique recombinant forms (URF) [2, 5,6,7]. In East Africa, HIV-1 prevalence rates are higher at above 5% [8, 9], with subtypes A1, C and D and their recombinants co-circulating. Multiple factors account for the heterogeneity of HIV prevalence and disease progression rates across sub-Saharan Africa, including sociocultural factors such as male circumcision practices, socioeconomic factors, coinfections, microbiota, host genetic and viral factors [10,11,12,13,14,15]. Considering that curtailing HIV transmission and providing optimal treatment remain formidable public health challenges worldwide, a better understanding of the viral factors that enhance transmission efficiency or the rate of disease progression following infection may be required for novel interventions such as vaccines .
It is well documented that HIV-1 subtypes display differences in the rate of disease progression. For example, in West Africa, individuals infected with non-A subtypes were more likely to develop AIDS compared to those with subtype A, whereas infection with the recombinant HIV-1 A3/CRF02_AG was associated with increased risk of AIDS and AIDS-related death compared to the sub-subtype A3 [17, 18]. In East Africa, HIV-1 subtype D has been associated with faster disease progression than subtype A, whereas in southern Africa, subtype C was associated with significantly slower rates of CD4 T-cell decline and higher frequencies of long-term non-progression compared to subtypes A or D [19,20,21,22]. Interestingly, studies have suggested that the more prevalent or efficiently transmitted subtypes are not necessarily the ones associated with faster disease progression, suggesting that the biological factors that determine transmission efficiency and in vivo virulence can be uncoupled and are distinct. For example, infections with HIV-1 subtype A were associated with higher transmission rates than subtype D [20, 23], which may have resulted in the increasing prevalence of subtype A over D within the study population over a period of time . Similarly, it has been proposed that a transmission-virulence evolutionary trade-off may explain the overall increasing relative prevalence of the more transmissible subtype A over the more virulent subtype D strain in Uganda; and subtype C is speculated to have become predominant in sub-Saharan Africa partly as a result of its lower virulence compared to subtypes A and D . Other studies have argued that socio-historical rather than evolutionary factors have been the main determinants of the rapid expansion of subtype C in sub-Saharan Africa .
Overall, the viral factors that may underlie variation in transmission efficiency and virulence are poorly understood. However, subtype-specific functional differences in various HIV-1 genetic loci are now well documented, suggesting that some of these may account for the distinct transmission patterns, epidemic spread, and rates of disease progression [27,28,29,30,31,32,33,34]. Gag-protease mediated VRC has been shown to differ between HIV-1 subtypes and correlates strongly with whole virus isolates, making it an important determinant of the rate of disease progression [35,36,37,38]. We and others have also shown that HIV-1 transmission selects for consensus-like viral variants and that similarity to consensus correlates with lower VRC [39, 40]. Overall, accumulating evidence suggests that the HIV-1 transmission genetic bottleneck favors consensus-like viruses that may have lower VRC, but that increased VRC in vitro is an important determinant of a faster rate of HIV disease progression. In this study, using chronically infected antiretroviral-naïve patient samples, we sought to further examine the diversity of circulating HIV-1 strains in West and East Africa, two regions of sub-Saharan African with distinct patterns of low versus moderately high prevalence of HIV and to investigate whether there are differences in gag-protease-mediated VRC that may explain the reported differences in prevalence and subtype-specific rates of disease progression. Furthermore, we also wanted to identify viral genetic and host genetic determinants of differences in VRC, considering that this information may have important implications for biomedical prevention and treatment strategies against HIV.
Materials and methods
This was a retrospective cross-sectional study, where ART-naive plasma samples collected between 2000 and 2010 from chronically infected individuals in West and East Africa were analyzed. In West Africa, study participants were from previously described cohorts from the Cameroon (n = 169) , Nigeria (n = 31)  and Senegal (n = 96) . Age, sex, and viral load data were available for most participants from Cameroon, and some of the participants had CD4 + T cell count data. Participants from Senegal were all female and only CD4 + T cell count data were available, while no demographic or clinical information was available for participants from Nigeria. In East Africa, study participants were from the well characterized IAVI Protocol-C cohort  from Kenya (n = 73), Rwanda (n = 61) and Uganda (n = 107) . However, samples from East Africa were preselected for subtypes A, AD and D based on previous pol gene sequencing . Age, sex, viral load, CD4 + T cell count and HLA class I data were available for participants from East Africa.
Viral RNA extraction
Viral RNA was extracted from chronically infected patient plasma samples using the QIAamp® Viral RNA Mini extraction kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. Samples with a viral load of less than 5,000 copies/mL were concentrated by spinning 500 µL of plasma at 14,000 × g for 2 h at 4 °C. Approximately 350 µL of supernatant was removed and discarded. The pellet was then re-suspended in the remaining volume of plasma and used for RNA extraction.
Amplification and viral sequencing
Extracted viral RNA was converted into cDNA and then double-stranded DNA by a one-step RT-PCR using the Superscript III One-Step RT-PCR, Platinum Taq High Fidelity kit (Invitrogen, San Diego, USA). Each reaction was made up of 20 µL of 2X reaction mix, 0.8 µL of each of the 10 µM forward Gag + 1 (5′-GAG GAG ATC TCT CGA CGC AGG AC-3′; HXB2 numbering 675–697) and reverse 3′RVP (5′-GGA GTG TTA TAT GGA TTT TCA GGC CCA ATT-3′; 2725–2696) primers, 0.8 µL of Superscript III RT/Platinum Taq High Fidelity Enzyme Mix and 4 µL of RNA. A second round of PCR was performed to amplify the 1,760 bp gag-protease, HXB2 positions 790–2550  using the TaKaRa Ex Taq HS enzyme kit (Shiga, Japan). Each reaction consisted of 5 µL of 10Ex Taq Buffer, 4 µL 2.5 mM dNTP, 0.8 µl each of 10 µM forward primer HXB2 695–794 (5′-GAC TCG GCT TGC TGA AGC GCG CAC GGC AAG AGG CGA GGG GCG GCG ACT GGT GAG TAC GCC AAA AAT TTT GAC TAG CGG AGG CTA GAA GGA GAG AGA TGG G-3′), and reverse primer HXB2 2,704–2,605 (5′-GGC CCA ATT TTT GAA ATT TTT CCT TCC TTT TCC ATT TCT GTA CAA ATT TCT ACT AAT GCT TTT ATT TTT TCT TCT GTC AAT GGC CAT TGT TTA ACT TTT G-3′) that are the NL4-3 sequences immediately flanking gag-protease, 0.25 µL of ExTaq enzyme and 2 µL of the RT-PCR product. Amplification of gag-protease was confirmed by agarose gel electrophoresis and amplicons were then sequenced as previously described . The resulting amplicon was diluted 1:15 in distilled water and population sequenced using the ABI PRISM Big Dye Terminator Ready reaction mix V3 (Applied Biosystems, Waltham, USA) and run on an ABI 3130xl Genetic Analyzer (Applied Biosystems, Waltham). Sequences were assembled and edited using the Sequencher 5.1 (Genecodes, USA) software program and were aligned to the HIV-1 subtype B reference strain HXB2 (GenBank accession no. K03455). Full-length gag-protease sequences were submitted under the accession numbers MW123927- MW124300.
Edited sequences were exported as a FASTA file and then aligned and translated using the Gene-Cutter tool from the Los Alamos National Laboratory (LANL) database (https://www.hiv.lanl.gov/content/sequence/HIV/mainpage.html). Sequences were uploaded to the REGA HIV-1 subtyping tool v3.46  and COMET subtyping tool  independently for subtype classification. Maximum-likelihood phylogenetic trees were constructed using PhyML  to determine evolutionary relationships between patient sequences and subtype reference strains obtained from the Los Alamos HIV sequence database https://www.hiv.lanl.gov/content/sequence/GENE_CUTTER/cutter.html. In instances where the subtyping methods were not in agreement, the phylogenetic tree was used for the final subtype assignment. Figtree v1.4.4 was used to infer the branching topology.
HIV-1 subtype recombination analysis
Inter-subtype recombination breakpoints were identified using the Jumping Profile Hidden Markov Model (jpHMM) online tool  and the Simplot  program to predict inter-subtype breakpoint coordinates. Default settings were used for jpHMM, while Simplot was optimized to a window size of 350, base pair step size of 30 and a consensus value of 50%. Each region between subtype breakpoints was extracted, trimmed, and realigned together with subtype reference sequences. Phylogenetic trees were then constructed to verify the subtype for each sequence fragment between predicted breakpoints.
Preparation of patient-derived Gag-protease NL4-3 virus stocks
The HIV-1 patient-derived gag-protease amplicons from the second round nested PCR were then co-transfected with linearized NL4-3Δgag-protease  by electroporation into CEM–GXR cells to generate chimeric viruses by homologous recombination as described previously . The electroporated CEM–GXR cells maintained in R10 medium were incubated at 37 °C and 5% CO2 for 12 days. The percentage of virus-infected cells was monitored by flow cytometry using a FACSCalibur (BD Biosciences, New Jersey, USA) since CEM-GXR cells express green fluorescent protein (GFP) when HIV-infected . Upon reaching a threshold infection rate of 25–30%, the cells were pelleted, and virus-containing supernatant was collected, aliquoted and stored at − 80 °C for subsequent cell culture experiments. To validate the viral stocks, a random subset was amplified and sequenced to confirm that the virus-derived sequences matched the original plasma HIV RNA sequences.
Titration of virus stocks
As previously described , 1 × 106 CEM–GXR cells were infected with 0.4 ml of virus stock and incubated for 48 h. The percentage of infected cells was then measured by flow cytometry, followed by analysis using Flowjo software (BD Biosciences, New Jersey, USA), to determine GFP expression (i.e., the percentage of infected cells) at 48 h post-infection. This percentage infected cells were used to calculate the viral volume required to achieve 0.3% (i.e., multiplicity of infection [MOI] of 0.003) infected cells after 48 h.
Replication capacity assay
The viral replication assay was performed in duplicate as previously described . Briefly, CEM-GXR cells were infected, using the same method as described in the titration assay, with the volume of virus stock calculated by the titration assay made up to a final volume of 0.4 ml using R10 medium. On days 2–6 following infection, the percentage infected cells were quantified by flow cytometry. A negative (R10 media alone) and positive (wild-type NL4-3 virus) control was included in each assay. The VRC of the chimeric viruses determined by gag-protease function was calculated by the mean slope of exponential growth from 3 to 6 days post infection using the semi-log method in Excel. These values were then normalized to those of the wild-type NL4-3 positive control.
VRCs of patient-derived chimeric viruses were compared based on HIV-1 subtype classification and geographical region using either the Student's t-test or Mann–Whitney U test if two groups were compared, or ANOVA with Tukey post-hoc tests where more than two groups were compared. ANOVA was used to test for differences in VRC across each group of the major HLA class I genes (A, B and C) for participants infected with subtype A1 from East Africa, while the Student's t-test was used to compare differences in VRC values between participants that expressed specific class I HLA alleles and those that did not express such alleles. HLA class I alleles expressed in a minimum of 5 individuals were included. These analyses were performed using GraphPad Prism v.8.4.3 (GraphPad Software, San Diego, California, USA). A multivariable linear regression model performed using Stata 15 (StataCorp LLC, Texas, USA), was used to assess the relationship between subtype and VRC after adjusting for potential confounders. A variable was included in the final model if its inclusion resulted in a 10% or greater change in the coefficients of the subtype category variable. Multiple imputation using chained equations, was used to assess the sensitivity of results to missing data. Codon-by-codon Mann–Whitney U tests with q values (available at https://bblab-hivresearchtools.ca/django/tools/codon_by_codon/) were used to identify specific amino acid variants associated with increased or decreased gag-protease-driven VRC.
Samples collected from HIV-1 chronically infected participants from West and East Africa were analyzed in this study. Of the 296 samples collected from West Africa, 213 (72%) successfully yielded gag-protease amplicons while 160 (66%) of the 241 samples from East Africa yielded PCR products. A total of 178 (84%) and 114 (71%) of amplified samples successfully yielded chimeric viruses and were assayed for VRC for West and East Africa, respectively. Table 1 shows a summary of the study participants demographic and clinical characteristics.
HIV-1 subtype diversity and distribution
A high level of subtype diversity was observed in West Africa, consistent with previous reports [6, 49]. CRF02_AG remained the most common circulating subtype in the West African isolates at 65%, followed by subtypes G and A3 at 7 and 5% respectively, while the remaining pure subtypes and CRFs were present at < 5% prevalence (Fig. 1a, b). Subtype assignment is shown in (Additional file 1: Tables S1 and S2). To facilitate comparison with West African samples and extend our previous analysis of gag-protease subtype in East Africa , we preselected subtypes A, D and AD recombinants based on previous pol subtyping data  since subtypes A, D and their recombinants have previously been demonstrated to dominate the epidemic in this region. The gag-protease subtyping data for the East African samples is shown in Fig. 2a, b. Though subtypes A1 and D were identified in both regions, both subtypes were significantly more prevalent in East Africa. Phylogenetic analysis showed clustering of sequences by both subtype and region (Additional file 1: Figure S1), indicating region-specific viral evolution and adaptation within subtypes.
Gag-protease inter-subtype recombination
Inter-subtype recombination analysis within the gag-protease region indicated that all HIV-1 CRFs and other recombinant forms in West and East Africa were recombinants of A1 (Fig. 3). In West Africa, the most prevalent inter-subtype recombinants were those of subtypes A1 and G, and less common inter-subtype recombinants were combinations of A1 and J, or A1, G and J mosaics, with CRF02_AG being the most common followed by CRF11_cpx. In the gag region, the remaining West African CRFs were close genetic relatives of subtype A1 (CRF01_AE, CRF09_cpx, CRF22_01A1, CRF36_cpx, CRF37_cpx, CRF45_cpx) or A1, G (CRF06_cpx and CRF43_02G) recombinants, with 96% (146/152) of these inter-subtype recombinants having a breakpoint within the p6 region in gag. In East Africa, we identified inter-subtype recombination between subtype A1 and C or A1 and D. All 5 A1C recombinants identified were from Rwanda and had a subtype A component towards the 3′ end of the sequences, while 80% (n = 12) of all A1D sequences from Uganda had a subtype D 3′ end (which included the entire p6 region). This recombination pattern is consistent with previous work from our group  indicating an evolutionary preference in A1D inter-subtype recombinant Gag sequences for subtype D towards the 3′ end of the sequence.
Gag-protease driven replication capacity of West and East African subtypes
Overall, prevalent viral isolates from West Africa had significantly higher VRC compared to those from East Africa (Fig. 4a). A subsequent comparison of the most prevalent circulating subtype in each region showed that CRF02_AG, which was most predominant in West Africa had a higher VRC than subtype A1 which was the most predominant subtype in the East African samples (Fig. 4b), indicating that the prevalent subtypes were driving the overall difference in VRC observed between both regions. Consistent with this idea, intra-subtype comparison of A1 and D between regions showed no difference in VRC according to region (Fig. 4c, d). In West Africa, pure subtype A1 and its close genetic relative A3/CRF22_01A1 had the lowest gag-protease driven VRC when compared to subtypes D, G, CRF02_AG, and CRF11_cpx (Fig. 4e). In East Africa, pure subtype A1 also had the lowest gag-protease driven VRC when compared to other subtypes, and a hierarchy of A/A1C < A1D < D was observed (Fig. 4f), which is consistent with previous reports by our group .
Multivariate regression analysis (West and East Africa)
Men have been shown to have higher viral loads than women [23, 50]. Viral load as well as CD4 + T cell count are significantly associated with gag-protease driven VRC [35, 37] Therefore, a multivariable linear regression was used to assess the relationship between subtype and VRC after adjusting for potential confounders, including participant’s sex, age, viral load and CD4 + T cell counts (Table 1).
In West Africa, the most frequently observed subtype was CRF02_AG which was observed in 139 (65%) of participants. In the univariate analysis, subtypes A1, A1G recombinants, A3 and CRF22_01A1 were all associated with significantly lower VRC when compared to CRF 02_AG. Less than 25% of participants had complete information on all the following parameters: CD4 + count, viral load, and VRC, rendering the construction of a multivariable or imputation model challenging. In view of country-specific missing data patterns, we constructed separate multivariable models for Cameroon and Senegal (but not Nigeria due to the small sample size) to gain insight into the effect of subtype on VRC while adjusting for confounders (Table 2). The country-specific multivariate models collectively supported the significant findings in the univariate analysis: A1 (0.49, p < 0.0001), A1G recombinants (0.35, p = 0.01), A3 (0.2, p = 0.03) and CRF22_01A1 (0.4, p < 0.0001) had significantly lower mean VRC compared to CRF02_AG.
In East Africa, a total of 114 participants had complete information on VRC and were included in the regression models. The A1D and D subtypes were associated with significantly higher VRC values, after adjusting for age, sex, country, CD4 + count and viral load. Using multiple imputation with multivariate regression the dataset was imputed 50 times and estimates were combined on the imputed dataset, yielding revised regression estimates. After multiple imputation A1D had a 0.24 significantly higher mean VRC compared to A1 (p < 0.0001), while subtype D had a 0.29 significantly higher mean VRC compared to A1 (p < 0.0001) (Table 3) after adjusting for age, sex, country, CD4 + T cell count and viral load.
Association of HLA class I alleles with gag-protease-driven replication capacity
Genome wide association studies have demonstrated that the human leukocyte antigen class I (HLA-I) alleles are the most significant genetic determinant of clinical outcome in HIV-1 infection . To investigate the impact of HLA-I allele expression on the Gag-protease driven VRC of subtype A1 patient-derived isolates from East Africa, VRC data were grouped according to HLA-I alleles expressed by the hosts for each of the class I loci (Fig. 5). VRC did not differ significantly across alleles in any of the class I loci, however analysis of individual HLA-I alleles showed that HLA-A*23:01 and HLA-C*07:01 were associated with significantly lower VRC (Student's t-test, p = 0.0014 and p = 0.002) respectively, while HLA-B*07:02 was associated with significantly higher VRC (Student's t-test, p = 0.004).
Gag p6 subtype predicts VRC
Recombination patterns in West and East African sequences showed a preference for subtypes D, G or J rather than the predominant subtype A in the p6 region of Gag, highlighting an important evolutionary trait that may impact VRC. In support of this, a multiple sequence alignment of URFs from East Africa was performed, and the p6 region of the alignment was extracted and confirmed by phylogenetic analysis. Corresponding VRC data were categorized based on subtype of the p6 region. The results showed that in East Africa (Fig. 6a), p6 regions with subtype A1 had significantly lower (p = 0.0021) VRC than those with subtype D. A corresponding analysis could not be done for West African sequences as we lacked sufficient sequences with a subtype A1 p6 region for the West African CRFs.
To investigate subtype-specific amino acid sequences in the p6 region of Gag, consensus sequences of individual subtypes irrespective of region were generated A1 (n = 104), A1D (n = 16), A3 (n = 10), CRF02_AG (n = 138), CRF11_cpx (n = 5), CRF22_01A1 (n = 8), D (n = 43) and G (n = 15) and aligned with the HXB2 subtype B sequence as a reference. While amino acid sequence in the PTAP domain remained conserved irrespective of subtype, there were marked subtype-specific variations in the LYPLASL domain in the Alix budding motif. The subtypes with above median VRC are conserved for the LY residues except for subtype G which has LP as the consensus at this position in our sequences, while the subtypes with below median VRC have other amino acids at the same position (Fig. 6b), suggesting, along with previous reports [38, 52], that variation in the Alix motif could impact VRC.
An exploratory codon-by-codon analysis was performed to identify amino acid polymorphic variants associated with differences in VRC for the predominant subtypes in West (CRF02_AG) and East (A1) Africa using methods previously described . The analysis identified several amino acid polymorphisms that associated with altered VRC at p < 0.05, although after correcting for multiple codon comparisons none of the associations met the significant threshold of q < 0.2 (Tables 4 and 5).
In CRF02_AG sequences, A83V, a previously inferred escape variant  in known CTL epitopes in the p17 region of gag was the only non-consensus amino acid variant associated with a lower VRC (Table 4). This result was consistent with previous reports linking A83V to reduced viral fitness in subtypes B and CRF01_AE [54, 55]. In subtype A sequences, the only non-consensus polymorphism associated with decreased VRC was P497L, while non-consensus amino acid variants K69Q, S125N, L147I, A158V and R380K, were associated with significantly higher VRC (Table 5). However, the consensus amino acids at positions 12 and 28 were associated with increased VRC, which suggests that mutations at these codons, within known CTL epitopes, confer a fitness cost. Consistent with this, polymorphisms at codon 28 have been reported to alter VRC in other subtypes [40, 56]. Further work is required to validate these findings since the relatively small number of sequences available for analysis resulted in limited statistical power.
A hallmark of the HIV-1 epidemic in sub-Saharan Africa is the region-specific differences in prevalence such that southern African countries tend to have the highest prevalence rates followed by East African countries, with West and Central African countries having the lowest prevalence rates [57,58,59]. It is also noteworthy that HIV-1 subtypes are unevenly distributed across the continent, with subtype C predominant in southern Africa, subtypes A, D and C common in East Africa and almost all subtypes present in West and Central Africa [59, 60]. The possible contribution of viral genetic and functional characteristics to the uneven distribution of subtypes is unresolved but differences in transmission efficiency and rate of disease progression have been confirmed in epidemiological and clinical studies, suggesting that viral factors may partially explain the heterogeneity in prevalence and uneven spread of HIV-1 subtypes within the continent [18, 19, 26, 38, 61,62,63]. Here, using samples from HIV-1 chronically infected individuals from 3 West African countries, we confirm findings from previous studies that the epidemic in this region is genetically heterogenous with CRF02_AG (65%) the predominant circulating strain but with other subtypes and recombinant forms each contributing a significant proportion of the epidemic in the region [7, 9, 59]. Our study was not designed to explore the genetic diversity of the HIV-1 epidemic in East Africa since the samples analyzed for this study were preselected to be either subtypes A, D or their recombinants based on previous pol sequencing , however, numerous other studies have demonstrated a far less diverse epidemic in East Africa dominated by subtypes A, D, C and their recombinant forms [19, 20, 38, 64, 65].
Recombination analyses of the West African isolates identified several CRFs, consistent with the multiple pure subtypes that have been identified within the region . In West Africa, most recombinants comprised of A1 and G genetic fragments and breakpoint analysis suggested that recombination patterns were not random, with a strong preference for subtype G or J over A1 at the 3′-end of most recombinant sequences, with a recombination hotspot around the Gag p6 region. A similar trend was noted for the East African isolates where there was also a preference for subtype D over A for the Gag p6 region in recombinant sequences. Overall, the data suggested that subtype A1 in both West and East Africa has a high propensity for inter-subtype recombination with non-A subtypes at the 3′ end, preferred particularly from around the Gag p6 region. We hypothesize that this could be a mechanism to facilitate the adaptation of HIV-1 subtype A into a fitter virus. Indeed, it has been reported that throughout sub-Saharan Africa, pure forms of subtype A are in decline, with a concomitant increase in its recombinants .
The notion that recombination may be a non-random mechanism to enhance viral fitness is backed up by our VRC data. In West and East Africa, gag-protease-driven VRC differed by HIV-1 subtype with a clear hierarchy where A1/A3 had generally lower VRC compared to pure subtypes D and G or A1 recombinants (although this was not the case for A1/C recombinants and CRF22_01A1). Our data is thus consistent with previous studies of Gag-protease functional differences among HIV-1 subtypes and recombinants that may also extend to other HIV-1 proteins [33, 35, 53, 67, 68]. Enhanced gag-protease-driven VRC has been associated with faster rate of CD4 + T cell decline and disease progression , and our data is consistent with previously reported subtype-specific differences in disease progression [17,18,19, 22].
Interestingly, in East Africa, it is subtype A1, with a lower replication capacity than subtype D that has been reported to be more transmissible; despite the latter being associated with faster disease progression [20, 21, 38, 64, 65]. This observation, together with reports that the HIV-1 transmission bottleneck may favor consensus-like viruses, which have been associated with lower replication capacity [35, 37, 39, 40], led us to our hypothesis that West African viral isolates, where HIV prevalence is low may have higher replication capacity compared to the viruses from East Africa where prevalence is higher. Although we found that West African viruses have higher overall viral replication capacity compared to East African viruses, these differences are explained by the subtypes prevalent in these regions and are not region-specific per se. Moreover, in contrast to the situation in East Africa where a virus strain (subtype A) with lower replication capacity predominates, CRF02_AG is the predominant strain in West Africa and yet it has a higher viral replication capacity than other strains prevalent in the region. Overall, our data suggest that viral replication capacity alone does not explain the epidemiological success of a virus strain. Our data highlights the need for further studies to better understand the viral determinants that may underlie transmissibility and virulence.
Recombination patterns suggested that the 3′ region of the Gag protein is a recombination hotspot particularly for A1 recombinants. Specifically, the non-A subtype was preferred for the p6 region of Gag. The p6 contains two well studied segments, the PTAP and LYPLASL domains which interact with the host cell TSG101 and Alix factors respectively to facilitate efficient virion maturation and budding from host cells [52, 70]. Whereas the PTAP motif showed no amino acid variation across subtypes, the LYPLASL motif showed variation across subtypes which corresponded to either low or high VRC and a comparative analysis of East African recombinants with A1 and D components within the p6 region showed a lower VRC for those with subtype A1 (p < 0.05). We and others have shown previously that variation in the Alix budding motif in Gag affects VRC [38, 71]. Further work to interrogate the effect of these subtype-specific amino acid variations within the LYPLASL motifs on VRC is warranted.
Finally, our study identified HLA-A class I alleles that have differential impact on HIV-1 VRC for HIV-1 subtype A, the predominant subtype in East Africa which suggest the need for further studies to identify regions of viral vulnerability for HIV-1 vaccine design strategies. Furthermore, viral replication capacities did not differ across HLA class I loci for HIV-1 subtype A in East Africa, even for HLA-B alleles, contrary to earlier observations, particularly in subtype B and C infections, of strong selective pressure by certain HLA-B alleles that may impact replication capacity [13, 35, 56, 72]. However, individually, HLA-A*23:01 and HLA-C*07:01 were associated with lower replication capacity, suggesting an underappreciated role in immune selection pressure on HIV-1 subtype A by certain HLA-A and HLA-C alleles. This finding appears consistent with previous findings of subtype-specific differences in HLA-driven viral evolution that may have consequences for natural and vaccine-mediated immunity . Overall, the data highlight the need for further studies to identify mechanisms of immune control or regions of viral vulnerability by HLA alleles common in regions of the world with non-B and C subtypes as this knowledge may be useful for universal virus attenuation-based vaccine strategies.
Limitations of the current study include missing clinical data such as CD4 + T cell counts and HLA-I genotypes, for the West African samples, thus limiting the extent of statistical analysis. East African samples were preselected for subtypes A, D and their recombinants based on prior pol genotyping, therefore not allowing comprehensive inter-subtype comparison of the strains that constitute the East Africa epidemic. It should also be noted that all replication assays were performed in the subtype B pNL4-3∆gag-protease backbone, however, we have previously demonstrated that VRC data generated from this assay is generally reflective of whole virus isolates in HIV-1 subtype C [35, 38]. A higher sample size would also have benefited the robustness of the codon-by-codon analysis to identify amino acids associated with altered VRC for different subtypes and recombinants.
Overall, in this study we show that the most prevalent HIV-1 viruses from West Africa display higher replication capacities than those from East Africa, consistent with the hypothesis that higher prevalence is associated with lower replication capacity of circulating strains. Subtype-specific differences in replication capacity in agreement with previous studies were noted, and consistent with reported differences in the rate of clinical disease progression. Our study identified inter-subtype recombination patterns as a driver of virus replication capacity differences and identified HLA class I alleles and specific amino acids that may alter virus replication, information that may be relevant for HIV-1 vaccine design strategies.
Availability of data and materials
The datasets used and analysed during this study are available from the corresponding author.
Viral replication capacity
Circulating recombinant Form
Unique recombinant Form
Jumping profile Hidden Markov model
Green fluorescent protein
Human leucocyte antigen
Complementary deoxyribonucleic acid
Reverse transcriptase polymerase chain reaction
Los Alamos National Laboratory
Multiplicity of infection
Human immunodeficiency virus type 1
UNAIDS. UNAIDS DATA 2019. 2019.
Lihana RW, Ssemwanga D, Abimiku A, Ndembi N. Update on HIV-1 diversity in Africa: a decade in review. AIDS Rev. 2012;14(2):83–100.
Morison L, Buve A, Zekeng L, Heyndrickx L, Anagonou S, Musonda R, et al. HIV-1 subtypes and the HIV epidemics in four cities in sub-Saharan Africa. AIDS. 2001;15(Suppl 4):S109–16.
Santos AF, Soares MA. HIV Genetic Diversity and Drug Resistance. Viruses. 2010;2(2):503–31.
Djomand G, Quaye S, Sullivan PS. HIV epidemic among key populations in west Africa. Curr Opin HIV AIDS. 2014;9(5):506–13.
Tongo M, Martin DP, Zembe L, Mpoudi-Ngole E, Williamson C, Burgers WA. Characterization of HIV-1 gag and nef in Cameroon: further evidence of extreme diversity at the origin of the HIV-1 group M epidemic. Virol J. 2013;10:29.
Bbosa N, Kaleebu P, Ssemwanga D. HIV subtype diversity worldwide. Curr Opin HIV AIDS. 2019;14(3):153–60.
Kharsany AB, Karim QA. HIV infection and AIDS in Sub-Saharan Africa: current status challenges and opportunities. Open AIDS J. 2016;10:34–48.
Fettig J, Swaminathan M, Murrill CS, Kaplan JE. Global epidemiology of HIV. Infect Dis Clin North Am. 2014;28(3):323–37.
Serwadda D, Wawer MJ, Musgrave SD, Sewankambo NK, Kaplan JE, Gray RH. HIV risk factors in three geopgraphic strata of rural Rakai District. Uganda AIDS. 1992;6:983–90.
Auvert B, Buonamico G, Lagarde E, Williams B. Sexual behavior, heterosexual transmission, and the spread of HIV in sub-Saharan Africa: a simulation study. Comput Biomed Res. 2000;33(1):84–96.
Weiss HA, Quigley MA. Male circumcision and risk of HIV infection in sub-Saharan Africa: a systematic review and meta-analysis. AIDS. 2000. https://doi.org/10.1097/00002030-200010200-00018.
Kiepiela P, Leslie AJ, Honeyborne I, Ramduth D, Thobakgale C, Chetty S, et al. Dominant influence of HLA-B in mediating the potential co-evolution of HIV and HLA. Nature. 2004;432(7018):769–75.
Fellay J, Shianna KV, Ge D, Colombo S, Ledergerber B, Weale M, et al. A whole-genome association study of major determinants for host control of HIV-1. Science. 2007;317(5840):944–7.
Gosmann C, Anahtar MN, Handley SA, Farcasanu M, Abu-Ali G, Bowman BA, et al. Lactobacillus-Deficient Cervicovaginal Bacterial Communities Are Associated with Increased HIV Acquisition in Young South African Women. Immunity. 2017;46(1):29–37.
Barouch DH. Challenges in the development of an HIV-1 vaccine. Nature. 2008;455(7213):613–9.
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.
Kanki PJ, Hamel DJ, Sankale JL, Hsieh C, Thior I, Barin F, et al. Human immunodeficiency virus type 1 subtypes differ in disease progression. J Infect Dis. 1999;179(1):68–73.
Kaleebu P, Ross A, Morgan D, Yirrell D, Oram J, Rutebemberwa A, et al. Relationship between HIV-1 Env subtypes A and D and disease progression in a rural Ugandan cohort. AIDS. 2001;15(3):293–9.
Kiwanuka N, Laeyendecker O, Robb M, Kigozi G, Arroyo M, McCutchan F, et al. Effect of human immunodeficiency virus Type 1 (HIV-1) subtype on disease progression in persons from Rakai, Uganda, with incident HIV-1 infection. J Infect Dis. 2008;197(5):707–13.
Ssemwanga D, Nsubuga RN, Mayanja BN, Lyagoba F, Magambo B, Yirrell D, et al. Effect of HIV-1 subtypes on disease progression in rural Uganda: a prospective clinical cohort study. PLoS ONE. 2013;8(8):e71768.
Venner CM, Nankya I, Kyeyune F, Demers K, Kwok C, Chen PL, et al. Infecting HIV-1 Subtype Predicts Disease Progression in Women of Sub-Saharan Africa. EBioMedicine. 2016;13:305–14.
Amornkul PN, Karita E, Kamali A, Rida WN, Sanders EJ, Lakhi S, et al. Disease progression by infecting HIV-1 subtype in a seroconverter cohort in sub-Saharan Africa. AIDS. 2013;27(17):2775–86.
Conroy SA, Laeyendecker O, Redd AD, Collinson-Streng A, Kong X, Makumbi F, et al. Changes in the distribution of HIV type 1 subtypes D and A in Rakai District, Uganda between 1994 and 2002. AIDS Res Hum Retroviruses. 2010;26(10):1087–91.
Blanquart F, Grabowski MK, Herbeck J, Nalugoda F, Serwadda D, Eller MA, et al. A transmission-virulence evolutionary trade-off explains attenuation of HIV-1 in Uganda. Elife. 2016;5:e20492.
Faria NR, Vidal N, Lourenco J, Raghwani J, Sigaloff KCE, Tatem AJ, 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.
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;83(11):5592–605.
Velazquez-Campoy A, Todd MJ, Vega S, Freire E. Catalytic efficiency and vitality of HIV-1 proteases from African viral subtypes. PNAS. 2001. https://doi.org/10.1073/pnas.111152698.
Bachu M, Yalla S, Asokan M, Verma A, Neogi U, Sharma S, et al. Multiple NF-kappaB sites in HIV-1 subtype C long terminal repeat confer superior magnitude of transcription and thereby the enhanced viral predominance. J Biol Chem. 2012;287(53):44714–35.
Ball SC, Abraha A, Collins KR, Marozsan AJ, Baird H, Quinones-Mateu ME, et al. Comparing the ex vivo fitness of CCR5-tropic human immunodeficiency virus type 1 isolates of subtypes B and C. J Virol. 2003;77(2):1021–38.
Mann JK, Byakwaga H, Kuang XT, Le AQ, Brumme CJ, Mwimanzi P, Omarjee S, Martin E, Lee GQ, Baraki B, Danroth R. Ability of HIV-1 Nef to downregulate CD4 and HLA class I differs among viral subtypes. Retrovirology. 2013;10(1):1–6.
Montano MA, Nixon CP, Essex M. Dysregulation through the NF-kB enhancer and TATA box of the human immunodeficiency virus type 1 subtype E promoter. J Virol. 1998;72(10):8446–52.
Iordanskiy S, Waltke M, Feng Y, Wood C. Subtype-associated differences in HIV-1 reverse transcription affect the viral replication. Retrovirology. 2010;7:85.
Umviligihozo G, Cobarrubias KD, Chandrarathna S, Jin SW, Reddy N, Byakwaga H, et al. Differential Vpu-mediated CD4 and tetherin downregulation functions among major HIV-1 Group M subtypes. J Virol. 2020. https://doi.org/10.1128/JVI.00293-20.
Wright JK, Brumme ZL, Carlson JM, Heckerman D, Kadie CM, Brumme CJ, et al. Gag-protease-mediated replication capacity in HIV-1 subtype C chronic infection: associations with HLA type and clinical parameters. J Virol. 2010;84(20):10820–31.
Claiborne DT, Prince JL, Scully E, Macharia G, Micci L, Lawson B, et al. Replicative fitness of transmitted HIV-1 drives acute immune activation, proviral load in memory CD4+ T cells, and disease progression. Proc Natl Acad Sci U S A. 2015;112(12):E1480–9.
Prince JL, Claiborne DT, Carlson JM, Schaefer M, Yu T, Lahki S, et al. Role of transmitted Gag CTL polymorphisms in defining replicative capacity and early HIV-1 pathogenesis. PLoS Pathog. 2012;8(11):e1003041.
Kiguoya MW, Mann JK, Chopera D, Gounder K, Lee GQ, Hunt PW, et al. Subtype-specific differences in Gag-protease-driven replication capacity are consistent with intersubtype differences in HIV-1 disease progression. J Virol. 2017. https://doi.org/10.1128/JVI.00253-17.
Deymier MJ, Ende Z, Fenton-May AE, Dilernia DA, Kilembe W, Allen SA, et al. Heterosexual transmission of subtype C HIV-1 selects consensus-like variants without increased replicative capacity or interferon-alpha resistance. PLoS Pathog. 2015;11(9):e1005154.
Naidoo VL, Mann JK, Noble C, Adland E, Carlson JM, Thomas J, et al. Mother-to-child HIV transmission bottleneck selects for consensus virus with lower gag-protease-driven replication capacity. J Virol. 2017. https://doi.org/10.1128/JVI.00518-17.
Ajoge HO, Gordon ML, de Oliveira T, Green TN, Ibrahim S, Shittu OS, et al. Genetic characteristics, coreceptor usage potential and evolution of Nigerian HIV-1 subtype G and CRF02_AG isolates. PLoS ONE. 2011;6(3):e17865.
Hamel DJ, Sankale JL, Eisen G, Meloni ST, Mullins C, Gueye-Ndiaye A, et al. Twenty years of prospective molecular epidemiology in Senegal: changes in HIV diversity. AIDS Res Hum Retroviruses. 2007;23(10):1189–96.
Price MA, Rida W, Kilembe W, Karita E, Inambao M, Ruzagira E, et al. Control of the HIV-1 load varies by viral subtype in a large cohort of african adults with incident HIV-1 infection. J Infect Dis. 2019;220(3):432–41.
Tözsér J, Bláha I, Copeland TD, Wondrak EM, Oroszlan S. Comparison of the HIV-1 and HIV-2 proteinases using oligopeptide substrates representing cleavage sites in Gag and Gag-Pol polyproteins. FEBS Lett. 1991;281(1–2):77–80. https://doi.org/10.1016/0014-5793(91)80362-7.
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 Res. 2009;37:W647–51.
Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, Ingersoll R, Sheppard HW, Ray SC. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999;73(1):9.
Miura T, Brockman MA, Brumme ZL, Brumme CJ, Pereyra F, Trocha A, et al. HLA-associated alterations in replication capacity of chimeric NL4-3 viruses carrying gag-protease from elite controllers of human immunodeficiency virus type 1. J Virol. 2009;83(1):140–9.
Brockman MA, Tanzi GO, Walker BD, Allen TM. Use of a novel GFP reporter cell line to examine replication capacity of CXCR4- and CCR5-tropic HIV-1 by flow cytometry. J Virol Methods. 2006;131(2):134–42.
Edson Delatorre DM, Bello G. Spatiotemporal dynamics of the HIV-1 subtype G epidemic in West and Central Africavv. PLoS ONE. 2014;9(6):e98908.
MonicaGandhi PB, Miotti P, Quinn TC, Veronese F, Greenblatt RM. Does patient sex affect human immunodeficiency virus levels? Clin Infect Dis. 2001;35:313–22.
Pereyra F, Jia X, McLaren PJ, Telenti A, de Bakker PI, Walker BD, et al. The major genetic determinants of HIV-1 control affect HLA class I peptide presentation. Science. 2010;330(6010):1551–7.
Barbara Muller TP, Hans-Georg K. The late-domain-containing protein p6 is the predominant phosphoprotein of human immunodeficiency virus type 1 particles. J Virol. 2002;76(3):10.
Ojwach DBA, Reddy T, Novitsky V, Brumme ZL, Brockman MA, Ndung’u T, Mann JK. Pol-driven replicative capacity impacts disease progression in hiv-1 subtype C infection. J Virol. 2018;92(19):eoo811.
Boutwell CL, Carlson JM, Lin TH, Seese A, Power KA, Peng J, et al. Frequent and variable cytotoxic-T-lymphocyte escape-associated fitness costs in the human immunodeficiency virus type 1 subtype B Gag proteins. J Virol. 2013;87(7):3952–65.
Takahashi N, Matsuoka S, Thi Minh TT, Ba HP, Naruse TK, Kimura A, et al. Human leukocyte antigen-associated gag and nef polymorphisms in HIV-1 subtype A/E-infected individuals in Vietnam. Microbes Infect. 2019;21(2):113–8.
Brockman MA, Brumme ZL, Brumme CJ, Miura T, Sela J, Rosato PC, et al. Early selection in Gag by protective HLA alleles contributes to reduced HIV-1 replication capacity that may be largely compensated for in chronic infection. J Virol. 2010;84(22):11937–49.
Shao Y, Williamson C. The HIV-1 epidemic: low- to middle-income countries. Cold Spring Harb Perspect Med. 2012;2(3):a007187.
Abdool KQ. The global HIV epidemic: current status and challenges. Curr HIV/AIDS Rep. 2013;10(2):111–2.
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;19(2):143–55.
JoriHemelaar EG, Ghys PD, Osmanov S. Global and regional distribution of HIV-1 genetic subtypes and recombinants in 2004. AIDS. 2006;20:11.
Blackard JT, Renjifo B, Fawzi W, Hertzmark E, Msamanga G, Mwakagile D, et al. HIV-1 LTR subtype and perinatal transmission. Virology. 2001;287(2):261–5.
Odaibo GN, Heyndrickx L, Vereecken K, Houwer K, Jassens W. Mother-to-child transmission of different hiv-1 subtypes among arv naive iinfected pregnant woment in Nigeria. Rev Inst Med Trop S Paulo. 2006;48(2):77–80.
Kiwanuka N, Laeyendecker O, Quinn TC, Wawer MJ, Shepherd J, Robb M, et al. HIV-1 subtypes and differences in heterosexual HIV transmission among HIV-discordant couples in Rakai. Uganda AIDS. 2009;23(18):2479–84.
Vasan A, Renjifo B, Hertzmark E, Chaplin B, Msamanga G, Essex M, et al. Different rates of disease progression of HIV type 1 infection in Tanzania based on infecting subtype. Clin Infect Dis . 2006;42(6):843–52.
Baeten JM, Chohan B, Lavreys L, Chohan V, McClelland RS, Certain L, et al. HIV-1 subtype D infection is associated with faster disease progression than subtype A in spite of similar plasma HIV-1 loads. J Infect Dis. 2007;195(8):1177–80.
Sharp PM, Hahn BH. Origins of HIV and the AIDS pandemic. Cold Spring Harb Perspect Med. 2011;1(1):a006841.
Jaclyn K, et al. Ability of HIV-1 Nef to downregulate CD4 and HLA class I differs among viral subtypes. Retrovirology. 2013;10(100):1–16.
Jin SW, Mwimanzi FM, Mann JK, Bwana MB, Lee GQ, Brumme CJ, et al. Variation in HIV-1 Nef function within and among viral subtypes reveals genetically separable antagonism of SERINC3 and SERINC5. PLoS Pathog. 2020;16(9):e1008813.
Wright JK, Novitsky V, Brockman MA, Brumme ZL, Brumme CJ, Carlson JM, et al. Influence of Gag-protease-mediated replication capacity on disease progression in individuals recently infected with HIV-1 subtype C. J Virol. 2011;85(8):3996–4006.
Solbak SM, Reksten TR, Hahn F, Wray V, Henklein P, Henklein P, et al. HIV-1 p6 - a structured to flexible multifunctional membrane-interacting protein. Biochim Biophys Acta. 2013;1828(2):816–23.
Aralaguppe SG, Winner D, Singh K, Sarafianos SG, Quinones-Mateu ME, Sonnerborg A, et al. Increased replication capacity following evolution of PYxE insertion in Gag-p6 is associated with enhanced virulence in HIV-1 subtype C from East Africa. J Med Virol. 2017;89(1):106–11.
Kawashima Y, Pfafferott K, Frater J, Matthews P, Payne R, Addo M, et al. Adaptation of HIV-1 to human leukocyte antigen class I. Nature. 2009;458(7238):641–5.
McKinnon LR, Capina R, Peters H, Mendoza M, Kimani J, Wachihi C, et al. Clade-specific evolution mediated by HLA-B*57/5801 in human immunodeficiency virus type 1 clade A1 p24. J Virol. 2009;83(23):12636–42.
This work was supported by the Poliomyelitis Research Foundation (PRF), the South African Department of Science and Innovation through the National Research Foundation (South African Research Chairs Initiative), the Victor Daitz Foundation and the International AIDS Initiative (UKZNRSA1001). Additional funding was provided by the Sub-Saharan African Network for TB/HIV Research Excellence (SANTHE), a DELTAS Africa Initiative [grant # DEL-15-006]. The DELTAS Africa Initiative is an independent funding scheme of the African Academy of Sciences (AAS)’s Alliance for Accelerating Excellence in Science in Africa (AESA) and supported by the New Partnership for Africa’s Development Planning and Coordinating Agency (NEPAD Agency) with funding from the Wellcome Trust [grant # 107752/Z/15/Z] and the UK government. The views expressed in this publication are those of the author(s) and not necessarily those of AAS, NEPAD Agency, Wellcome Trust, or the UK government. We are also grateful to staff and volunteers from IAVI’s protocol C sites in Africa. This project was also made possible in part by the generous support of the American people through the United States Agency for International Development (USAID). The contents are the responsibility of the study authors and do not necessarily reflect the views of USAID or the United States Government.
Ethics approval and consent to participate
Ethical approval for this study was obtained from the Biomedical Research Ethics Committee of the University of KwaZulu-Natal (BE 313/16). Written informed consent was obtained from all study participants at the time of recruitment and patient identity remained confidential.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Farinre, O., Gounder, K., Reddy, T. et al. Subtype-specific differences in Gag-protease replication capacity of HIV-1 isolates from East and West Africa. Retrovirology 18, 11 (2021). https://doi.org/10.1186/s12977-021-00554-4
- HIV-1 subtype
- HIV-1 recombinants
- Sub-Saharan Africa
- HIV-1 replication capacity
- HLA class I alleles