Europe PMC

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

Abstract 


Subtype B HIV-1 has been the primary driver of the HIV-1 epidemic in the United States (U.S.) for over forty years and is also a prominent subtype in the Americas, Europe, Australia, the Middle East and North Africa. In this study, the neutralization profiles of contemporary subtype B Envs from the U.S. were assessed to characterize changes in neutralization sensitivities over time. We generated a panel of 30 contemporary pseudoviruses (PSVs) and demonstrated continued diversification of subtype B Env from the 1980s up to 2018. Neutralization sensitivities of the contemporary subtype B PSVs were characterized using 31 neutralizing antibodies (NAbs) and were compared with strains from earlier in the HIV-1 pandemic. A significant reduction in Env neutralization sensitivity was observed for 27 out of 31 NAbs for the contemporary as compared to earlier-decade subtype B PSVs. A decline in neutralization sensitivity was observed across all Env domains; the NAbs that were most potent early in the pandemic suffered the greatest decline in potency over time. A meta-analysis demonstrated this trend across multiple subtypes. As HIV-1 Env diversification continues, changes in Env antigenicity and neutralization sensitivity should continue to be evaluated to inform the development of improved vaccine and antibody products to prevent and treat HIV-1.

Free full text 


Logo of plospathLink to Publisher's site
PLoS Pathog. 2023 Dec; 19(12): e1011780.
PMCID: PMC10727358
PMID: 38055771

Evolution of HIV-1 envelope towards reduced neutralization sensitivity, as demonstrated by contemporary HIV-1 subtype B from the United States

Lindsay Wieczorek, Conceptualization, Data curation, Formal analysis, Project administration, Supervision, Visualization, Writing – original draft, Writing – review & editing, 1 , 2 Eric Sanders-Buell, Data curation, Methodology, Project administration, Supervision, Writing – review & editing, 1 , 2 Michelle Zemil, Data curation, Investigation, Methodology, Project administration, Supervision, Writing – review & editing, 1 , 2 Eric Lewitus, Data curation, Formal analysis, Visualization, Writing – review & editing, 1 , 2 Erin Kavusak, Investigation, Writing – review & editing, 1 , 2 Jonah Heller, Investigation, Writing – review & editing, 1 , 2 Sebastian Molnar, Investigation, Methodology, Writing – review & editing, 1 , 2 Mekhala Rao, Investigation, Writing – review & editing, 1 , 2 Gabriel Smith, Investigation, Writing – review & editing, 1 , 2 Meera Bose, Investigation, Writing – review & editing, 1 , 2 Amy Nguyen, Investigation, 1 , 2 Adwitiya Dhungana, Investigation, Writing – review & editing, 1 , 2 Katherine Okada, Investigation, Writing – review & editing, 1 , 2 Kelly Parisi, Investigation, 1 , 2 Daniel Silas, Investigation, 1 , 2 Bonnie Slike, Investigation, Writing – review & editing, 1 , 2 Anuradha Ganesan, Resources, Writing – review & editing, 2 , 3 , 4 Jason Okulicz, Resources, Writing – review & editing, 3 , 5 Tahaniyat Lalani, Resources, Writing – review & editing, 2 , 3 , 6 Brian K. Agan, Resources, Writing – review & editing, 2 , 3 Trevor A. Crowell, Resources, Writing – review & editing, 1 , 2 Janice Darden, Writing – review & editing, 2 , 7 Morgane Rolland, Supervision, Writing – review & editing, 1 , 2 Sandhya Vasan, Project administration, Writing – review & editing, 1 , 2 Julie Ake, Funding acquisition, Project administration, Resources, Writing – review & editing, 1 Shelly J. Krebs, Writing – review & editing, 1 Sheila Peel, Funding acquisition, Resources, Writing – review & editing, 7 Sodsai Tovanabutra, Project administration, Resources, Supervision, Writing – review & editing, 1 , 2 and Victoria R. Polonis, Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – review & editingcorresponding author 1 ,*
Penny L. Moore, Editor

Associated Data

Supplementary Materials
Data Availability Statement

Abstract

Subtype B HIV-1 has been the primary driver of the HIV-1 epidemic in the United States (U.S.) for over forty years and is also a prominent subtype in the Americas, Europe, Australia, the Middle East and North Africa. In this study, the neutralization profiles of contemporary subtype B Envs from the U.S. were assessed to characterize changes in neutralization sensitivities over time. We generated a panel of 30 contemporary pseudoviruses (PSVs) and demonstrated continued diversification of subtype B Env from the 1980s up to 2018. Neutralization sensitivities of the contemporary subtype B PSVs were characterized using 31 neutralizing antibodies (NAbs) and were compared with strains from earlier in the HIV-1 pandemic. A significant reduction in Env neutralization sensitivity was observed for 27 out of 31 NAbs for the contemporary as compared to earlier-decade subtype B PSVs. A decline in neutralization sensitivity was observed across all Env domains; the NAbs that were most potent early in the pandemic suffered the greatest decline in potency over time. A meta-analysis demonstrated this trend across multiple subtypes. As HIV-1 Env diversification continues, changes in Env antigenicity and neutralization sensitivity should continue to be evaluated to inform the development of improved vaccine and antibody products to prevent and treat HIV-1.

Author summary

HIV-1 has evolved and diversified over the last forty years of the pandemic. These changes have shifted viral properties in such a way that the effectiveness of current therapeutic HIV-specific broadly neutralizing antibodies may be reduced over time. Understanding these changes will allow for the development of improved vaccine and antibody products. In this study, we generated a panel of 30 HIV-1 strains from the United States from samples collected between 2017 and 2018. We characterized the sensitivity of these contemporary HIV-1 strains using 31 antibodies that can inhibit, or neutralize, infection and compared the results to those of viruses from earlier in the pandemic. Our results show that antibodies that were effective against HIV-1 strains from earlier in the pandemic are not as effective against current HIV-1 strains. This trend in reduced sensitivity over time is occurring globally, as similar results were observed in Thailand and Africa. These results will inform the development of more effective antibodies and vaccines to elicit them.

Introduction

The HIV-1 pandemic has been ongoing since 1981 and has resulted in an estimated 84.2 million infections [1]. The global spread of HIV-1 occurred as distinct viral subtypes expanded into different geographical regions [2]. Group M HIV-1 is divided into nine subtypes, including A-D, F-H, J and K, multiple circulating recombinant forms (CRFs), including the prominent CRF01_AE and CRF02_AG, as well as unique recombinant forms [3]. Much of the subtype distribution of group M HIV-1 is preserved today, however diversification has continued over time [47].

HIV-1 subtype B is the primary driver of the HIV-1 epidemics in North and South America, Europe, Australia, the Middle East and North Africa [2]. While HIV-1 incidence rates in the United States (U.S.) have declined by as much as 73% between 1984 and 2019, an estimated 34,800 people acquire HIV-1 yearly and reduced incidence is not consistent across all populations [8]. Over the last four decades, considerable efforts have been made to develop biologic interventions to reduce HIV-1 transmission and disease progression. To date, HIV-1 vaccine candidates have not proven efficacious and have failed to elicit neutralizing antibodies effective against circulating, tier 2 HIV-1 envelope (Env) strains [9]. In the face of these challenges, passive immunization strategies for the prevention and treatment of HIV-1 are also being pursued [10,11].

Broadly neutralizing antibodies (bNAbs) can prevent HIV-1 acquisition through several mechanisms [12,13] and have been shown to be protective in vivo [1417]. While plasma antibodies with the breadth and potency of bNAbs have been shown to develop naturally, they only occur in 20–30% of individuals after 2–3 years of living with HIV-1 [1822]. Prophylactic and therapeutic use of synthetic antibodies may therefore prove effective in the absence of protective active immunity. However, the successful use of these bNAbs may be limited by their potency and breadth against circulating HIV-1 strains, as demonstrated by recent Antibody Mediated Prevention (AMP) trials.

Two AMP trials have evaluated the efficacy of bNAb prophylaxis on HIV-1 acquisition. HVTN 703/HPTN 081 and HVTN704/HPTN 085, were designed to evaluate the protective effect of passive immunization of the HIV-1 bNAb VRC01, which targets the Env CD4 binding site (CD4bs). VRC01 bNAb administration was found to be 75% effective at preventing infection by HIV-1 strains that were susceptible to neutralization by VRC01, with an in vitro sensitivity defined as IC80 <1 μg/ml [23]. These studies demonstrated that the in vitro neutralization sensitivity of the HIV-1 Env from the bNAb recipient dictated the preventative efficacy of the bNAb [23]. These results suggest that bNAbs can be effective at preventing infection when used against circulating, neutralization-sensitive HIV-1 strains and establish baseline potency for potential monoclonal antibodies (mAbs) for HIV-1 prophylaxis.

Understanding the neutralization sensitivity of circulating HIV-1 is therefore critical for prophylactic mAb and vaccine development. Since the beginning of the HIV-1 pandemic, studies have demonstrated a diversification of HIV-1 subtypes and reduction of the neutralization sensitivity of HIV-1 Env [47,24,25]. In this study, we sought to create a contemporary (2017–2018) panel of subtype B Envs from the U.S. and to characterize their neutralization sensitivity to mAbs. We compared neutralization sensitivity of contemporary Envs to Envs from earlier decades and demonstrated significant reductions in subtype B Env neutralization sensitivity over time. Additionally, we demonstrate that the most significant reduction in neutralization sensitivity occurred at target epitopes that were most sensitive to neutralization at the beginning of the pandemic. Results from this study could inform development of improved antibody and vaccine products to prevent and treat HIV-1 and the data strongly support the development of next-generation bNAbs for improved efficacy.

Results

Development of a contemporary HIV-1 subtype B Env panel

Contemporary subtype B HIV-1 envelope genes (envs) were cloned from 30 post residual HIV-1 test plasma collected between 2017 and 2018 from the routine universal screening of U.S. military personnel as part of the RV551 protocol designed to surveil the molecular virology and immunobiology of HIV-1 incidence and prevalence. Twelve individuals were at early stages of infection (Fiebig I-IV) and 19 individuals were at intermediate or late stages of infection (Fiebig V-VI). We sequenced 371 envs from the 30 individuals and evaluated sequence diversity within and between study samples (Fig 1A). In general, individuals in Fiebig stage V or VI showed higher Env diversity than participants in Fiebig stages I-IV (Fig 1A). However, higher Env diversity was observed for two individuals with early-stage (Fiebig II) infections. Sequences from one individual, designated as B027, demonstrated the presence of multiple founder variants and another individual, B016, was determined to have a dual infection. HIV-1 Env amino acid sequences from this contemporary U.S. subtype B panel were then aligned with subtype B Env sequences available in the Los Alamos National Laboratory (LANL) HIV database that were isolated between 2011 and 2020 in the U.S. (Fig 1B). Subtype B Envs from this panel intersperse with subtype B Envs from the past decade (represented as gray lines only). These results suggest the subtype B Envs included in this panel encompass the diversity of subtype B HIV-1 Envs observed in the U.S.

An external file that holds a picture, illustration, etc.
Object name is ppat.1011780.g001.jpg
Phylogenetic relationship of contemporary U.S. subtype B Envs.

(A) Contemporary subtype B HIV-1 Envs were cloned from post residual HIV test samples from 30 U.S. military personnel and sequence diversity was evaluated. Fiebig stage of Env isolation is shown in blue (Fiebig I-IV) and red (Fiebig V-VI); high Env diversity is marked with an asterisk. (B) The phylogenetic relationship was evaluated between the contemporary U.S. subtype B Envs from this study (colored symbols) and subtype B Env sequences collected between 2011–2020 in the U.S. (gray lines) downloaded from the Los Alamos National Laboratory (LANL) HIV sequence database. (C) Divergence from the subtype B Env amino acid most recent common ancestor (MRCA) sequence and (D) pairwise diversity was determined for the consensus Envs from this study (n = 32, magenta), U.S. Envs from 2017–2018 (n = 26, pink), 2011–2016 (n = 239, light green), and prior to 2011 (n = 823, dark green) in the LANL database. Significant differences were determined between consensus Envs from this study and the LANL subset by a Mann Whitney U test. (E) Shannon entropy was estimated at each HXB2 site for the different env alignments, shown as matched colored lines. Pink highlights represent the Env variable loops V1 –V5, from left to right.

Increased sequence diversity of contemporary subtype B Env

HIV-1 Env majority consensus sequences were generated for 30 individuals. Consensus sequences were constructed separately for each of individual B016’s dual infection viruses. We calculated the divergence from the subtype B most recent common ancestor (MRCA) and the pairwise diversity for the Envs from this panel (2017–18, n = 32) and for three sets of publicly available Env sequences collected in the U.S. prior to 2011 (n = 823), between 2011–2020 (n = 266) and in 2017–2018 (n = 26). While Envs from the RV551 panel were significantly more divergent from the subtype B Env MRCA (Fig 1C) and had greater pairwise diversity (Fig 1D) than the subsets of LANL sequences, these sequences were within the expected range of divergence and diversity, i.e. the divergence and diversity values found for RV551 were nested in the distribution corresponding to the larger datasets of circulating sequences (smaller subsets of sequences taken from a set of sequences can have varying (lower or larger) point estimates).

Shannon entropy was estimated at each HXB2 site for the different env alignments. The average entropy across sites was 0.63 (hypervariable regions = 2.20) for Envs from the RV551 panel compared to 0.66 (2.40) for LANL sequences collected prior to 2011, 0.67 (2.39) between 2011–2020, and 0.59 (2.16) between 2017–2018 (Fig 1E). The variability seen among RV551 sequences corresponded to HIV-1 Env variability patterns.

Neutralization sensitivity of subtype B HIV-1 Env pseudoviruses

Contemporary U.S. subtype B Envs from the RV551 protocol were then cloned into expression vectors. Thirty-one HIV-1 pseudovirus (PSV) stocks were produced, including one PSV from each individual; the infectious PSV with Env that was closest to the individual’s consensus sequence was then studied. PSV coreceptor usage was determined; 1 PSV utilized both the CXCR4 and CCR5 coreceptors, while the remaining 30 PSVs utilized only the CCR5 coreceptor. The CXCR4 utilizing PSV was derived from an individual at Fiebig stage IV and was excluded from further analysis.

Neutralization sensitivity of each contemporary U.S. subtype B PSV was assessed using 31 neutralizing antibodies (NAbs) or reagents targeting the Env V1V2, V3, CD4 binding site (CD4bs) and membrane proximal external region (MPER) domains. To demonstrate overall sensitivity of contemporary subtype B PSVs to each NAb, we calculated geometric mean IC50 and IC80 scores and the percentage of viruses with >50% or >80% neutralization observed for each NAb (Fig 2A and 2B). Greater than 50% neutralization was observed for 29 NAbs, while greater than 80% neutralization was observed for only 21 NAbs. NAb IC80 values <1 μg/ml were observed for only three NAbs, V3-specific 10–1074 and PGT128 NAbs and 10E8-iMab (ibalizumab; Fig 2A and 2B). 10E8-iMab is a bispecific antibody that was engineered to target the HIV-1 Env MPER at the 10E8 contact site and the host cellular CD4 molecule [10,26]. 10E8-iMab potency against the contemporary subtype B PSVs ranged from IC80 of 0.0052 to 0.0834 μg/ml, with a geometric mean IC80 of 0.0159 μg/ml; 100% of the PSVs tested were neutralized by 10E8-iMab (Fig 2B).

An external file that holds a picture, illustration, etc.
Object name is ppat.1011780.g002.jpg
Neutralization sensitivities of contemporary subtype B HIV-1 Envs.

The neutralization sensitivity of the contemporary, CCR5-utilizing, subtype B pseudoviruses (PSVs) was determined using 31 neutralizing antibodies (NAbs) or reagents. NAbs are colored by specificity, including: V1V2 loop (orange), V3 loop (blue), CD4 binding site (CD4bs, green), membrane proximal external region (MPER, purple), glycan (CHO, yellow) and bridging region (BR, gray). Geometric mean (GM) (A) IC50 and (B) IC80 scores were calculated for each NAb against all PSVs and compared with the percent of PSVs neutralized by each NAb. The red dotted line at 50 indicates the negative cutoff value for this assay; the gray dotted line at 1 indicates the estimated NAb IC80 value for treatment efficacy, as determined by the AMP studies. Neutralization potency of clinically relevant antibodies and antibody combinations was evaluated against early-stage, subtype B HIV-1 Envs. (C) NAbs were tested individually (shown as circles), or (D) in combination (shown as triangles; open if no bispecific antibody was present and gray if bispecific antibody was present). Combinations of NAbs included three NAbs that targeted different domains, as indicated by the x-axis title and the colored legend.

We then combined three NAbs of different specificities to improve potency over the use of a single NAb. NAbs currently under clinical development [10,2732], including V1V2 NAbs PG9 and PGDM1400, V3 NAbs PGT121 and 10–1074, CD4bs NAbs VRC07 and 3BNC117 and MPER NAb 10E8-iMab, were used to evaluate the neutralization sensitivity of early-stage contemporary subtype B PSVs against all antibody cocktail combinations at the same total starting concentration as used for testing single mAbs. The range of IC80s is shown for each NAb (Fig 2C). The NAb combinations with consistent neutralization potencies of IC80 <0.1ug/mL contained 10E8-iMab. When 10E8-iMab was used at 1/3 of the dose in the mAb cocktail, the potency slightly decreased compared to using 10E8-iMab alone.

PSV neutralization sensitivity was also compared for the contemporary Envs isolated from early-stage infection and those isolated from late-stage infection (S1 Fig). No significant difference between the geometric mean IC50 scores was observed between PSVs from early versus late-stage infection.

Reduced neutralization sensitivity of subtype B HIV-1 Env over time

The neutralization sensitivities of all contemporary subtype B PSVs were compared to historic subtype B strains from prior decades. For this analysis, we utilized neutralization data for historic PSVs available in the CATNAP (Compile, Analyze and Tally Nab Panels) Server on the Los Alamos National Laboratory (LANL) HIV database. For all meta-data analysis, we utilized data from all HIV strains that matched our decade, location, and subtype selection criteria; HIV strains were removed if they were known to utilize the CXCR4 coreceptor. Limited information was available in the LANL HIV database regarding the infection stage of Env isolation or PSV neutralization sensitivity. The number of HIV strains included in all meta-data analyses are shown in S1 Table. To ensure concordance between data generated in our MHRP laboratory and previously published CATNAP data, we performed a bridging analysis using 3 PSVs and 18 NAbs and observed a highly significant correlation between IC50 endpoints (S2 Fig, p<0.0001, r = 0.9043).

Contemporary subtype B PSVs had significantly higher IC50 values than viruses from the 1980s (22/31 NAbs), 1990s (16/31 NAbs), 2000s (10/31 NAbs), as shown in Fig 3. The CATNAP HIV database had limited neutralization data available for other viruses from the 2010s. Significant differences between IC50 values of subtype B PSVs from different decade periods are indicated. Overall, contemporary subtype B PSVs were less sensitive to NAb neutralization than historic subtype B PSVs.

An external file that holds a picture, illustration, etc.
Object name is ppat.1011780.g003.jpg
Changes in neutralization sensitivities of subtype B HIV-1 Envs over time.

Neutralization sensitivities of contemporary U.S. subtype B pseudoviruses (PSVs) (2017–18) were compared with U.S. subtype B PSVs from earlier decades. Each graph demonstrates the neutralization profile of one neutralizing antibody (NAb), as indicated in the title above each graph. Symbols are colored by the NAb specificity, with V1V2-specific NAbs in orange, V3-specific NAbs in blue, CD4bs NAbs in green, MPER NAbs in purple, glycan-specific (CHO) NAbs in yellow, and bridging region-specific NAbs in gray. Significant differences were determined by Mann Whitney U test; * = p<0.01, ** = p<0.001, *** = p<0.0001.

We then evaluated Env domain specific changes in subtype B PSV sensitivity over time. For each NAb, we created a geometric mean score of IC50 values for all PSVs from a given decade period (Fig 4A–4D). Similarly, we calculated the percentage of PSVs from a given decade period that were neutralized by each NAb (Fig 4F–4I). For this analysis, we included the contemporary panel of 2017–18 PSVs with other PSVs from the 2010s, when available. NAbs were then grouped by target Env domain and the domain trends over time are shown. We observed significant and continuously reduced neutralization sensitivity to V3 glycan NAbs over time (Fig 4B and 4G); similar trends were observed for the CD4bs and MPER targeting NAbs (Fig 4C, 4D, 4H, and 4I). While slight variations over time were observed for the V1V2 NAbs, the trend was not significantly directional (Fig 4A and 4F). Significant differences over time were determined by rank comparison, as shown, as well as by trend analysis for ordered differences. NAbs targeting the V3, CD4bs, and MPER demonstrated a significant trend in increasing order of median IC50s (p = 0.0001, 0.0135, and 0.0249, respectively) and decreasing order of median breadth (p = 0.0001, 0.0598, and 0.0098, respectively) over time.

An external file that holds a picture, illustration, etc.
Object name is ppat.1011780.g004.jpg
Changes in domain-specific neutralization sensitivities of subtype B HIV-1 Envs over time.

For each neutralizing antibody (NAb), the geometric mean (GM) IC50 and percentage of pseudoviruses (PSVs) neutralized, was determined for all PSVs from each decade period, as shown in Fig 3. The trends for NAb potency (A-D) and breath (F-I) over time are shown by NAb specificity, including (A and F) V1V2 specific NAbs, (B and G) V3 NAbs, (C and H) CD4bs NAbs and (D and I) MPER NAbs. This neutralization analysis was repeated using 33 subtype B HIV+ plasma samples against PSVs from each decade period; the (E) GM ID50 and (J) percentage of PSVs neutralized, was determined for all PSVs from each decade. Significant differences were determined by Wilcoxon matched-pairs signed rank test; * = p<0.01, ** = p<0.001, *** = p<0.0001.

We evaluated this trend of reduced Env neutralization sensitivity over time using polyclonal neutralizing reagents. We tested the neutralization potency of 33 plasma samples from individuals with subtype B HIV-1 against PSVs from the different decade periods, including 5 PSVs from the 1980s (SF162 was excluded from mean, median and statistical calculations), 8 PSVs from the 1990s, 6 PSVs from the 2000s, and 22 PSVs from the 2017–18 contemporary panel, including 11 from Fiebig I-IV and 11 from Fiebig V-VI (S3 Fig). The 11 PSVs with the highest titer were chosen for each group. We created a geometric mean ID50 score and determined the number of PSVs neutralized for each plasma across each PSV decade period (Fig 4E and 4J, respectively). Similar to the trends observed for the mAbs, we observed significantly decreasing PSV sensitivity (as indicated by decreasing geometric mean [GM] ID50, S3A Fig) to the polyclonal plasma and a decreasing order of median IC50s (p<0.0001) and breadth (p<0.0001, S3B Fig)) over time. Additionally, we stratified the HIV+ plasma samples by decade period of sample collection and observed a reduction in plasma potency and breadth over time (S3A and S3B Fig). A heat map showing all of the individual plasma reactivities is shown in S3C Fig; the SF162 data points were not included in calculations, as this represents an outlier as an exceptionally sensitive tier 1 virus.

NAbs with the greatest potency against 1980s PSVs appeared to have the greatest rate of neutralization resistance over time. We therefore evaluated the fold change in NAb potency against the PSVs from the 1980s and 2010s in comparison with the initial NAb potency against 1980s PSVs (Fig 5A). A positive fold change in Env sensitivity to a given NAb indicates increasing sensitivity over time (1980s PSV IC50 > 2010s PSV IC50), while a negative fold change in Env sensitivity indicates decreasing sensitivity over time (1980s PSV IC50 < 2010s PSV IC50). The dashed vertical line on the x-axis at 1 indicates no change was observed over time. We identified a significant direct correlation between changes in U.S. subtype B Env neutralization sensitivity over time and initial 1980s Env sensitivity (Fig 5A; p = 0.0059, r = 0.4990). This indicates that NAbs with the greatest initial neutralization potency experienced the greatest longitudinal reduction in neutralization potency (increased resistance).

An external file that holds a picture, illustration, etc.
Object name is ppat.1011780.g005.jpg
Associations between initial neutralizing antibody potency and rate of declining neutralization potency over time.

For each neutralizing antibody (NAb), the geometric mean (GM) IC50 was determined for all PSVs from an earlier and later decade period, as shown in the x-axis title. Fold change in NAb potency was determined between the earlier and later decade periods. A positive fold change in Env sensitivity to a given NAb indicates increasing sensitivity (earlier PSV GM IC50 > later PSV GM IC50), while a negative fold indicates decreasing sensitivity over time (earlier PSV GM IC50 < later PSV GM IC50). The vertical dashed line at 1 indicates no observed change. The correlation between the fold change in NAb potency over time and the potency of the NAb against PSVs from the earlier decade period are shown for (A) subtype B from the United States, (B) CRF01_AE from Thailand, (C) subtype A2 from East Africa, (D) subtype D from Uganda, and (E) subtype C from South Africa. Individual NAbs are shown as unique symbols colored by target specificity, including NAbs targeting the V1V2 (orange), V3 (blue), CD4bs (green) and MPER (purple) domains. Correlations were evaluated by Spearman correlation analysis. The summary of fold change in NAb potency over a ten-year period is shown for (F) all NAbs from all regions; data are stratified by initial NAb potency and significant differences were determined by Mann Whitney U test; ** = p<0.001, *** = p<0.0001.

To determine whether this trend was applicable to other subtypes, we evaluated neutralization data from the LANL CATNAP database and our previous work [33], to characterize four additional subtypes (Fig 5B–5E). Due to the limits in data availability, not all NAbs are demonstrated and the fold change in Env sensitivity was determined over a two- or three-decade period, as indicated in each chart. Similar trends in the association between initial NAb potency and rate of potency decline over time were observed for CRF01_AE from Thailand (Fig 5B; p = 0.0063, r = 0.7055), subtype A1 from East Africa (Kenya, Rwanda, Tanzania and Uganda; Fig 5C; p<0.0001, r = 0.8087), subtype D from Uganda (Fig 5D; p = 0.0087, r = 0.5571) and subtype C from South Africa (Fig 5E; p<0.0001, r = 0.7545). Changes in NAb potency over a ten-year period (1990s to 2000s) were determined for all NAbs and regions; the results are summarized in Fig 5F. NAbs with an initial potency of <0.5 μg/ml IC50 demonstrated a significant reduction in potency, whereas NAbs with an initial potency of >5 μg/ml IC50 demonstrated a significant increase in potency over time, indicating that changes in both sensitivity and resistance were observed. Minimal change in potency was observed for NAbs with IC50s between 0.5–5 μg/ml. Of note, the Ugandan subtype D viruses showed a less frequent decrease in sensitivity, as compared to the other 4 subtypes.

Subtype B Env features associated with reduced neutralization sensitivity

Lastly, we evaluated changes in Env features that may impact the observed reduced neutralization sensitivity. For this analysis, we used Env sequences for the PSVs used to compare neutralization sensitivity over time, shown in Fig 3. Utilizing the subtype B HXB2 strain as the reference sequence, we calculated the number of mismatched amino acids in the contact sites for four NAbs, including PGDM1400, 10–1074, VRC01 and 10E8. Envs were grouped by decade to demonstrate the accumulation of contact site mismatched amino acids over time. The number of mismatched amino acids significantly increased over time for PGDM1400 and VRC01 (Fig 6A and 6C); a trend towards increasing number of mismatched amino acids was observed for 10–1074 (Fig 6B) and no trend was observed for 10E8 (Fig 6D).

An external file that holds a picture, illustration, etc.
Object name is ppat.1011780.g006.jpg
Changes in subtype B Env neutralizing antibody contact site deviations over time and impact on Env neutralization sensitivity.

Subtype B Env sequences for pseudoviruses included in the neutralization analysis for each decade period were analyzed to determine the number of mismatched amino acids from the HXB2 reference strain that occurred in contact sites for neutralizing antibodies (NAbs) (A) PGDM1400, (B) 10–1074, (C) VRC01, and (D) 10E8. The relationship between the number of mismatched residues in the NAb contact site and Env NAb neutralization sensitivity were determined for (E) PGDM1400, (F) 10–1074, (G) VRC01, and (H) 10E8. Significant differences were determined by Mann Whitney U test; * = p<0.01, ** = p<0.001, *** = p<0.0001. Correlations were evaluated by Spearman correlation.

We then evaluated the relationship between the number of mismatched residues in the NAb contact site and Env neutralization sensitivity to that NAb. While the data were widely distributed, we observed a significant direct correlation between the number of mismatched residues and reduced neutralization potency for NAbs PGDM1400 (Fig 6E; p = 0.0032, r = 0.2409) and 10–1074 (Fig 6F; p = 0.0012, r = 0.2626); no significant associations were observed for NAbs VRC01 or 10E8.

To visualize potential changes in NAb contact site amino acids over time, we generated logo plots for each NAb using Env sequences downloaded from the LANL sequence database for U.S. subtype B Envs from the 1980s (N = 75) and 2010s (N = 785) (S4 Fig). The probability of each amino acid, indicated by one-letter symbols, at each position is indicated by the height of the letter. A gray bar indicates a gap or missing amino acid compared to the alignment, at that position. For PGDM1400, we can see a loss of Q at position 170. For 10–1074, we can see changes in potential N-linked glycosylation sites (PNGSs), including the N332 glycan. From our neutralization analysis, all PSVs from the 1980s had PNGS N332, while 15% of PSVs from our contemporary panel lacked the PNGS N332, making them resistant to NAb 10–1074 (S5 Fig). For VRC01, we can see variations across the contact sites, while the contact sites for 10E8 remain mostly unchanged. These results suggest that the accumulation of mismatched residues in individual NAb contact sites only partially contributes to the reduced neutralization sensitivity we observed for certain NAbs.

Discussion

The diversity and mutability of HIV-1 remain as significant challenges for the development of effective vaccines and other products to prevent and treat HIV-1 [34]. Given the diversity in HIV-1 Env between subtypes, a median of 25% difference at the amino acid level [35], subtype-specific vaccines that are matched to regional circulating strains have been pursued. HIV-1 Envs from each subtype have similar antigenic features and neutralization-sensitive targets [3638]. However, even within subtypes, a median Env amino acid difference of 15% exists [35]. Our study, along with others [47], demonstrates a continuing diversification of HIV-1 Env within subtype over time. Population level adaptation of HIV-1 has been shown to associate with reduced neutralization sensitivity and Env immunogenicity [47,24,25]. These adaptations may be impacted by humoral immune responses to natural infection, as well as partially effective vaccines and other products [39].

In this study, we focused on subtype B HIV-1 Env, as subtype B has been the primary driver of the HIV-1 epidemic in the United States (U.S.) for over forty years and is also a prominent subtype in North and South America, Europe, Australia, Middle East and North Africa. Bouvin-Pley et al. previously demonstrated reduced neutralization sensitivity of subtype B PSVs between 1987 and 2010 [6]. Significant decreases in sensitivity were observed for 8 out of the 13 NAbs tested, however no differences were observed in gp120 of V3 length or in number of PNGS [6]. In a further study, antigenic changes associated with reduced neutralization sensitivity also correlated with increased infectivity and reduced sensitivity to fusion inhibitors and CCR5 antagonists [40]. Similarly, increased diversification and reduced neutralization sensitivity have been reported for other subtypes. Rademeyer et al. showed that the subtype C epidemic in South Africa diversified significantly between 1998 and 2010 [25]. Diversification was not associated with changes in variable loop length or glycan density, however, significant decreases in neutralization sensitivity were observed for NAbs 4E10, VRC01 and PG9 [25]. Stefic et al. demonstrated a similar trend for the CRF02_AG epidemic in France between 1997 and 2012 [7]. For CRF02_AG, significant direct correlations between amino acid diversification and reduced neutralization sensitivity were observed for CD4bs NAbs 3BNC117 and NIH45-46.

Our analysis of U.S. subtype B HIV-1 encompassed a forty-year period and employed a larger number of NAbs that allowed us to identify trends associated with reduced Env neutralization sensitivity. We observed an association between the number of mismatched residues in the NAb contact site and reduced neutralization sensitivity. In some cases, this reduced sensitivity corresponded to loss of PNGS or other NAb contact site deviations, but not consistently. Amino acid divergence may also be associated with tertiary shifts in HIV-1 Env structure that could mitigate neutralization at the most sensitive Env targets. However, this study was not designed to evaluate this in depth.

Similar to the generalized amino acid divergence observed across Env, we observed reduced neutralization sensitivity of contemporary subtype B HIV-1 across Env domains over time. The degree of reduced neutralization sensitivity that we observed was associated with the potency of the specific NAb against early-pandemic HIV-1 strains. This supports the importance of bNAb potency in effective in vivo viral inhibition and in effecting population level changes in Env. Interestingly, while we observed a significant reduction in potency for NAbs with IC50s <0.5 μg/ml, we also observed a significant increase in potency for NAbs with IC50s >5 μg/ml. Env mutations that can circumvent potent NAbs may allow for changes in Env neutralization sensitivity as distal sites. As Env neutralization sensitivity can vary by HIV-1 subtype, patterns of reduced neutralization sensitivity will develop and vary accordingly. Most strikingly, subtype B and CRF01_AE HIV-1 have distinct neutralization profiles; while subtype B Envs are most sensitive to V3-specific NAbs, CRF01_AE Envs are most sensitive to V1V2-specific NAbs. Indeed, we observed the greatest decrease in NAb potency for V3-specific NAbs for subtype B and V1V2-specific NAbs for CRF01_AE. This trend of declining NAb potency was further demonstrated for three additional subtypes, using data from the LANL CATNAP database, as well as data from our laboratory. However, infection stage or neutralization sensitivity (tier) of the virus strains used in this meta-analysis were not factored into the data selection criteria, as the information was not available for many strains.

The recent AMP studies demonstrated that bNAbs could effectively prevent HIV-1 acquisition if the HIV-1 strain to which the individual is exposed is potently neutralized by the bNAb in vitro (IC80 [concentration of antibody that neutralizes 80% of the virus] <1 μg/ml). These results provide a baseline for effective prophylactic bNAb potency that was previously not defined [23]. The neutralization data generated in this study were produced using the same in vitro neutralization assay platform used to assess the AMP participant HIV-1 strains. Our results show that NAb potency of <1 μg/ml IC80 may be difficult to achieve against contemporary U.S. subtype B viruses. However, our NAb selection was based on availability at initiation of the study, and more potent reagents are currently in development or available. Only 5 of the 31 NAbs evaluated in this study achieved potency <1 μg/ml IC80 against more than 50% of the contemporary U.S. subtype B strains; only 3 NAbs (V3 specific PGT128 and 10–1074, and 10E8-iMAb) achieved a GM potency <1 μg/ml IC80 overall. 10E8-iMab potently neutralized the contemporary subtype B PSVs with a GM IC80 of 0.0139 μg/ml and 100% breadth. Wagh et al. previously demonstrated that the potency of 10E8-iMab exceeded <1 μg/ml IC80 against panels of subtype A, C, and D PSVs, with 96–100% neutralization coverage across the subtypes [41]. This antibody is now advancing in clinical trials; these types of antibodies are promising as they achieve the cross-subtype potency needed for globally effective products.

While reduced neutralization sensitivity over time was observed for HIV-1 Env, the changes reported in this study occurred over decades. We observed development of NAb resistance over time, as defined by the limits of our neutralization assay, but the contemporary Envs were not completely neutralization resistant. Additionally, this study only characterized HIV-1 Envs from a small convenience sample set from 30 individual infections. Changes in Env antigenicity may impact accessibility of contact sites for known bNAbs. However, these changes in antigenicity may also open opportunities for next generation bNAb targets on Env and development of novel bNAbs from contemporary infections should continue to be advanced. While changes in antigenicity over time may have reduced Env neutralization sensitivity, it is unknown how these changes have impacted other antibody functions, such as antibody-dependent cellular cytotoxicity (ADCC) and antibody-dependent cellular phagocytosis (ADCP), which are important parameters for additional study.

This work further highlights the need for continued surveillance of circulating HIV-1 strains and the inclusion of contemporary HIV-1 in vaccine analysis viral panels. As results from the AMP studies highlighted, the effectiveness of a given antibody-based treatment or vaccine is dependent on the sensitivity of the circulating Envs to the synthetic or elicited antibodies. Development and co-administration of next generation antibody-based prophylaxis may be required to achieve the level of breadth and potency needed to provide protective coverage [41]. Continued detailed characterization of circulating HIV-1 Envs should inform improved targeted development of immune-based interventions for the prevention and treatment of HIV-1.

Materials and methods

Cohort descriptions and ethics statement

Contemporary subtype B envs were cloned from plasma samples selected from the RV551 study, a non-human subjects research protocol approved by the Walter Reed Army Institute of Research (WRAIR) Human Subjects Protection Branch and the US Medical Research Development Command Office of Human Research Oversight. This study adhered to the policies for protection of human subjects, as prescribed in Army Regulations 70–25: Use of Volunteers as Subjects of Research. Formal consent was obtained in writing. RV551 was designed to assess the molecular virology and immunobiology of incident and prevalent contemporary U.S. Military HIV-1 subtype B infections. Samples were obtained from post residual HIV-1 test specimens submitted to the HIV Diagnostic and Reference Laboratory between 2017 and 2018 as part of universal screening for HIV-1 in U.S. military service members. Samples were anonymized and linked to HIV meta test data by a dis-interested data steward; Fiebig stage was determined by analysis of meta test data.

HIV-1 Env sequencing

Env genes were sequenced via Single Genome Amplification (SGA). HIV-1 RNA was extracted from plasma using the QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany). Complementary DNA (cDNA) was synthesized using Superscript III (Invitrogen Corp., Carlsbad, CA) as instructed by the manufacturer using the primer UNINEF-7’ (5’-GCACTCAAGGCAAGCTTTATTGAGGCTT-3’; HXB2 9632–9605 nt). Using near-endpoint dilution of cDNA template, a first round PCR was performed to amplify a near full length or 3’ half-length genome by using MSF12b/UNINEF-7’ (5’- GAAGCYATGCATGGACAAGTRGA-3’; 623–649 nt) or PolJv2/UNINEF-7’ (5’- GAAGCYATGCATGGACAAGTRGA-3’; HXB2 4371–4392 nt). From the first round PCR product, the env gene was amplified using the primers BH4minus (5’-TAGGCATTTCCTATGGCAGGAAGAAG; HXB2 5958–5983 nt) and BH2NOENZ2ACE (5’-GTCTCGAGATACTGCTCCTACTC; HXB2 8904–8882 nt). Both rounds of PCR were performed with Q5 polymerase (New England Biolabs, Ipswich, MA, USA). Amplicons were gel purified using the QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany) and Sanger sequenced on a 3730xl capillary sequencher (Applied Biosystems Foster City, CA, USA). HIV-1 subtype was determined using the NCBI genotyping tool (https://www.ncbi.nlm.nih.gov/projects/genotyping/formpage.cgi).

Phylogenetic and sequence analysis

For the phylogenetic analysis: RV551 amino acid sequences were aligned to 2011–2020 circulating subtype B Env sequences using the add function in the multiple sequence alignment program, MAFFT v7.475 [42]. The phylogeny was constructed from the alignment of subtype B and RV551 sequences using the IQ-TREE 2 software package for phylogenetic inference using the maximum likelihood (ML) criterion [43], with an HIV-1 between-patient matrix [44]. A sequence corresponding to the most recent common ancestor (MRCA) of subtype B Envs was reconstructed using FastML v.3.11 [PMID: 22661579] with branch-length optimization and a gamma distribution on subtype B Envs in the LANL database collected prior to 2010. The distribution of distances of RV551 sequences to the subtype B MRCA sequence was determined by computing pairwise distances to the consensus of sequences for each RV551 plasma sample assuming an HIV-1 between-sample matrix substitution model; the distribution of distances to the subtype B MRCA was also determined for subtype B Env sequences sampled between 1980–2010, between 2011–2020, in the U.S. between 2017–2018 and for the RV551 sequences. Pairwise distances were calculated across sequences within each subset. Shannon entropy, H = –∑pilog2pi, was calculated at each site for each subset. For the NAb contact site mismatch analysis, Env amino acid sequences for PSVs included in the neutralization analysis for each decade period were analyzed to determine the number of mismatched amino acids in the NAb contact sites. To visualize Env amino acid differences between decade periods, logo plots were generated using U.S. subtype B Envs selected from the LANL HIV sequence database from 1980s (N = 75) and 2010s (N = 785).

Generation of envelope gene clones for PSV production

Envelope glycoprotein (gp160) genes were amplified by single genome amplification and cloned into eukaryotic expression vectors for use as PSVs, as previously described [45,46], Briefly, The amplicon was gel purified using the QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany), cloned into pcDNA3.1/V5-His-TOPO eukaryotic cloning vector (Invitrogen, Waltham, MA, USA), and transformed into either STBL2 or STBL4 cells (Invitrogen, Waltham, MA, USA). Plasmid DNA was purified from bacterial cultures using the Qiagen Plasmid Maxi Kit (Qiagen, Hilden, Germany). Env sequence was confirmed by comparison with original SGA-derived amplicon. There were no amino acid differences between the cloned env genes and the uncloned env amplicon. Env function was confirmed by infection of the TZMbl reporter cell line (ATCC, Manassas, VA, USA).

PSV production and titration

PSVs were prepared as previously described [46]. Briefly, PSVs were produced by co-transfection of HEK293T/17 cells (ATCC, Manassas, VA, USA) with a pSG3ΔEnv DNA plasmid encoding the HIV-1 backbone (24 μg DNA) and a plasmid encoding the HIV-1 env genes (8 μg DNA). Plasmids were combined with X-tremeGENE 9 DNA transfection reagent (Sigma-Aldrich, St. Louis, MO, USA) at a 1:3 ratio in serum-free DMEM, as per manufacturer recommendations. The transfection mix was added to a T75 flask containing plated HEK293T/17 cells in culture medium. The flask was incubated for 6 hours at 37°C. Then the media was replaced. PSVs were harvested 72 hours after transfection, filtered through a 0.45 μm PVDF membrane (MilliporeSigma, Burlington, MA, USA) and stored in 3 ml aliquots at -80°C. PSVs were then titrated on TZMbl cells to determine an acceptable dilution for neutralization assays. PSVs were serially diluted, then added to 384-well black-walled plates (Thermo Fisher Scientific, Waltham, MA, USA) in a final 12.5 μL/well volume. An additional 12.5 μL of supplemented DMEM was added to each well. TZM-bl cells were added to plates in 25 μL of supplemented DMEM with 40 μg/mL of DEAE-Dextran (Sigma, St. Louis, MO, USA). PSVs and cells were incubated together for 48 hours at 37°C. The supernatant was then removed from each well to leave a total of 25 μL of supernatant per well. Infectivity was measured by luminescence after adding 25 μL of the Bright-Glo Luciferase Assay System (Promega, Madison, WI, USA) to each well; luminescence was quantified using the EnVision luminometer (Perkin Elmer, Waltham, MA, USA).

Coreceptor determination by the GHOST cell assay

The GHOST cell infection assay was used to determine PSV coreceptor usage. Parental, CXCR4-expressing, or CCR5-expressing GHOST cells were cultured in 96-well plates at 1.6x104 cells/well and incubated at 37°C overnight. Cells were infected with undiluted PSV in the presence of 20 μg/ml of polybrene infection reagent (MilliporeSigma, Billerica, MA, USA) for 4 hours. Cells were washed and cultured for 2 days, then harvested and fixed with 4% paraformaldehyde. The percentage of infected cells expressing green fluorescent protein was measured by flow cytometry analysis using an LSRII cytometer and FACSDIVA software (Becton-Dickinson, San Jose, CA, USA). Post-acquisition analysis was conducted with FlowJo software (FlowJo, LLC, Ashland, OR, USA). PSVs were designated as using the CXCR4 or CCR5 receptor if the ratio to control negative (RTCN) was >10; RTCN = (MFI x %pos)infected / (MFI x %pos)uninfected. Assay controls included an uninfected negative control, murine leukemia virus as a positive control, and known CCR5 and CXCR4 utilizing PSVs as cell line controls.

Antibodies and plasma

The following reagents were obtained through the HIV Reagent Program, Division of AIDS, NIAID, NIH (Manassas, VA, USA): 10–1074, 10E8, 2191, 2219, 3074, 35O22, 3869, 3BNC117, 7H6, CH02, CH03, CH04, HJ16, NIH45-64, PGT121, PGT126, PGT128, PGT145, sCD4, VRC01, and VRC03. The following reagents were purchased from Polymun (Klosterneuburg, Austria): 2F5, 2G12, 447-52D, 4E10, b12, PG16, PG9. VRC07 and 10E8-iMab were graciously provided by Yaoxing Huang and David Ho from the Aaron Diamond AIDS Research Center, Columbia University Medical Center. PGDM1400 and PGT130 were graciously provided by Dr. Devin Sok at the International AIDS Vaccine Initiative. 830A was graciously provided by Dr. Susan Zolla-Pazner at the Icahn School of Medicine at Mount Sinai, New York, NY.

Additional samples for longitudinal virus characterization were also obtained from the U.S. Military HIV Natural History Study (NHS), a 37-year prospective, continuous enrollment cohort study of active duty service members and Military Health System beneficiaries infected with HIV [47]. The NHS protocol has been IRB-approved centrally at the Uniformed Services University of the Health Sciences, with review at participating sites and WRAIR. All NHS subjects provide written informed consent. All selected plasma samples were from chronic, subtype B HIV infections.

PSV neutralization

Neutralizing antibody titers were determined using TZM-bl cells [48] in a high-throughput assay utilizing robotic liquid handling, as previously described [33]. Plasma samples were prepared by heat inactivation at 56°C for 45 min, followed by 6 min centrifugation at 3000 rpm to remove particulates. Plasma samples were prepared at a starting dilution of 1:20 in supplemented DMEM; NAbs were prepared at a starting concentration of 25 μg/mL in supplemented DMEM. PSVs were diluted to a concentration to yield 100,000 relative light units (RLUs) and at least 3 times above the cell-only controls, as determined by the PSV titration assay. Neutralization reagents were serially diluted using a Biomek NXP liquid handler (Beckman Coulter, Indianapolis, IN, USA). Titered NAbs and plasma samples were transferred to 384-well culture plates and incubated with an equal volume of PSV for 1 hour at 37°C. TZMbl cells (3x103 cell/well) mixed with 40 μg/mL DEAE-dextran (Sigma, St. Louis, MO, USA) were added to each well and incubated for an additional 48 hours. RLUs were detected with the EnVision luminometer (Perkin Elmer, Waltham, MA, USA) using the Bright-Glo Luciferase Assay System (Promega Corporation, Madison, WI, USA). Percentage neutralization (percentage reduction of relative light units in the presence of plasma) was calculated for each plasma dilution. Neutralization dose–response curves were fitted by nonlinear regression using the LabKey Server, and the final titer is reported as the concentration of NAb (IC50/IC80) or reciprocal of the dilution of plasma (ID50/ID80) necessary to achieve 50% or 80% neutralization.

Meta-data selection

Neutralization data for historic subtype B, and non-subtype B PSVs were selected from the CATNAP Server on the LANL HIV database. For all meta-data analysis, we utilized data from all HIV strains that matched our decade, location, and subtype selection criteria; HIV strains were removed if they were known to utilize the CXCR4 coreceptor. The number of HIV strains included in all meta-data analyses are shown in S1 Table.

Statistical analysis

Significant differences between neutralization titers were determined by Mann Whitney U test or Wilcoxon matched-pairs signed rank, as indicated; * = p<0.01, ** = p<0.001, *** = p<0.0001. The Jonckheere’s trend test was used to determine statistically significant trends over time. Significant correlations between data sets were determined by Spearman correlation; where relevant, trend lines were shown.

Supporting information

S1 Fig

Neutralization sensitivity of early-stage versus late-stage subtype B HIV-1 Envs.

Pseudoviruses were categorized by infection stage to demonstrate differences in neutralization sensitivity between Envs from early-stage infections (Fiebig I-IV, black bars), as compared to late-stage infections (Fiebig V-VI, red bars). The red dotted line at 50 indicates the negative cutoff value for this assay. Significant differences were determined by Mann Whitney U test; * = p<0.01.

(TIFF)

S2 Fig

Correlation between neutralization endpoints.

Neutralization endpoints were compared between data generated in our MHRP (U.S. Military HIV Research Program) laboratory and data available on the CATNAP (Compile, Analyze and Tally Nab Panels) Server on the Los Alamos National Laboratory (LANL) HIV database. Neutralization log10 IC50 values were compared for the same 3 pseudoviruses, including SF162, MN and AC10, against 18 NAbs. The correlation between data sets was determined by Spearman correlation; the trend line is shown.

(TIFF)

S3 Fig

Changes in HIV+ plasma neutralization sensitivities of subtype B HIV-1 Envs over time.

Neutralization potencies of 33 subtype B HIV+ plasma samples collected in the 1990s, 2000s and 2010s were determined against subtype B pseudoviruses (PSVs) with Envs isolated in the 1980s, 1990s, 2000s and 2010s. The A) geometric mean (GM) score of plasma ID50 values and B) breadth (% of neutralized viruses) for each plasma:virus decade period were determined. C) The heat map shows the log10 ID50 value of each plasma sample against each PSV. SF162 PSV, determined to be exceptionally neutralization-sensitive, was excluded from the analysis but the data are included here for reference.

(TIFF)

S4 Fig

Broadly neutralizing antibody contact sites for 1980s and 2010s HIV-1 Envs.

Subtype B HIV-1 Env sequences from the 1980s (N = 75) and 2010s (N = 785) were from the Los Alamos National Laboratory HIV sequence database. Amino acid contact sites of four broadly neutralizing antibodies, including PGDM1400 (top left), 10–1074 (top middle), 10E8 (top right), and VRC01 (bottom), are shown as logo plots for each decade period. The probability of each amino acid, indicated by one-letter symbols, at each position is indicated by the height of the letter. A gray bar indicates a gap, or missing amino acid compared to the alignment, at that position.

(TIFF)

S5 Fig

Impact of N332 glycan on 10–1074 broadly neutralizing antibody neutralization.

Pseudoviruses from the 1980s and 2017–18 periods were evaluated for presence (N332+) or absence (N332-) of asparagine at Env position 332. Env sensitivity to the 10–1074 broadly neutralizing antibody is demonstrated for N332+ and N332- pseudoviruses. NA indicates that no Envs of that genotype were available. Significant differences were determined by Mann Whitney U test; *** = p<0.0001.

(TIFF)

S1 Table

Number of viruses utilized from each decade for mAb neutralization analyses.

The number of pseudoviruses tested from each decade using each NAb are tabulated for the analyses presented in the indicated figures.

(TIFF)

Acknowledgments

We would like to thank our collaborators for their kind contributions of reagents. 10E8-iMab was graciously provided by Yaoxing Huang and David Ho from the Aaron Diamond AIDS Research Center, Columbia University Medical Center. PGDM1400 and PGT130 were graciously provided by Dr. Devin Sok at the International AIDS Vaccine Initiative. 830A was graciously provided by Dr. Susan Zolla-Pazner at the Icahn School of Medicine at Mount Sinai. We also thank the protocol staff and research teams. We would like to thank Xiuping Chu and Katie Lammersen for their assistance with plasma retrieval and Fiebig stage data.

Disclaimer

The views expressed are those of the authors and do not reflect the official views of the Walter Reed Army Institute of Research, Uniformed Services University of the Health Sciences, the Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc., the National Institutes of Health or the Department of Health and Human Services, Brooke Army Medical Center, the U.S. Army Medical Department, the U.S. Army Office of the Surgeon General, the Department of Defense, Defense Health Agency, the Departments of the Army, Navy or Air Force, and U.S. Government. Mention of trade names, commercial products, or organizations does not imply endorsement by the U.S. Government. Some of the authors of this work are employees of the United States Government. This work was prepared as part of their official duties. Title 17 U.S.C. § 105 provides that ‘Copyright protection under this title is not available for any work of the United States Government’. Title 17 U.S.C. § 101 defines a U.S. Government work as a work prepared by a military service member or employee of the U.S. Government as part of that person’s official duties.

Funding Statement

This work was supported by a cooperative agreement (W81XWH-18-2-0040) between the Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc., and the U.S. Department of Defense. Coauthors affiliated with the Henry M. Jackson Foundation and the U.S. Military HIV Research Program received salaries from this funding mechanism. Support for portions of this work (IDCRP-000-37) was provided by the Infectious Disease Clinical Research Program (IDCRP), a Department of Defense program executed through the Uniformed Services University of the Health Sciences, Department of Preventive Medicine and Biostatistics through a cooperative agreement with The Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc. (HJF). These portions have been funded by the National Institute of Allergy and Infectious Diseases, National Institutes of Health, under Inter-Agency Agreement Y1-AI-5072, and the Defense Health Program, U.S. DoD, under award HU0001190002. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

All relevant data are within the manuscript and its Supporting Information files.

References

1. UNAIDS. Global HIV statistics 2022 fact sheet. 2022.
2. Bbosa N, Kaleebu P, Ssemwanga D. HIV subtype diversity worldwide. Curr Opin HIV AIDS. 2019;14(3):153–60. Epub 2019/03/19. 10.1097/COH.0000000000000534 . [Abstract] [CrossRef] [Google Scholar]
3. Sharp PM, Hahn BH. Origins of HIV and the AIDS pandemic. Cold Spring Harb Perspect Med. 2011;1(1):a006841. Epub 2012/01/10. 10.1101/cshperspect.a006841 ; PubMed Central PMCID: PMC3234451. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
4. Bunnik EM, Euler Z, Welkers MR, Boeser-Nunnink BD, Grijsen ML, Prins JM, et al.. Adaptation of HIV-1 envelope gp120 to humoral immunity at a population level. Nat Med. 2010;16(9):995–7. Epub 2010/08/31. 10.1038/nm.2203 . [Abstract] [CrossRef] [Google Scholar]
5. Stefic K, Bouvin-Pley M, Braibant M, Barin F. Impact of HIV-1 Diversity on Its Sensitivity to Neutralization. Vaccines (Basel). 2019;7(3). Epub 2019/07/28. 10.3390/vaccines7030074 ; PubMed Central PMCID: PMC6789624. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
6. Bouvin-Pley M, Morgand M, Moreau A, Jestin P, Simonnet C, Tran L, et al.. Evidence for a continuous drift of the HIV-1 species towards higher resistance to neutralizing antibodies over the course of the epidemic. PLoS Pathog. 2013;9(7):e1003477. Epub 2013/07/16. 10.1371/journal.ppat.1003477 ; PubMed Central PMCID: PMC3701719. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
7. Stefic K, Bouvin-Pley M, Essat A, Visdeloup C, Moreau A, Goujard C, et al.. Sensitivity to Broadly Neutralizing Antibodies of Recently Transmitted HIV-1 Clade CRF02_AG Viruses with a Focus on Evolution over Time. J Virol. 2019;93(2). Epub 2018/11/09. 10.1128/JVI.01492-18 ; PubMed Central PMCID: PMC6321924. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
8. CDC. Estimated HIV incidence and prevalence in the United States, 2015–2019. HIV Surveillance Supplemental Report 2021. CDC reports. 2021;26(1). [Google Scholar]
9. Kim J, Vasan S, Kim JH, Ake JA. Current approaches to HIV vaccine development: a narrative review. J Int AIDS Soc. 2021;24 Suppl 7(Suppl 7):e25793. Epub 2021/11/23. 10.1002/jia2.25793 ; PubMed Central PMCID: PMC8606871. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
10. Walsh SR, Seaman MS. Broadly Neutralizing Antibodies for HIV-1 Prevention. Front Immunol. 2021;12:712122. Epub 2021/08/07. 10.3389/fimmu.2021.712122 ; PubMed Central PMCID: PMC8329589. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
11. Ali MG, Zhang Z, Gao Q, Pan M, Rowan EG, Zhang J. Recent advances in therapeutic applications of neutralizing antibodies for virus infections: an overview. Immunol Res. 2020;68(6):325–39. Epub 2020/11/09. 10.1007/s12026-020-09159-z ; PubMed Central PMCID: PMC7648849. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
12. McCoy LE. The expanding array of HIV broadly neutralizing antibodies. Retrovirology. 2018;15(1):70. Epub 2018/10/18. 10.1186/s12977-018-0453-y ; PubMed Central PMCID: PMC6192334. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
13. Wibmer CK, Moore PL, Morris L. HIV broadly neutralizing antibody targets. Curr Opin HIV AIDS. 2015;10(3):135–43. Epub 2015/03/12. 10.1097/COH.0000000000000153 ; PubMed Central PMCID: PMC4437463. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
14. Rudicell RS, Kwon YD, Ko SY, Pegu A, Louder MK, Georgiev IS, et al.. Enhanced potency of a broadly neutralizing HIV-1 antibody in vitro improves protection against lentiviral infection in vivo. J Virol. 2014;88(21):12669–82. Epub 2014/08/22. 10.1128/JVI.02213-14 ; PubMed Central PMCID: PMC4248941. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
15. Hessell AJ, Rakasz EG, Poignard P, Hangartner L, Landucci G, Forthal DN, et al.. Broadly neutralizing human anti-HIV antibody 2G12 is effective in protection against mucosal SHIV challenge even at low serum neutralizing titers. PLoS Pathog. 2009;5(5):e1000433. Epub 2009/05/14. 10.1371/journal.ppat.1000433 ; PubMed Central PMCID: PMC2674935. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
16. Julg B, Sok D, Schmidt SD, Abbink P, Newman RM, Broge T, et al.. Protective Efficacy of Broadly Neutralizing Antibodies with Incomplete Neutralization Activity against Simian-Human Immunodeficiency Virus in Rhesus Monkeys. J Virol. 2017;91(20). Epub 2017/08/05. 10.1128/JVI.01187-17 ; PubMed Central PMCID: PMC5625479. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
17. Parsons MS, Kristensen AB, Selva KJ, Lee WS, Amarasena T, Esterbauer R, et al.. Protective efficacy of the anti-HIV broadly neutralizing antibody PGT121 in the context of semen exposure. EBioMedicine. 2021;70:103518. Epub 2021/08/14. 10.1016/j.ebiom.2021.103518 ; PubMed Central PMCID: PMC8361295. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
18. Lewitus E, Townsley SM, Li Y, Donofrio GC, Dearlove BL, Bai H, et al.. HIV-1 infections with multiple founders associate with the development of neutralization breadth. PLoS Pathog. 2022;18(3):e1010369. Epub 2022/03/19. 10.1371/journal.ppat.1010369 ; PubMed Central PMCID: PMC8967031. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
19. Townsley SM, Donofrio GC, Jian N, Leggat DJ, Dussupt V, Mendez-Rivera L, et al.. B cell engagement with HIV-1 founder virus envelope predicts development of broadly neutralizing antibodies. Cell Host Microbe. 2021;29(4):564–78.e9. Epub 2021/03/05. 10.1016/j.chom.2021.01.016 ; PubMed Central PMCID: PMC8245051. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
20. Euler Z, van den Kerkhof TL, van Gils MJ, Burger JA, Edo-Matas D, Phung P, et al.. Longitudinal analysis of early HIV-1-specific neutralizing activity in an elite neutralizer and in five patients who developed cross-reactive neutralizing activity. J Virol. 2012;86(4):2045–55. Epub 2011/12/14. 10.1128/JVI.06091-11 ; PubMed Central PMCID: PMC3302374. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
21. Hraber P, Seaman MS, Bailer RT, Mascola JR, Montefiori DC, Korber BT. Prevalence of broadly neutralizing antibody responses during chronic HIV-1 infection. Aids. 2014;28(2):163–9. Epub 2013/12/24. 10.1097/QAD.0000000000000106 ; PubMed Central PMCID: PMC4042313. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
22. Simek MD, Rida W, Priddy FH, Pung P, Carrow E, Laufer DS, et al.. Human immunodeficiency virus type 1 elite neutralizers: individuals with broad and potent neutralizing activity identified by using a high-throughput neutralization assay together with an analytical selection algorithm. J Virol. 2009;83(14):7337–48. Epub 2009/05/15. 10.1128/JVI.00110-09 ; PubMed Central PMCID: PMC2704778. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
23. Corey L, Gilbert PB, Juraska M, Montefiori DC, Morris L, Karuna ST, et al.. Two Randomized Trials of Neutralizing Antibodies to Prevent HIV-1 Acquisition. N Engl J Med. 2021;384(11):1003–14. Epub 2021/03/18. 10.1056/NEJMoa2031738 ; PubMed Central PMCID: PMC8189692. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
24. Bouvin-Pley M, Morgand M, Meyer L, Goujard C, Moreau A, Mouquet H, et al.. Drift of the HIV-1 envelope glycoprotein gp120 toward increased neutralization resistance over the course of the epidemic: a comprehensive study using the most potent and broadly neutralizing monoclonal antibodies. J Virol. 2014;88(23):13910–7. Epub 2014/09/19. 10.1128/JVI.02083-14 ; PubMed Central PMCID: PMC4248973. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
25. Rademeyer C, Korber B, Seaman MS, Giorgi EE, Thebus R, Robles A, 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. Epub 2016/07/21. 10.1371/journal.ppat.1005742 ; PubMed Central PMCID: PMC4951126. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
26. Huang Y, Yu J, Lanzi A, Yao X, Andrews CD, Tsai L, et al.. Engineered Bispecific Antibodies with Exquisite HIV-1-Neutralizing Activity. Cell. 2016;165(7):1621–31. Epub 2016/06/18. 10.1016/j.cell.2016.05.024 ; PubMed Central PMCID: PMC4972332. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
27. Priddy FH, Lewis DJM, Gelderblom HC, Hassanin H, Streatfield C, LaBranche C, et al.. Adeno-associated virus vectored immunoprophylaxis to prevent HIV in healthy adults: a phase 1 randomised controlled trial. Lancet HIV. 2019;6(4):e230–e9. Epub 2019/03/20. 10.1016/S2352-3018(19)30003-7 ; PubMed Central PMCID: PMC6443625. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
28. Julg B, Stephenson KE, Wagh K, Tan SC, Zash R, Walsh S, et al.. Safety and antiviral activity of triple combination broadly neutralizing monoclonal antibody therapy against HIV-1: a phase 1 clinical trial. Nat Med. 2022;28(6):1288–96. Epub 2022/05/14. 10.1038/s41591-022-01815-1 ; PubMed Central PMCID: PMC9205771 interests. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
29. Stephenson KE, Julg B, Tan CS, Zash R, Walsh SR, Rolle CP, et al.. Safety, pharmacokinetics and antiviral activity of PGT121, a broadly neutralizing monoclonal antibody against HIV-1: a randomized, placebo-controlled, phase 1 clinical trial. Nat Med. 2021;27(10):1718–24. Epub 2021/10/09. 10.1038/s41591-021-01509-0 ; PubMed Central PMCID: PMC8516645. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
30. Caskey M, Schoofs T, Gruell H, Settler A, Karagounis T, Kreider EF, et al.. Antibody 10–1074 suppresses viremia in HIV-1-infected individuals. Nat Med. 2017;23(2):185–91. Epub 2017/01/17. 10.1038/nm.4268 ; PubMed Central PMCID: PMC5467219. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
31. Gaudinski MR, Houser KV, Doria-Rose NA, Chen GL, Rothwell RSS, Berkowitz N, et al.. Safety and pharmacokinetics of broadly neutralising human monoclonal antibody VRC07-523LS in healthy adults: a phase 1 dose-escalation clinical trial. Lancet HIV. 2019;6(10):e667–e79. Epub 2019/09/02. 10.1016/S2352-3018(19)30181-X . [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
32. Caskey M, Klein F, Lorenzi JC, Seaman MS, West AP Jr., Buckley N, et al.. Viraemia suppressed in HIV-1-infected humans by broadly neutralizing antibody 3BNC117. Nature. 2015;522(7557):487–91. Epub 2015/04/10. 10.1038/nature14411 ; PubMed Central PMCID: PMC4890714. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
33. Kuriakose Gift S, Wieczorek L, Sanders-Buell E, Zemil M, Molnar S, Donofrio G, et al.. Evolution of Antibody Responses in HIV-1 CRF01_AE Acute Infection: Founder Envelope V1V2 Impacts the Timing and Magnitude of Autologous Neutralizing Antibodies. J Virol. 2023;97(2):e0163522. Epub 2023/02/08. 10.1128/jvi.01635-22 ; PubMed Central PMCID: PMC9973046. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
34. Elangovan R, Jenks M, Yun J, Dickson-Tetteh L, Kirtley S, Hemelaar J. Global and Regional Estimates for Subtype-Specific Therapeutic and Prophylactic HIV-1 Vaccines: A Modeling Study. Front Microbiol. 2021;12:690647. Epub 2021/08/03. 10.3389/fmicb.2021.690647 ; PubMed Central PMCID: PMC8320730. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
35. Korber B, Gaschen B, Yusim K, Thakallapally R, Kesmir C, Detours V. Evolutionary and immunological implications of contemporary HIV-1 variation. Br Med Bull. 2001;58:19–42. Epub 2001/11/21. 10.1093/bmb/58.1.19 . [Abstract] [CrossRef] [Google Scholar]
36. Dreja H, O’Sullivan E, Pade C, Greene KM, Gao H, Aubin K, et al.. Neutralization activity in a geographically diverse East London cohort of human immunodeficiency virus type 1-infected patients: clade C infection results in a stronger and broader humoral immune response than clade B infection. J Gen Virol. 2010;91(Pt 11):2794–803. Epub 2010/08/06. 10.1099/vir.0.024224-0 . [Abstract] [CrossRef] [Google Scholar]
37. Hraber P, Korber BT, Lapedes AS, Bailer RT, Seaman MS, Gao H, et al.. Impact of clade, geography, and age of the epidemic on HIV-1 neutralization by antibodies. J Virol. 2014;88(21):12623–43. Epub 2014/08/22. 10.1128/JVI.01705-14 ; PubMed Central PMCID: PMC4248897. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
38. Bandawe GP, Moore PL, Werner L, Gray ES, Sheward DJ, Madiga M, et al.. Differences in HIV type 1 neutralization breadth in 2 geographically distinct cohorts in Africa. J Infect Dis. 2015;211(9):1461–6. Epub 2014/11/16. 10.1093/infdis/jiu633 ; PubMed Central PMCID: PMC4462656. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
39. Herbeck JT, Peebles K, Edlefsen PT, Rolland M, Murphy JT, Gottlieb GS, et al.. HIV population-level adaptation can rapidly diminish the impact of a partially effective vaccine. Vaccine. 2018;36(4):514–20. Epub 2017/12/16. 10.1016/j.vaccine.2017.12.004 ; PubMed Central PMCID: PMC6701864. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
40. Bouvin-Pley M, Beretta M, Moreau A, Roch E, Essat A, Goujard C, et al.. Evolution of the Envelope Glycoprotein of HIV-1 Clade B toward Higher Infectious Properties over the Course of the Epidemic. J Virol. 2019;93(6). Epub 2018/12/21. 10.1128/JVI.01171-18 ; PubMed Central PMCID: PMC6401426. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
41. Wagh K, Seaman MS, Zingg M, Fitzsimons T, Barouch DH, Burton DR, et al.. Potential of conventional & bispecific broadly neutralizing antibodies for prevention of HIV-1 subtype A, C & D infections. PLoS Pathog. 2018;14(3):e1006860. Epub 2018/03/06. 10.1371/journal.ppat.1006860 ; PubMed Central PMCID: PMC5854441 its variants, for HIV-1 prevention and therapy. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
42. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80. Epub 2013/01/19. 10.1093/molbev/mst010 ; PubMed Central PMCID: PMC3603318. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
43. Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al.. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol Biol Evol. 2020;37(5):1530–4. Epub 2020/02/06. 10.1093/molbev/msaa015 ; PubMed Central PMCID: PMC7182206. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
44. Nickle DC, Heath L, Jensen MA, Gilbert PB, Mullins JI, Kosakovsky Pond SL. HIV-specific probabilistic models of protein evolution. PLoS One. 2007;2(6):e503. Epub 2007/06/07. 10.1371/journal.pone.0000503 ; PubMed Central PMCID: PMC1876811. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
45. Salazar-Gonzalez JF, Salazar MG, Keele BF, Learn GH, Giorgi EE, Li H, et al.. Genetic identity, biological phenotype, and evolutionary pathways of transmitted/founder viruses in acute and early HIV-1 infection. J Exp Med. 2009;206(6):1273–89. Epub 2009/06/03. 10.1084/jem.20090378 ; PubMed Central PMCID: PMC2715054. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
46. Brown BK, Wieczorek L, Sanders-Buell E, Rosa Borges A, Robb ML, Birx DL, et al.. Cross-clade neutralization patterns among HIV-1 strains from the six major clades of the pandemic evaluated and compared in two different models. Virology. 2008;375(2):529–38. Epub 2008/04/25. 10.1016/j.virol.2008.02.022 . [Abstract] [CrossRef] [Google Scholar]
47. Agan BK, Ganesan A, Byrne M, Deiss R, Schofield C, Maves RC, et al.. The US Military HIV Natural History Study: Informing Military HIV Care and Policy for Over 30 Years. Mil Med. 2019;184(Suppl 2):6–17. Epub 2019/11/30. 10.1093/milmed/usy430 ; PubMed Central PMCID: PMC6886375. [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
48. Montefiori DC. Evaluating neutralizing antibodies against HIV, SIV, and SHIV in luciferase reporter gene assays. Curr Protoc Immunol. 2005;Chapter 12:12.1.1-.1.7. Epub 2008/04/25. 10.1002/0471142735.im1211s64 . [Abstract] [CrossRef] [Google Scholar]

Articles from PLOS Pathogens are provided here courtesy of PLOS

Citations & impact 


Impact metrics

Jump to Citations

Alternative metrics

Altmetric item for https://www.altmetric.com/details/157287354
Altmetric
Discover the attention surrounding your research
https://www.altmetric.com/details/157287354

Smart citations by scite.ai
Smart citations by scite.ai include citation statements extracted from the full text of the citing article. The number of the statements may be higher than the number of citations provided by EuropePMC if one paper cites another multiple times or lower if scite has not yet processed some of the citing articles.
Explore citation contexts and check if this article has been supported or disputed.
https://scite.ai/reports/10.1371/journal.ppat.1011780

Supporting
Mentioning
Contrasting
1
0
0

Article citations

Data 


Data behind the article

This data has been text mined from the article, or deposited into data resources.

Similar Articles 


To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.


Funding 


Funders who supported this work.

Division of Intramural Research, National Institute of Allergy and Infectious Diseases (1)

NIAID NIH HHS (1)

U.S. Department of Defense (3)