- Open Access
Neutralization diversity of HIV-1 Indian subtype C envelopes obtained from cross sectional and followed up individuals against broadly neutralizing monoclonal antibodies having distinct gp120 specificities
Retrovirology volume 18, Article number: 12 (2021)
The potential use of the broadly neutralizing monoclonal antibodies (bnAbs) towards prophylaxis and treatment to HIV-1 is currently being explored. While a number of promising bnAbs have been discovered and a few of them have progressed towards clinical development, their extent of neutralization coverage with respect to global HIV-1 variants given the existence of genetically distinct subtypes and recombinants circulating globally is not clearly known. In the present study, we examined the variation in the neutralization susceptibility of pseudoviruses expressing 71 full length primary HIV-1 subtype C envs obtained from limited cross-sectional individuals over different time points against four bnAbs that target gp120 with distinct specificities: VRC01, CAP256-VRC26.25, PGDM1400 and PGT121.
We found significant variations in the susceptibility of Indian clade C to these four bnAbs. These variations were found to be distinct to that observed in African subtype C based on the existing datasets and concordant with their sequence diversity. Trend analysis indicated an increasing neutralization resistance observed over time with CAP25-VRC26.25, PGDM1400 and PGT121 when tested on pseudoviruses expressing envs obtained from 1999 to 2016. However, inconsistent trend in neutralization susceptibility was observed, when pseudoviruses expressing envs obtained from three followed up individuals were examined. Finally, through predictive analysis of the 98 Indian subtype C including those assessed in the present study by employing additive model implemented in CombiNAber (http://www.hiv.lanl.gov), we observed two possibilities where combinations of three bnAbs (VRC01/CAP56-VRC26.25/PGT121 and PGDM1400/CAP256-VRC26.25/PGT121) could achieve near 100% neutralization coverage.
Our findings not only indicate disparate intra-clade C genetic vis-à-vis neutralization diversities but also warrant the need for more comprehensive study using additional isolates towards comparing inter and intra-clade neutralization diversities which will be necessary for selecting the bnAb combinations suitable for optimal coverage of the region-specific HIV-1 circulating subtypes. Expanding these efforts is imperative for designing efficacious bnAb based intervention strategies for India as well as subtype C in general.
The development of a preventive vaccine for HIV-1 with a potential to tackle the enormous genetic diversity remains a challenge. This genetic variation is primarily accounted for by the HIV-1 envelope gene (env) which is responsible for greater than 35% difference in amino acid sequences between different subtypes . A number of broadly neutralizing monoclonal antibodies (bnAbs) capable of cross neutralizing diverse genetic subtypes have been discovered from elite neutralizers since 2009, which have fuelled in recent times their potential use in prevention and therapy over and above the existing antiretroviral (ARV) drugs [2,3,4,5,6,7,8,9]. Few of these bnAbs have also advanced into clinical development based on their promising antiviral activity observed in animals and human clinical trials [10,11,12,13,14,15,16]. In general, breadth and potency of bnAbs are defined by their extent in neutralizing a panel of pseudoviruses expressing envs derived from various subtypes and recombinants representing different geographic regions. Hence, the suitability of bnAbs that exhibited the highest breadth to overcome region-specific diversity is not clearly known. Previous studies have shown that although bnAbs isolated from individuals infected with one particular subtype are generally effective at neutralizing viruses belonging to other subtypes, antibody potency is often found to be correlated with matched subtypes [17,18,19,20,21]. Moreover, viral and host diversities have been found to have an impact even within a matched subtype, as demonstrated by the fact that the subtype-matched neutralization advantage was more apparent in regions with distinct viral diversities [19, 22, 23].
HIV-1 subtype C accounts for approximately half of the global infections , which predominates in India and South Africa. Recently very few studies have attempted to understand a few selected bnAbs for their extent to neutralize HIV-1 subtype C of African origin either singly and/or in combinations [25,26,27], however, those studies did not include pseudoviruses expressing envs representing Indian subtype C. It is possible that an intra-clade C specific neutralization patterns may exist. For example, Rademeyer et al.  have shown that subtype C viruses of African origin have become more resistant to VRC01, PG9 and 4E10 compared to CAP256-VRC26.25 and PGT128. In the present study, we examined the variation in neutralization sensitivity of pseudoviruses expressing 71 complete HIV-1 subtype C primary envs collected over a period of time (1999–2014) against four bNAbs (VRC01, CAP256-VRC26-25, PGDM1400 and PGT121) that are not only the broadest and most potent amongst the best bnAbs discovered with distinct epitope specificities on viral envelope to date but also amongst ones that are farthest in under clinical development and testing [29,30,31].
Evidence of variation in susceptibilities of HIV-1 India clade C to bnAbs targeting gp120 having distinct target epitope specificities
We first examined the degree of neutralization sensitivity of HIV-1 Indian clade C to four bnAbs: VRC01 (having CD4-binding site specificity), PGDM1400 and CAP256-VRC26.25 (both having V1V2 specificity) and PGT121 (target V3 supersite). A total of 71 pseudoviruses expressing complete env gene (gp160) obtained from 40 disease stage specific samples from 28 HIV-1 infected individuals from Southern and Western India during the years 1999–2014 (the complete details of these env clones along with their genetic properties are given in Additional file 1: Table S1) were assessed for their degree of neutralization sensitivity towards four selected bNAbs of distinct specificities in TZM-bl neutralization assay. Several of these envs have been reported earlier [32,33,34,35,36,37,38]. As shown in Table 1, the degree of susceptibility of Env-pseudotyped viruses to the bnAbs varied. The neutralization data demonstrated that while VRC01 exhibited most breadth (neutralized 62/71 viruses), PGDM1400 exhibited least (46/71) of all the four bnAbs assessed. The neutralization coverage was assessed by measuring the ability of each of the bnAb to demonstrate 50% virus neutralization by IgG concentration less than or equal to 5 µg/ml (IC50). In addition, CAP256-VRC26.25 was found to neutralize the virus panel with maximal potency (IC50 of 0.0033 with 76% breadth). Interestingly, while neutralization resistance of Env-pseudotyped viruses to all four bnAbs was found to be associated with known resistance signatures (Additional file 2: Table S2), for a few envs, no known resistance signatures were evident, suggesting an alternate mechanism of neutralization resistance. Overall, our data showed variation in neutralization breadth and potency of VRC01, CAP256-VRC26.25, PGDM1400 and PGT121 to HIV-1 Indian clade C envelopes that were tested in this study.
Comparing neutralization diversity of Indian clade C with other globally circulating subtypes and recombinants
We next examined whether inter and intra clade env diversity have any association with altered neutralization phenotype. For this first we retrieved full length gp160 sequences of Env-pseudotyped viruses reported in the CATNAP database (www.hiv.lanl.gov) and compared diversity by building a phylogenetic tree (Fig. 1a). As expected, our data very clearly demarcated the distinctness of Indian clade C with non-Indian clade C and other subtypes (Fig. 1a; top left). This observation also corroborates with our earlier finding . The neutralization scores (IC50 values) against VRC01, CAP256-VRC26.25, PGDM1400 and PGT121 of most of the same retrieved gp160 sequences of all subtypes were next obtained from the CATNAP database for the purpose of comparing with that obtained in the current study. A total of 1020 gp160 sequences were used for the analysis with the following distribution: Subtype C (Pan Africa):290, Subtype C (India): 98 that includes 71 used in the present study, Subtype A: 76, Subtype B: 255, Subtype D: 42, CRF01_AE: 70, CRF07_B/C recombinants: 49, other subtypes and recombinant forms: 140 were assessed for their phylogenetic relatedness using IQtree (HIVdb model, non-parametric 1000 fast bootstrap with aSH-LRT test). Four heatmaps based on their responses to VRC01, CAP256-VRC26.25, PGDM1400 and PGT121 were built based on the phylogenetic tree to reflect the clustering based on their IC50 values obtained against each of these Env-pseudotyped viruses. As shown in Fig. 1a, We observed that while VRC01 appeared to be most broad, it was found to be least potent amongst all. CAP256-VRC26.25 appeared to be most potent, however our data also indicates that it has been mostly assessed against HIV-1 subtype C (including one used in the present study). These qualitative observations also reflected upon statistical assessment of viruses from India as indicated in Fig. 1b.
Finally, we compared the overall neutralization sensitivity against these four bnAbs of year matched (1999–2011), randomly selected an equal number of viral gp160 datasets (n = 85) from CATNAP database with the ones assessed in this current study. The IC50 value of 5 µg/ml was considered as a threshold of neutralization sensitivity. As shown in Fig. 1c, the neutralization susceptibility of Indian clade C envelopes to VRC01 were found to be distinct from that observed with the African clade C as well as other subtypes. Interestingly, while the degree of CAP256-VRC26.25 susceptibility between Indian clade C and other subtypes were observed to be comparable, this was found to be significantly different from African clade C. Overall, our data highlights disparate sensitivity of Indian clade C compared to other geographically and phylogenetically distinct HIV-1 subtypes against these four bnAbs.
Trend analysis of variation in neutralization sensitivity over time at population and individual levels
We next examined the variation in neutralization sensitivity of HIV-1 env sequences sampled overtime to these four bnAbs. For this, we grouped IC50 values (in three clusters) based on corresponding year of sampling (1990–2000, 2001–2010 and 2011–2016) of the viral env sequences. As shown in Fig. 2, the PGDM1400 sensitivity of Env-pseudotyped viruses obtained from both Indian and African population was observed to be decreasing over time. For PGT121 and VRC01 bnAbs, a similar trend was observed in the Indian population. When compared, a gradual increase in resistant neutralization phenotype of African subtypes to PGT121 and VRC01 was also observed over time, it was, however, not found to be significant. Interestingly, an increasing trend in neutralization resistance of Indian clade C envelopes to CAP256-VRC26.25 was observed in contrast to that with envelopes obtained from the African population. Overall, we observed that although a number of Indian clade C envelopes tested in this study were not large enough for an absolute conclusion, our data clearly indicated differences in neutralization susceptibility trend against these four bnAbs overtime between Indian clade C and other subtypes.
We next examined the trend of neutralization susceptibility of pseudoviruses expressing primary HIV-1 clade C envs obtained from three individuals at different time points to VRC01, CAP256-VRC26.25, PGDM1400 and PGT121. Several of these envs were reported earlier by our group . The primary envs (complete gp160 sequences) from the donor NARI IVC-2 were obtained at following time point: between 0–6, 12, 18, 24, 36 and 57 months; from donors NARI IVC-3 and NARI IVC-4 at 0–6, 12 and 18 months respectively. As demonstrated in Fig. 3, the trend of susceptibility of pseudoviruses expressing envs isolated from these three donors overtime to VRC01, CAP256-VRC26.25, PGDM1400 and PGT121 varied. Out of the three donors, in one (NARI IVC-11), a clear trend of increasing susceptibility of Env-pseudotyped viruses to all the four bnAbs were observed. When tested against PGT121, for NARI IVC-2, pseudoviruses expressing envs obtained from baseline to two years, were found to be resistant, while pseudoviruses expressing envs obtained from subsequent time points till 57 months were found to be highly sensitive with IC50 over 0.01 µg/ml. For all the resistant envs, N332 glycan residue was found to be absent. Interestingly, we observed that the env (2–11. J16) obtained from 57-month time point was sensitive to PGT121 despite lacking the N332 glycan. While this could possibly be compensated by the presence of a glycan at the 334 position (Additional file 3: Figure S1) as was reported for PGT128 with similar specificity , four pseudoviruses expressing envs isolated from NARI IVC-2 were found to be resistant to PGT121 despite containing N334. While pseudoviruses expressing envs obtained from all time points from this donor were found to be highly sensitive to CAP256-VRC26.25 with IC50 < 0.001 µg/ml, no increasing or decreasing trend in susceptibility to VRC01 and PGDM1400 of all the envelopes were observed. For the third donor, NARI IVC-3, envelopes obtained from all the time points were found to be resistant to V1V2 -specific bnAbs CAP256-VRC26.25 and PGDM1400, while they were found to be sensitive to VRC01 and PGT121, however no clear trend in neutralization sensitivity was observed over time. While resistance to different bnAbs was found to be associated with an absence of known target motifs, we found evidence of neutralization resistance in presence of known epitopes that are targets of some bnAbs (Additional file 2: Table S2). Overall, we found that while increasing neutralization resistance trend was observed at the population level, when cross-sectional virus isolates across years were analyzed such a trend was not evident with longitudinal viral envs obtained from three followed up Indian donors.
Predicting optimal bnAb combination using experimental data
Finally, based on the neutralization data of all the Indian clade C envelopes assessed in this study and neutralization scores of additional Indian clade C envs retrieved from CATNAP database (http://www.hiv.lanl.gov), we attempted to predict the optimal combination amongst these four bnAbs that would likely confer maximal neutralization coverage with the highest potency. Towards achieving this, we employed an additive model implemented in the CombiNAber tool (https://www.hiv.lanl.gov/content/sequence/COMBINABER/combinaber.html) to all possible combinations of these four bnAbs (VRC01, CAP256-VRC26.25, PGDM1400 and PGT121). A total of 98 Indian clade C envelopes (71 from the current study and existing neutralization data of 27 additional Indian clade C retrieved from HIV database) for preparing for coverage analysis. As shown in Fig. 4a, the cumulative coverage analysis indicated that each of VRC01, CAP256-VRC26.25, PGDM1400 and PGT121 bnAbs individually could provide 80%, 70%, 65% and 74%coverage respectively at IgG concentration of equal to or below 5 µg/ml. When assessed for the minimum combination that would provide maximal neutralization coverage, two combinations of three bnAbs each (VRC01 + CAP56-VRC26.25 + PGT121 and PGDM1400 + CAP256-VRC26.25 + PGT121) appears to provide 100% cumulative coverage at total IgG concentration of 5 µg/ml or less as predicted by CombiNAber. In parallel, we also examined whether these two bnAb combinations of three bnAbs could achieve the maximal neutralization coverage at low doses. Towards this, the median IC50 values of Indian clade C envelopes (n = 98) obtained against single and different combinations of all the four bnAbs were obtained and the coverage versus potency assessment was statistically validated by Mann–Whitney test on the dataset obtained by implementation of the additive model by the CombiNAber tool with a target limit of IgG concentration of 5 µg/ml. As shown in Fig. 4b, we observed that a combination of CAP256-VRC26.25 + PGDM1400 and PGT121 could achieve near 100% breadth with the highest potency (Fig. 4b). This was found to be strikingly different from what we observed with African subtype C, where two combinations of each of these four bnAbs appeared to provide near 100% neutralization coverage with maximal potency. For all the different combination we assessed, we compared the median IC50 obtained in our study with that reported for African clade C envelopes (n = 250) by retrieving neutralization IC50 values from CATNAP database (http://www.hiv.lanl.gov). As expected, we found that the same combinations of three bnAbs (VRC01 + CAP256-VRC26.25 + PGT121 and PGDM1400 + CAP256-VRC26.25 + PGT121) demonstrated maximal coverage with geometric cumulative mean IC50 of 0.0109 and 0.01137 µg/ml respectively and are observed to be comparable to what was observed with the combination of all four bnAbs. Of interest, when compared with African clade C envelopes, we observed that a combination of two bnAbs that commonly include CAP256-VRC26.25 (such as CAP256-VRC26.25 + PGDM1400, CAP256-VRC26.25 + VRC01 and CAP256-VRC26.25 + PGT121) appear to provide comparable coverage with significantly lower doses than what we observed with the three bnAb combinations with Indian clade C viruses (Fig. 4b). This is likely because CAP256-VRC26.25 was found to demonstrate better neutralization coverage of African clade C with increased potency compared to what we observed with Indian clade C viruses.
Understudied and unexplored population and clade-specific HIV-1 sequence diversity presents a formidable challenge that is likely to play a major hindrance in designing of effective antibody-mediated intervention strategies. Subtype C being the most predominant subtype circulating majorly in Africa and India, a greater understanding of the antigenic and neutralization diversity of HIV-1 subtype C env would facilitate understanding of optimal bnAb combination that could potentially overcome region-specific the intra-clade C diversity. Our recent work has demonstrated genotypic distinctness of HIV-1 subtype C envelope (gp120) sequences reported from Indian and African epidemics . In that study, we observed these sequences to be distinct in terms of variable loop lengths (V1, V2 and V4), an abundance of potential N-linked glycosylation sites and entropy at bnAb contact sites (including PGDM1400 and CAP256-VRC26.25). Through prediction algorithms, we explored the impact of these diverging signatures on sensitivity to bnAbs and predicted appreciable differences. Through this present study, we attempted to experimentally elucidate the distinct neutralization phenotype of primary HIV-1 subtype C obtained from Indian patients to four key bnAbs that not only have demonstrated considerable breadth, but also amongst those which are farthest into clinical development and testing. When analyzed against pseudoviruses expressing 71 primary envs, we observed that while VRC01 demonstrated maximum breadth, PGDM1400 showed the least amongst the four. Conversely, CAP256-VRC26.25 was found to neutralize over 76% viruses with maximal potency amongst others. Interestingly, while the neutralization resistance of the viruses to PGDM1400 was majorly found to be associated with known substitutions as found in the HIV database (http://www.hiv.lanl.gov), for some resistant envelopes, no known resistant motifs were found, indicating plausible underlying unknown mechanism associated with PGDM1400 resistance. Furthermore, while CAP256-VRC26.25 and PGDM1400 target epitopes on the V1V2 loop, Indian clade C envelopes were found to be more susceptible to CAP256-VRC26.25 with fewer resistant signatures compared to PGDM1400. This is an interesting and perhaps an important observation which warrants the need to screen for additional significantly large population-based study to understand the frequency of occurrence of signature motifs and/or substitutions associated with neutralization resistance. This is particularly important while choosing bnAbs targeting V1V2 region for clinical use as V1V2 region represents the most hypervariable region and subtle changes in loop length, glycosylation patterns and substitutions will compromise with the antibody-mediated virus neutralization efficacy . Indeed, a significant difference in neutralization breadth and potency was observed when we compared the sensitivity of Indian subtype C with that of African subtype C based on the existing data. Interestingly though, while comparable sensitivity to PGDM1400 was observed between Indian and African subtype C viruses, significant differences were observed when sensitivities to VRC01, CAP256-VRC26.25 and PGT121 were assessed. This observation further highlights the variation in population and geography based env genetic diversity associated with neutralization diversities.
Trend analysis at the population level with limited data sets indicated a sense of increasing resistance to all the four bnAbs tested in this study with envelopes isolated from Indian patients between the years 2011–2016. This could possibly be linked to the differences in host characteristics and varied genetic bottlenecks which contribute to altered virus and antibody evolution pathways in the course of natural infection. Our observation while found reasonably comparable to what was observed with African subtype C, except for CAP256-VRC26.25, we noted that most if not all the African C envelopes were obtained from acute/early infection. In this line, given this was an observational study, we also hypothesize that the stage of the disease, origin of the env, and the possible geographic differences (within and outside India) could potentially influence on the neutralization phenotypes to draw any robust conclusion. At the individual level, while an increasing susceptibility trend was observed with PGT121, no consistent trend was observed with other three bnAbs. While we tested envs obtained from follow up visits of only three individuals, it would be important to study envs obtained from more number of individuals in the same follow-up setting. Such exercise would possibly help understand optimal combination of bnAbs that would be suitable for prophylaxis or early treatment versus the combination that would be more suitable to treat at the late disease stage.
One of the approaches that has widely been considered to provide maximal neutralization breadth at low dose and overcoming virus escape (particularly in the context of antibody-mediated treatment) is a combination of bnAbs having distinct epitope specificities. Based on our experimental data, we observed that combination of two V1V2 directed (CAP256-VRC26.25, PGDM1400) and one V3 supersite-directed (PGT121) bnAbs could achieve near 100% neutralization coverage with maximal potency. Intriguingly, this observation was found to be distinct to what we observed with African subtype C where combinations of two bnAbs could possibly provide near 100% neutralization coverage at the highest potency. Our observation further highlighted neutralization diversity between Indian and African subtype C. It is to be noted that while the majority of the African clade C env sequences analyzed in this study were derived from acute and early infection cases (http://www.hiv.lan.gov), majority of the Indian subtype C env sequences used in the present study were obtained from individuals post 6 months of infection and at late disease stages. Assessment of late disease stages from both regions would be important to design bnAb based therapeutic strategies. While not many Indian subtype C env sequences from acute/early infection are available, we believe that it will be very important to obtain such sequences that will be necessary to define, select and optimize bnAbs and their combination for developing antibody-mediated prevention strategies. Several factors might account for this variability. As a result of the continuous drift of the HIV-1 genome and epitopes, it is possible that divergent subtypes follow separate evolution pathways that may impact their degree of susceptibility to different bnAbs. This may also be due to the population-specific diverse and distinct selection pressure. Our hypothesis is supported by the fact that different classes of bnAbs isolated from HIV-1 infected individuals were influenced by both the infecting subtypes  and also by the ethnic origin [24, 42]. In addition, the possibility of increasing resistance to bnAbs as observed [25, 26] in recent studies over the course of the epidemic warrants the necessity for continued surveillance of virus evolution that would impact the ability of the promising bnAbs mediated neutralization coverage of region-specific HIV-1 diversities.
In conclusion, our data indicated disparate neutralization diversities of HIV-1 subtype C when assessed against four bnAbs with distinct epitope specificities and those that are amongst the furthest into the clinical development and testing. This observation necessitates the urgent need to compare a large number of env sequences across different geographies, risk groups and key populations to understand how disparate intra-clade genetic diversity could potentially impact virus neutralization. This will inform the choice of the bnAb combinations that would be suitable for optimal coverage of the region-specific HIV-1 circulating subtypes for prevention and treatment over and above the existing antiretroviral drugs.
The plasma and PBMC samples used for amplification of complete env genes by PCR were collected following respective institutional ethical approvals.
Plasmids, antibodies and cell lines
Plasmid DNA encoding HIV-1 gp160 sequences published earlier are listed in Table 1. Plasmid DNA encoding heavy and light variable IgG chain sequences of the following monoclonal antibodies: PGT121, PGDM1400 and VRC01 were obtained from the IAVI neutralizing antibody center at the Scripps Research and CAP256.VRC26.25 was kindly provided by Prof Lynn Morris, National Institute of Communicable Diseases (NICD), Johannesburg, South Africa under an ongoing collaborative project. TZM-bl cells were obtained from the NIH AIDS Research Program. HEK-293 T cell line was obtained from American Type Culture Collection (ATCC).
Amplification and cloning of gp160
Full length env (gp160) genes were PCR amplified from HIV + plasma samples with slight modification as described previously . Briefly, viral RNA was extracted using a high-pure viral RNA kit (Roche Inc.) by following manufacturer’s protocol, and cDNA was prepared by reverse transcription-PCR (RT-PCR) using a Superscript III first-strand synthesis kit (Invitrogen Inc.). The primer used for cDNA synthesis was EnvR1 5′-GCACTCAAGGCAAGCTTTATTGAGGCTT -3′ (HXB2: 9605–9632). Rev-Env gp160 cassette were amplified from the cDNA product using La Taq high fidelity DNA polymerase in the 1st round (Takara Bio Inc.) and PrimeSTAR GXL high fidelity DNA polymerase (Takara Bio Inc.) in the second round. The primers used for the 1st round were EnvF1: 5′- AGARGAYAGATGGAACAAGCCCCAG-3′ (HXB2: 5550–5574) and EnvRP2: 5′-GTGTGTAGTTCTGCCAATCAGGGAA -3′ (HXB2: 9157–9181) while for the second round were Env IF: 5′-CACCGGCTTAGGCATCTCCTATGGCAGGAAGAA -3′ (HXB2: 5950–5982) and Env IR: 5′-TATCGGTACCAGTCTTGAGACGCTGCTCCTACTC -3′ (HXB2: 8882–8915). PCR condition followed in the first round was initial denaturation of 94 °C for 2 min followed by 12 cycles of 94 °C for10 secs, 60 °C for 30 secs, 68 °C for 3 min, 23 cycles of 94 °C for 10 secs, 55 °C for 30 s, 68 °C for 3 min with final extension of 68 °C for 10 min. PCR condition followed in the second round was initial denaturation of 94 °C for 2 min followed by 12 cycles of 94 °C for 10 secs, 62 °C for 30 s, 68 °C for 3 min, 23 cycles of 94 °C for 10 secs, 60 °C for 30secs, 68 °C for 3 min with a final extension of 68 °C for 10 min. The gp160 amplicons were purified and ligated into pcDNA 3.1/V5-His-TOPO (Invitrogen Inc.) vector or pSVIII as described before .
The gp160 amino acid sequences generated in the present study, those published in the GenBank earlier as well as those derived from the CATNAP database (supplementary information) were aligned to each other using Muscle v3.8.1551. The resulting alignment was curated using BioEdit sequence alignment editor v7.2.5 . Variable loop regions for V1 (131–157: HXB2 numbering), V2 (158–196), V3 (296–331), V4 (386–417) and V5 (460–469) were retrieved from amino acid alignments with Extractalign implementation of the Emboss package . Each of the loop datasets were then processed with custom bash/awk scripts to generate length statistics. Potential N-linked glycosylation sites were predicted using N-Glycosite tool available at the LANL HIV database . Statistics regarding the antibody resistance-associated signature residues were generated as described recently . A phylogenetic tree was generated for all the HIV-1 viral clone sequences available in the CATNAP database along with those generated in the present study (Total N = 1020) with iqtree under ‘HIVb’ model with estimated Ƴ parameters and number of invariable sites . The robustness of the tree topology was further assessed by SH-aLRT as well as 1000 ultrafast bootstrap replicates implemented in iqtree. The tree as well as neutralization IC50 (µg/ML) values for each of the sequences against bnAbs PGDM1400, PGT121, VRC01 and CAP256-VRC26.25 were plotted with ‘ggtree’ package in R [47, 48]. IC50 values for year matched randomly selected equal number (N = 85) were compared for HIV-1 clade C viral clones reported from India (including this study), Africa as well as other subtypes using Mann–Whitney test. Trend analysis for change in viral sensitivity to bnAbs VRC01, PGDM1400, CAP256-VRC26.25 and PGT121 was further performed. To achieve this, IC50 values for all Indian viral clones were grouped into three-time periods as per their year of sampling (1990–2000, 2001–2010 and 2011–2016) and compared with Jonckheere Terpstra test (JTT). CombiNAber (http://www.hiv.lanl.gov) , was used to predict optimal antibody combination for maximum neutralization of viral clones from India (N = 98) and from Africa (N-250) using their experimentally determined IC50 values. Target IgG concentration of 5 µg/mL was considered as a threshold followed by implementation of the additive model in the CombiNAber tool to obtain geometric mean IC50 values for each combination. Sensitivity to each of these bnAbs combinations was further compared between Indian and African viruses using Mann–Whitney test. statistical computing software (v3.4.0) and R studio v1.0.143 [49, 50]. Statistical analysis for neutralization breadth and potency were done using GraphPad Prism version 5.01 for Windows, GraphPad Software, San Diego California USA. Scatter plots for bnAb potency distribution across geographically distinct populations as well as those from longitudinally collected samples were assessed statistically and plotted in GraphPad Prism v5.01.
Preparation of Env-pseudoviruses
Pseudotyped viruses were prepared by the co-transfection of env-expressing plasmid DNA along with the plasmid DNA expressing HIV-1 genes with premature stop codon for env (pSG3ΔEnv) into 293 T cells in 6-well tissue culture plates using FuGENE6 transfection reagent kit (Promega Inc.). Cell supernatants containing pseudotyped viruses were harvested at 48 h post transfection and subsequently stored at -80° C until use. The virus infectivity was measured using TZM-bl reporter cells by addition of pseudoviruses containing DEAE-dextran (25 µg/ml) in 96-well microtiter plates, and the viral titers were determined by measuring the luciferase activity using Britelite luciferase substrate (PerkinElmer Inc.) with a Victor X2 luminometer (PerkinElmer Inc.).
Neutralization assays were carried out using TZM-bl cells as described before . Briefly, Env-pseudotyped viruses were pre-incubated in 96-well tissue culture plates with various concentrations of bnAbs (IgG) for an hour at 37° C in a CO2 incubator under humidified conditions. Subsequently, 1 × 104 TZM-bl cells were added to the mixture in the presence of 25 µg/ml DEAE-dextran (Sigma, Inc.). The plates were further incubated for 48 h. The degree of virus neutralization was assessed by measuring reduction in relative luminescence units (RLU) in a luminometer (Victor X2; PerkinElmer Inc.).
Availability of data and materials
The dataset(s) supporting the conclusions of this article is (are) included within the article (and its additional file(s). The complete nucleotide sequences of the following envs: 4.2.J45b, 2–7.J1, 2–9.J20, 2–11.J16, NISA-N20-J10, NISA-N20.J14, NISA-N101.J12, INDO SA NLR 29.J80, INDO SA NLR 29.J11, PG37009v2.eJ9, PG37009v2.eJ38, PG37009v2.eJ58, PG37112v2.J5, PG37112v2.J9, PG37072.J12, PG37072.J16, PG37066.J1, PG37081.J36, PG37087.J39, PG37087.J44, PG37089.J17, PG37089.J20, PG37089.J83, PG37091.J41, PG37080.J6A, PG37080v1.J17, PG37080.J158 are being submitted to the GenBank.
Human Immunodeficiency Virus
Broadly neutralizing monoclonal antibodies
Shimodaira-Hasegawa approximate likelihood ratio test
Compile, Analyze and Tally NAb Panels
Stefic K, Bouvin-Pley M, Braibant M, Barin F. Impact of HIV-1 diversity on its sensitivity to neutralization. Vaccines (Basel). 2019;7:3.
Sok D, Burton DR. Recent progress in broadly neutralizing antibodies to HIV. Nat Immunol. 2018;19(11):1179–88.
Sok D, Le KM, Vadnais M, Saye-Francisco KL, Jardine JG, Torres JL, Berndsen ZT, Kong L, Stanfield R, Ruiz J, et al. Rapid elicitation of broadly neutralizing antibodies to HIV by immunization in cows. Nature. 2017;548(7665):108–11.
Walker LM, Huber M, Doores KJ, Falkowska E, Pejchal R, Julien JP, Wang SK, Ramos A, Chan-Hui PY, Moyle M, et al. Broad neutralization coverage of HIV by multiple highly potent antibodies. Nature. 2011;477(7365):466–70.
Walker LM, Phogat SK, Chan-Hui PY, Wagner D, Phung P, Goss JL, Wrin T, Simek MD, Fling S, Mitcham JL, et al. Broad and potent neutralizing antibodies from an African donor reveal a new HIV-1 vaccine target. Science. 2009;326(5950):285–9.
Doria-Rose NA, Bhiman JN, Roark RS, Schramm CA, Gorman J, Chuang GY, Pancera M, Cale EM, Ernandes MJ, Louder MK, et al. New Member of the V1V2-Directed CAP256-VRC26 lineage that shows increased breadth and exceptional potency. J Virol. 2016;90(1):76–91.
Wu X, Yang ZY, Li Y, Hogerkorp CM, Schief WR, Seaman MS, Zhou T, Schmidt SD, Wu L, Xu L, et al. Rational design of envelope identifies broadly neutralizing human monoclonal antibodies to HIV-1. Science. 2010;329(5993):856–61.
Kumar R, Qureshi H, Deshpande S, Bhattacharya J. Broadly neutralizing antibodies in HIV-1 treatment and prevention. Ther Adv Vacc Imunother. 2018;6(4):61–8.
Sok D, van Gils MJ, Pauthner M, Julien JP, Saye-Francisco KL, Hsueh J, Briney B, Lee JH, Le KM, Lee PS, et al. Recombinant HIV envelope trimer selects for quaternary-dependent antibodies targeting the trimer apex. Proc Natl Acad Sci USA. 2014;111(49):17624–9.
Liu Y, Cao W, Sun M, Li T. Broadly neutralizing antibodies for HIV-1: efficacies, challenges and opportunities. Emerg Microb Infect. 2020;9(1):194–206.
Caskey M, Klein F, Lorenzi JC, Seaman MS, West AP, Jr., Buckley N, Kremer G, Nogueira L, Braunschweig M, Scheid JF et al: Viraemia suppressed in HIV-1-infected humans by broadly neutralizing antibody 3BNC117. Nature 2015.
Caskey M, Klein F, Nussenzweig MC. Broadly neutralizing anti-HIV-1 monoclonal antibodies in the clinic. Nat Med. 2019;25(4):547–53.
Caskey M, Schoofs T, Gruell H, Settler A, Karagounis T, Kreider EF, Murrell B, Pfeifer N, Nogueira L, Oliveira TY, et al. Antibody 10–1074 suppresses viremia in HIV-1-infected individuals. Nat Med. 2017;23(2):185–91.
Mendoza P, Gruell H, Nogueira L, Pai JA, Butler AL, Millard K, Lehmann C, Suarez I, Oliveira TY, Lorenzi JCC, et al. Combination therapy with anti-HIV-1 antibodies maintains viral suppression. Nature. 2018;561(7724):479–84.
Bar-On Y, Gruell H, Schoofs T, Pai JA, Nogueira L, Butler AL, Millard K, Lehmann C, Suarez I, Oliveira TY, et al. Safety and antiviral activity of combination HIV-1 broadly neutralizing antibodies in viremic individuals. Nat Med. 2018;24(11):1701–7.
Scheid JF, Horwitz JA, Bar-On Y, Kreider EF, Lu CL, Lorenzi JC, Feldmann A, Braunschweig M, Nogueira L, Oliveira T, et al. HIV-1 antibody 3BNC117 suppresses viral rebound in humans during treatment interruption. Nature. 2016;535(7613):556–60.
Binley JM, Wrin T, Korber B, et al. Comprehensive cross-clade neutralization analysis of a panel of anti-human immunodeficiency virus type 1 monoclonal antibodies. J Virol. 2004;78:13232–52.
Bures R, Morris L, Williamson C, Ramjee G, Deers M, Fiscus SA, Abdool-Karim S, Montefiori DC. Regional clustering of shared neutralization determinants on primary isolates of clade C human immunodeficiency virus type 1 from South Africa. J Virol. 2002;76(5):2233–44.
Hraber P, Korber BT, Lapedes AS, Bailer RT, Seaman MS, Gao H, Greene KM, McCutchan F, Williamson C, Kim JH, et al. Impact of clade, geography, and age of the epidemic on HIV-1 neutralization by antibodies. J Virol. 2014;88(21):12623–43.
Kulkarni SS, Lapedes A, Tang H, Gnanakaran S, Daniels MG, Zhang M, Bhattacharya T, Li M, Polonis VR, McCutchan FE, et al. Highly complex neutralization determinants on a monophyletic lineage of newly transmitted subtype C HIV-1 Env clones from India. Virology. 2009;385(2):505–20.
Seaman MS, Janes H, Hawkins N, Grandpre LE, Devoy C, Giri A, Coffey RT, Harris L, Wood B, Daniels MG, et al. Tiered categorization of a diverse panel of HIV-1 Env pseudoviruses for assessment of neutralizing antibodies. J Virol. 2010;84(3):1439–52.
Burnie J, Guzzo C. The Incorporation of Host Proteins into the External HIV-1 Envelope. Viruses. 2019;11:1.
Rusert P, Kouyos RD, Kadelka C, Ebner H, Schanz M, Huber M, Braun DL, Hoze N, Scherrer A, Magnus C, et al. Determinants of HIV-1 broadly neutralizing antibody induction. Nat Med. 2016;22(11):1260–7.
Novitsky V, Smith UR, Gilbert P, McLane MF, Chigwedere P, Williamson C, Ndung’u T, Klein I, Chang SY, Peter T, et al. Human immunodeficiency virus type 1 subtype C molecular phylogeny: consensus sequence for an AIDS vaccine design? J Virol. 2002;76(11):5435–51.
Wagh K, Bhattacharya T, Williamson C, Robles A, Bayne M, Garrity J, Rist M, Rademeyer C, Yoon H, Lapedes A, et al. Optimal Combinations of Broadly Neutralizing Antibodies for Prevention and Treatment of HIV-1 Clade C Infection. PLoS Pathog. 2016;12(3):e1005520.
Lorenzi JCC, Mendoza P, Cohen YZ, Nogueira L, Lavine C, Sapiente J, Wiatr M, Mugo NR, Mujugira A, Delany S et al: Neutralizing Activity of Broadly Neutralizing anti-HIV-1 Antibodies against Primary African Isolates. J Virol 2020.
Bai H, Li Y, Michael NL, Robb ML, Rolland M. The breadth of HIV-1 neutralizing antibodies depends on the conservation of key sites in their epitopes. PLoS Comput Biol. 2019;15(6):e1007056.
Rademeyer C, Korber B, Seaman MS, Giorgi EE, Thebus R, Robles A, Sheward DJ, Wagh K, Garrity J, Carey BR, et al. Features of recently transmitted HIV-1 Clade C viruses that impact antibody recognition: implications for active and passive immunization. PLoS Pathog. 2016;12(7):e1005742.
Crowell TA, Colby DJ, Pinyakorn S, Sacdalan C, Pagliuzza A, Intasan J, Benjapornpong K, Tangnaree K, Chomchey N, Kroon E, et al. Safety and efficacy of VRC01 broadly neutralising antibodies in adults with acutely treated HIV (RV397): a phase 2, randomised, double-blind, placebo-controlled trial. The lancet HIV. 2019;6(5):e297–306.
Reeves DB, Huang Y, Duke ER, Mayer BT, Cardozo-Ojeda EF, Boshier FA, Swan DA, Rolland M, Robb ML, Mascola JR, et al. Mathematical modeling to reveal breakthrough mechanisms in the HIV Antibody Mediated Prevention (AMP) trials. PLoS Comput Biol. 2020;16(2):e1007626.
Mahomed S, Garrett N, Karim QA, Zuma NY, Capparelli E, Baxter C, Gengiah T, Archary D, Samsunder N, Rose ND, et al. Assessing the safety and pharmacokinetics of the anti-HIV monoclonal antibody CAP256V2LS alone and in combination with VRC07-523LS and PGT121 in South African women: study protocol for the first-in-human CAPRISA 012B phase I clinical trial. BMJ Open. 2020;10(11):e042247.
Gharu L, Ringe R, Bhattacharya J. HIV-1 clade C envelopes obtained from late stage symptomatic Indian patients varied in their ability towards relative CD4 usages and sensitivity to CCR5 antagonist TAK-779. Virus Res. 2011;158:216–24.
Gharu L, Ringe R, Bhattacharya J. Evidence of extended alternate coreceptor usage by HIV-1 clade C envelope obtained from an Indian patient. Virus Res. 2012;163(1):410–4.
Gharu L, Ringe R, Pandey S, Paranjape R, Bhattacharya J. HIV-1 clade C env clones obtained from an Indian patient exhibiting expanded coreceptor tropism are presented with naturally occurring unusual amino acid substitutions in V3 loop. Virus Res. 2009;144(1–2):306–14.
Ringe R, Thakar M, Bhattacharya J. Variations in autologous neutralization and CD4 dependence of b12 resistant HIV-1 clade C env clones obtained at different time points from antiretroviral naive Indian patients with recent infection. Retrovirology. 2010;7(1):76.
Mukhopadhyay S, Ringe R, Patil A, Paranjape R, Bhattacharya J. Characterization of Circulating HIV Type 1 env Genes in Plasma of Two Antiretroviral-Naive Slow Progressing Patients with Broad Neutralizing Antibody Response with Evidence of Recombination. AIDS Res Hum Retroviruses. 2012;28(7):739–45.
Patil S, Kumar R, Deshpande S, Samal S, Shrivastava T, Boliar S, Bansal M, Chaudhary NK, Srikrishnan AK, Murugavel KG, et al. Conformational epitope-specific broadly neutralizing plasma antibodies obtained from an HIV-1 Clade C-infected elite neutralizer mediate autologous virus escape through mutations in the V1 Loop. J Virol. 2016;90(7):3446–57.
Deshpande S, Patil S, Kumar R, Hermanus T, Murugavel KG, Srikrishnan AK, Solomon S, Morris L, Bhattacharya J. HIV-1 clade C escapes broadly neutralizing autologous antibodies with N332 glycan specificity by distinct mechanisms. Retrovirology. 2016;13(1):60.
Sutar J, Deshpande S, Mullick R, Hingankar N, Patel V, Bhattacharya J. Geospatial HIV-1 subtype C gp120 sequence diversity and its predicted impact on broadly neutralizing antibody sensitivity. BioRxiv. 2020. https://doi.org/10.1101/2020.09.09.289132.
Krumm SA, Mohammed H, Le KM, Crispin M, Wrin T, Poignard P, Burton DR, Doores KJ. Mechanisms of escape from the PGT128 family of anti-HIV broadly neutralizing antibodies. Retrovirology. 2016;13(1):8.
Ringe R, Phogat S, Bhattacharya J. Subtle alteration of residues including N-linked glycans in V2 loop modulate HIV-1 neutralization by PG9 and PG16 monoclonal antibodies. Virology. 2012;426(1):34–41.
Kist NC, Lambert B, Campbell S, Katzourakis A, Lunn D, Lemey P, Iversen AKN. HIV-1 p24Gag adaptation to modern and archaic HLA-allele frequency differences in ethnic groups contributes to viral subtype diversification. Virus Evol. 2020;6(2):085.
Hall T: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp 1999().
Rice P, Longden I, Bleasby A. EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000;16(6):276–7.
Zhang M, Gaschen B, Blay W, Foley B, Haigwood N, Kuiken C, Korber B. Tracking global patterns of N-linked glycosylation site variation in highly variable viral glycoproteins: HIV, SIV, and HCV envelopes and influenza hemagglutinin. Glycobiology. 2004;14(12):1229–46.
Nguyen LT, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32(1):268–74.
Yu G, Lam TT, Zhu H, Guan Y. Two methods for mapping and visualizing associated data on phylogeny using Ggtree. Mol Biol Evol. 2018;35(12):3041–3.
Yu G, Smith D, Zhu H, Guan Y, Lam TT. ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2017;8:28–36.
Team R: RStudio: Integrated Development for R. 2015.
Team RC: R: A language and environment for statistical computing. 2018.
We thank Prof Lynn Morris, National Institute for Communicable Diseases, Johannesburg, South Africa for providing us with the plasmids to produce CAP256-VRC26.25 bnAbs used in this study through an ongoing collaborative project supported by the Department of Biotechnology, Government of India (BT/PR24520/MED/29/1222/2017). Clinical samples collected under the IAVI Protocol G study was supported by the USAID. We thank the Protocol G study participants registered with YRG Care, Chennai, all of the research staff members at the Protocol G clinical center at YRG Care, Chennai, and all of the IAVI Protocol G team members. IAVI’s work was made possible by generous support from many donors, including the Bill & Melinda Gates Foundation, the Ministry of Foreign Affairs of Denmark, Irish Aid, the Ministry of Finance of Japan, the Ministry of Foreign Affairs of the Netherlands, the Norwegian Agency for Development Cooperation (NORAD), the United Kingdom Department for International Development (DFID), and the United States Agency for International Development (USAID). The full list of IAVI donors is available at http://www.iavi.org. The contents are the responsibility of the International AIDS Vaccine Initiative and do not necessarily reflect the views of USAID or the United States Government. We thank Prof Gagandeep Kang and Dr Subeer S. Majumdar, Translational Health Science & Technology Institute for support.
The authors wish to acknowledge the funding support from the Wellcome Trust/DBT India Alliance Team Science Grant (IA/TSG/19/1/600019), Science & Engineering Research Board, Department of Science & Technology, Government of India (CRG/2019/002939) and Department of Biotechnology, Government of India (BT/PR24520/MED/29/1222/2017). The funders played no role in the design of the study and collection, analysis, interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
The institutional ethical approvals and participant informed consents were obtained prior to collection of biological specimens.
Consent to publish
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. History of year of collection of clinical samples, disease stages and env genetic properties.
. Mapping key amino acid substitutions associated with resistance to the bnAbs.
Amino acid sequence alignment of V3 region of the autologous HIV-1 env obtained from NARI IVC-2, NARI IVC-3 and NARI IVC-11 collected at different time points in a span of 2006–2011 to denote the frequency of occurrence of N334 glycosylation.
About this article
Cite this article
Mullick, R., Sutar, J., Hingankar, N. et al. Neutralization diversity of HIV-1 Indian subtype C envelopes obtained from cross sectional and followed up individuals against broadly neutralizing monoclonal antibodies having distinct gp120 specificities. Retrovirology 18, 12 (2021). https://doi.org/10.1186/s12977-021-00556-2
- Clade C
- Neutralizing antibodies