Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Differentiating founder and chronic HIV envelope sequences

Abstract

Significant progress has been made in characterizing broadly neutralizing antibodies against the HIV envelope glycoprotein Env, but an effective vaccine has proven elusive. Vaccine development would be facilitated if common features of early founder virus required for transmission could be identified. Here we employ a combination of bioinformatic and operations research methods to determine the most prevalent features that distinguish 78 subtype B and 55 subtype C founder Env sequences from an equal number of chronic sequences. There were a number of equivalent optimal networks (based on the fewest covarying amino acid (AA) pairs or a measure of maximal covariance) that separated founders from chronics: 13 pairs for subtype B and 75 for subtype C. Every subtype B optimal solution contained the founder pairs 178–346 Asn-Val, 232–236 Thr-Ser, 240–340 Lys-Lys, 279–315 Asp-Lys, 291–792 Ala-Ile, 322–347 Asp-Thr, 535–620 Leu-Asp, 742–837 Arg-Phe, and 750–836 Asp-Ile; the most common optimal pairs for subtype C were 644–781 Lys-Ala (74 of 75 networks), 133–287 Ala-Gln (73/75) and 307–337 Ile-Gln (73/75). No pair was present in all optimal subtype C solutions highlighting the difficulty in targeting transmission with a single vaccine strain. Relative to the size of its domain (0.35% of Env), the α4β7 binding site occurred most frequently among optimal pairs, especially for subtype C: 4.2% of optimal pairs (1.2% for subtype B). Early sequences from 5 subtype B pre-seroconverters each exhibited at least one clone containing an optimal feature 553–624 (Ser-Asn), 724–747 (Arg-Arg), or 46–293 (Arg-Glu).

Citation: Murray JM, Maher S, Mota T, Suzuki K, Kelleher AD, Center RJ, et al. (2017) Differentiating founder and chronic HIV envelope sequences. PLoS ONE 12(2): e0171572. https://doi.org/10.1371/journal.pone.0171572

Editor: Luis Menéndez-Arias, Universidad Autonoma de Madrid Centro de Biologia Molecular Severo Ochoa, SPAIN

Received: July 21, 2016; Accepted: January 23, 2017; Published: February 10, 2017

Copyright: © 2017 Murray et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: The sequence data are available in the Dryad Data Depository. Data package title: Data from: Differentiating founder and chronic HIV envelope sequences Provisional DOI: doi:10.5061/dryad.r19c2 Data files: HIV envelope sequences Seroconverter HIV subtype B envelope sequences.

Funding: JM was funded by a UNSW Goldstar grant (RG134772, https://research.unsw.edu.au/unsw-internal-funding-opportunities). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

There has been a significant global effort to develop an effective vaccine for HIV. Vaccine trials to date have shown limited efficacy but what success there has been was associated with the ability of the vaccine to stimulate HIV envelope glycoprotein (Env) antibodies [1, 2]. Hence future vaccine candidates will most likely include a component that elicits antibodies specific to targets on Env. However, HIV-1 has an extremely high rate of sequence evolution and strains in different communities form distinct subtypes [3]; high strain diversity exists between individuals and even within any one patient [4]. The challenge for vaccines stimulating antibodies to Env is to target rare common epitopes between viral strains, and ideally between subtypes [1]. These targets should be representative of early founder virus clones that emerge through the strain-selecting bottleneck of transmission [5].

A newly infecting clone quickly undergoes sequence evolution impacted by a developing immune response, as assessed by differences between founder and chronic viral sequences [5–11]. Importantly for vaccines aiming to prevent transmission, HIV infection generally results from a single transmitted virus strain [5, 8], indicating the potential benefits of identifying key transmission-related features. Envelope features are of particular relevance as it is the only HIV protein exposed on the outer surface of an infectious particle, also making it the target for neutralising antibodies [12]. Despite extreme Env sequence diversity, all strains must preserve the functional properties of CD4 receptor binding and binding to chemokine coreceptors such as CCR5 and CXCR4 that facilitate entry. Almost all founder virus uses the CCR5 coreceptor [11], a trait mostly encoded in the V3 loop of Env [13]. Other HIV Env functional domains include motifs allowing interaction with cell adhesion and trafficking receptors, like the Integrin α4β7 glycoprotein. This glycoprotein acts as a gut-homing receptor for lymphocytes, targeting them to the extensive gut associated lymphatic tissues (GALT) that are important sites for explosive viral expansion in the early phase of infection [14–16].

Envelope is the target for broadly neutralizing antibodies (bNAb), and several conserved tertiary structures on Env have been identified as vulnerable sites for bNAb binding [17, 18]. It is unclear whether important new conformational targets for bNAb might be selected during transmission. Understanding what Env features are unique to the transmitted virus but which evolve to escape selection under immune pressure, may point to possible vaccine targets [11, 19].

One approach to identifying vaccine targets is to determine individual amino acids (AA) that differ significantly between chronic and founder sequences among the 857 positions of the Env glycoprotein gp160 [7]. This glycoprotein is cleaved to form the 511 AA CD4 binding subunit gp120, and the 346 AA gp41 that is required for fusion of the virus with the cell membrane. This non-covalently bound heterodimer self-associates into trimers that form the functional Env spike on the viral surface that determines viral tropism. By chance there will be many differences in these sequences, given the large variability in some regions of Env. Functionally related sites can also be compared between these groups, as well as the number of glycosylation sites [7, 9]. This direct comparison of individual positions or known functional sites has proven useful but is limited. How AA and regions within the linear Env genetic sequence determine function is related to their positions in the complex 3-dimensional Env trimer structure [20–22], and determining the interplay between AA in this structure is not straightforward. Any susceptibility in the transmission virus will result in the Env sequence evolving in a series of compensatory escape mutations, so that the trimer structure is altered in order to avoid or minimise the effectiveness of the developing immune response [11, 12]. This collection of escape mutations need not be contiguous in the sequence but will form a biologically related network of positions in Env.

One way of identifying biologically related areas within a highly structured protein involves calculating positions on a set of AA sequences that covary. Pairs of positions are said to covary if the AA combinations observed at these positions are sufficiently different from random combinations. For example if in 16 Env sequences at positions 223 and 432 there were 8 Phe-Lys pairs and 8 Tyr-Arg pairs then positions 223 and 432 would covary since these observed pairs are sufficiently different to the random combinations of 4 Phe-Lys, 4 Phe-Arg, 4 Tyr-Lys and 4 Tyr-Arg pairs. It would also suggest that these positions are linked in some functional manner. If a virus is to evolve under immune pressure, then a single mutation is usually insufficient and a number of compensatory, fitness-restoring mutations at other positions are required [23, 24]. Each of these compensatory mutations results in a covarying pair and the entire set of mutations results in a network of covarying pairs [25]. Analysing sequences using the covariance between amino acid positions is a useful approach in a variety of applications [26–29]. In the case of Env, amino acid sequence networks of covarying pairs contain possible vaccine targets, especially if we can identify those combinations that are transmission signatures present in founder sequences but absent in chronic sequences. Here we follow the approach developed in Murray et al. [30], where optimisation methods were used in combination with covariance calculations to determine the most prominent features that differentiate one hepatitis C virus group from another. We apply these methods to sets of founder and chronic HIV subtype B and subtype C Env sequences, with the aim of identifying features that distinguish founder from chronic virus and that represent possible vaccine targets.

Results

As described previously [9], HIV Env sequences for 133 transmission strain cases (78 subtype B, 55 subtype C) were obtained from Keele et al. and Abrahams et al. [5, 8]. Subtype C is the most common virus in Africa whereas subtype B is prevalent in developed countries such as Australia and the US. These founder sequences represent the inferred single clones that gave rise to productive infection in these individuals [5, 8]. A comparison group of 133 HIV Env sequences were derived from the plasma of individuals with chronic infection, obtained from the Los Alamos National Laboratory (LANL) HIV sequence database. Each randomly selected control sequence was investigated and submitted to rigorous exclusion criteria. Controls were frequency-matched on HIV-1 subtype and geographical location (consistent with the approach of Gnanakaran et al. [7]), with 78 subtype B chronic sequences selected from USA/Trinidad and Tobago and 55 subtype C sequences selected from South Africa/Malawi.

The sequences were aligned and numbered relative to the HXB2 reference strain [31]. The phylogenetic trees for these subtypes and how the Founder and Chronic sequences distribute are shown in Fig 1.

thumbnail
Download:
Fig 1. Unrooted phylogenetic trees for the subtype B and C HIV Env sequences.

Founder sequences are shown with red dashed branches, while chronic sequences are denoted by blue branches.

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

After alignment of the approximately 857 long AA Env sequences, covarying pairs were calculated as previously described [25, 30], to indicate AA positions in Env that are possibly connected through function, where a change in one position is likely to also result in a change in another position. In this study we analyse both subtype B and C viruses, hence it was necessary to perform the covariance calculations separately for each subtype. There were 2,495 covarying pairs over the 156 subtype B sequences and 3,021 subtype C covarying pairs over the 110 subtype C sequences. The maximum covariance values were similar for each of the subtypes, 14.0 for subtype B and 14.6 for subtype C, however there were more covarying pairs for subtype C that were close to this maximum value (Fig 2). The slightly higher covariance among subtype C, despite fewer sequences, tends to reflect the more clustered phylogenetic tree. This has been observed previously with differences in covariance between 1a and 1b HCV sequences mirroring the clustering of these sequences [30].

thumbnail
Download:
Fig 2. Conserved and covarying regions for subtypes B and C Env.

Regions are coloured as conserved across both subtypes (black), conserved within each subtype (dark blue), conserved except for a maximum of 2 individuals in that subtype (light blue), and covarying (magenta). Those covarying pairs with at least 20% of the maximum covariance value for subtype C, and 12.5% for subtype B, are connected with magenta lines. The different levels of covariance were determined to include approximately the same numbers of covarying pairs in each case: 78 for subtype B and 79 for subtype C. The signal, α4β7 binding site, constant (C1-C6) and variable (V1-V5) regions within gp120, and the gp41 domain are mapped onto the Env sequence.

https://doi.org/10.1371/journal.pone.0171572.g002

Covarying AA contain possible targets for a vaccine-stimulated response since they are not totally conserved, but are not so variable as to represent random changes. We mapped regions that were conserved, covarying and variable in Env over both subtypes (Fig 2). For subtype C, highly covarying pairs clustered in the signal, V2, and C5 regions of gp120 and in gp41. The pair with the highest covariance was 474–476, while the positions that appeared most were 192 (22 times), 476 (12), 27 (12), 11 (11), 706 (11), 388 (10), 595 (10). Covarying pairs were incident (one of the positions within the pair) 11 times to the α4β7 binding site (positions 179 to 181). For subtype B the pair with the highest covariance was 230–232, connecting within an N-linked glycosylation site (NXS/T at positions 230, 231, 232), and the other highly covarying pairs tended to be in close proximity as well. One of these covarying pairs (181–693), was also incident to the α4β7 binding site.

Optimal networks

The covarying pairs determine components within the sequences that vary but not in a random fashion. However as depicted in Fig 2, there can be many of these. Our aim was to determine what is special in founder viruses. We define a separating pair as a covarying pair of positions and its AA combination, where at that pair of positions, the AA combination is exhibited by some sequence in the founder group but by none in the chronic group. Hence we search among separating pairs to determine signatures of founder virus that are not expressed by any chronic virus in our samples. To determine those that are most pertinent to characterising founder virus we calculated optimal collections of separating pairs. Optimality was based on requiring the fewest pairs to separate these groups or that tended to maximise covariance (Methods). Separation of Founders from Chronics becomes more difficult as fewer pairs are used. Hence the solution with the fewest pairs will likely contain stronger predictors of what differs between these groups. This is similar in essence to identifying AA and their positions where comparisons between the groups gives the lowest p or q values [7]. Initially a single optimal network was identified for each of the subtype/covariance group combinations (covariance over all, founder, or chronic sequences), and objective function type (fewest pairs or optimising covariance). Hence, twelve optimal networks were identified and used in this part of the analysis.

Optimization enforces two limitations on this problem. Firstly it can only be applied to a single alignment rather than a bootstrapped collection of alignments. To compensate for this limitation we only included covariance calculated over positions that contained at most 10% gaps. This had the effect of excluding some of the more variable regions around indels. Secondly this binary integer programming problem belongs to the class of Non-Polynomially (NP) Hard problems that are characterised by small increases in problem size (here this is related to the number of covarying pairs) resulting in large increases in computational time. Partly for this reason we limited the size of the problem by only including pairs with a covariance value of at least 0.5, leading to the 2,495 covarying pairs for subtype B described above. Including pairs with a covariance of 0.1 would have led to an almost 10-fold increase in the number of pairs (21,937 for subtype B). The difficulty of solving these problems with lower covariance cut-off is highlighted by the related connected optimal network problem being computationally intractable in reasonable time even with the 0.5 cut-off [32]. Nevertheless to address sensitivity to this assumption, we also investigate in a later section a slightly simpler problem but with no covariance cut-off.

The optimal networks, determining features of Founder sequences that separate them from Chronic sequences, required between 13 and 15 AA pairs for subtype B (S1 Table). The fewer subtype C sequences required between 11 and 13 AA pairs to separate all Founder sequences from Chronic sequences. Some of the pairs within the optimal networks appeared multiple times for the 6 optimization problems in each subtype/group combination (3 minimizing the number of pairs, and 3 optimizing a measure of total covariance). For subtype B the pairs 278–620 (Ser Asp), and 750–836 (Asp Ile) each appeared 4 times, while for subtype C the most frequently appearing pairs (3 times each) were 170–192 (His Arg) and 588–662 (Arg Ala) (Table 1).

thumbnail
Download:
Table 1. Pairs observed multiple times (frequency f) in optimal networks for each subtype.

The number of individuals exhibiting each AA combination is denoted by n.

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

Several single positions within these Env sequences appeared frequently in optimal networks (Table 2). Positions 836 and 620 for subtype B, and positions 192 and 346 for subtype C were the most frequently appearing single positions. Position 346 for subtype C was also proximal to position 347 appearing 6 times for subtype B in the constant C3 region of gp120. The next most frequently appearing proximal positions were 178 for subtype B and 179 for subtype C that were within or next to the α4β7 binding site.

thumbnail
Download:
Table 2. Single AA in optimal networks determined on Founders, which appear at least 2 times (frequency f).

The region of Env is denoted ahead of a decimal point, and any recognised motif after the decimal.

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

Multiple optimal solutions for each subtype/group combination

The optimization procedure above determines an optimal solution but these are not unique. For example, the first row of S1 Table shows that an optimal solution for the problem where we minimize the number of pairs, contains 13 pairs for subtype B and 11 pairs for subtype C, when covariance is calculated over all sequences within each subtype. However a total of 13 different optimal networks that each differ by at least one AA pair for subtype B can be identified for this particular problem (Table 3). Similarly, for the subtype C sequences there are 75 optimal solutions. However these optimal solutions share several features. For subtype B there are 9 pairs that are present in each of the 13 optimal solutions, indicating possibly susceptible immune-evasion pathways. For subtype C, 644–781 (Lys Ala) is present in 74 solutions, while 2 pairs appear in 73 solutions, and 6 pairs in 72 solutions. It is interesting to note that there are no pairs for subtype C that are present in all 75 optimal solutions, suggesting a greater number of pathways by which this virus evolves and possibly making it a more difficult vaccine target.

thumbnail
Download:
Table 3. The pairs that are observed in a given number of optimal solutions for subtypes B and C.

Listed for each of the optimal separating pairs are the covarying positions, the amino acids for the sequences in the Founder separating pairs, and then the Env motifs for each of these positions. The optimization problem was solved using objective i) (minimizing the number of pairs) on a covariance network constructed using all sequences in each subtype.

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

Optimal networks with no covariance restriction

By restricting pairs to those with a covariance value of 0.5 or higher, we attempted to determine functionally relevant positions. However this excluded many other possibilities. Allowing all pairs to be considered in separating Founder from Chronics, increased the number of pairs more than 50-fold, making the calculation of optimal networks with criterion (ii) that incorporated the covariance value computationally impractical. However we could still determine optimal networks that achieved separation of Founders from Chronics with the fewest number of pairs. The optimal subtype B network contained 12 separating pairs: 46–293 (Arg Glu), 84–333 (Val Val), 269–767 (Asp Lys), 278–620 (Ser Asp), 291–758 (Ala Val), 293–375 (Val Thr), 336–535 (Glu Met), 345–842 (Ile Asn), 535–624 (Lys Asp), 553–624 (Ser Asn), 724–747 (Arg Arg), 750–836 (Asp Ile). Despite including many more pairs, this optimal network contained only one or two fewer optimal pairs than with the above S value cut-off. Moreover 750–836 (Asp Ile) appeared in all multiple optimal solutions (Table 3).

The subtype C optimal solution without covariance restriction consisted of 9 pairs: 9–515 (Asn Met), 12–350 (Gln Ser), 166–352 (Lys His), 330–620 (Tyr Thr), 352–379 (Tyr Gln), 448–821 (Ser Ala), 565–588 (Met Arg), 640–778 (Asn Val), 833–837 (Leu Cys).

Occurrence of optimal pairs in virus cloned prior to seroconversion

The above calculations determined aspects of founder sequences that might be targeted by vaccines. Given the large number of covarying pairs determined in this way, it is unlikely that all of these pairs would be robust predictors of susceptibility. To test this we sequenced subtype B Env virus from 5 individuals who were newly infected (Methods). The optimal pairs (for all subtype B solutions in S1 Table) showed little overlap with the clones derived from these pre-seroconverters (Table 4). Part of this limitation was due to the cloning process only including positions 44 to 752. Three of the five individuals and 5 of the 12 clones exhibited some of the optimal pairs but these were limited to: 343–621 (Gln Asp) and 354–636 (Pro Asp) for one clone; 553–624 (Ser Asn) for 2 individuals; 624–747 (Glu Arg) for 2 individuals, and 724–747 (Arg Arg) for one individual.

thumbnail
Download:
Table 4. Comparison of optimal Founder pairs with the AA combinations appearing for 5 subtype B pre-seroconverters with each of the clones sequenced.

https://doi.org/10.1371/journal.pone.0171572.t004

The two pairs 553–624 (Ser Asn) and 724–747 (Arg Arg) were also in the optimal network determined for subtype B when there was no covariance restriction. After including the 46–293 (Arg Glu) optimal pair, these three features were exhibited by at least one clone from each seroconverter (Table 4).

Discussion

By searching over covarying pairs that exhibit covariance above a minimal background level, we are attempting to determine aspects of Env that exhibit some change, and hence are likely to reflect, and to be susceptible to, an evolving immune response. Optimal pairs, selected either by being part of a network with the fewest pairs, or maximizing a measure of covariance, in some sense represent fundamental aspects of what separates Founder from Chronic virus, since separation cannot be achieved using fewer pairs or with a higher measure of covariance. This is analogous to determining AA with the lowest p or q values in a comparison at individual positions of Env [7]. As with these statistically significant AA, the optimal pairs may indicate components within transmitted viruses that determine selection at the transmission bottleneck and that can be targeted by bNAbs. Although there are a large number of different optimal pairs when collected over all problems (60 subtype B and 54 subtype C, determined from the unique pairs contained within S1 Table), these are not randomly distributed over Env (Fig 3). The majority are incident to gp41 (54 and 36 pairs respectively have at least one position in gp41), and occur at regions related to viral fusion and at AA internal to the membrane involved in virion structure and intracellular signalling. There are 20 and 22 optimal pairs respectively contained entirely within gp120 and that are incident to the C3 region, which supports the view that these calculations are not entirely driven by the variable regions. While the next highest region with incident pairs in subtype B is the C2 region (19 pairs) this was not frequent for subtype C (5 pairs) which instead had high incidence in V2 (19 pairs, 7 pairs for subtype B). Several identified covarying founder virus sites beyond 665 exist within or beneath the membrane layer and would not directly bind antibody.

thumbnail
Download:
Fig 3. The collection of optimal pairs displayed relative to the domains of Env.

A) over a linear representation and B) as networks.

https://doi.org/10.1371/journal.pone.0171572.g003

The gut-homing α4β7 integrin has been implicated in higher susceptibility to HIV infection [42] and α4β7 high CD4+ T cells have been suggested as preferential targets in mucosal transmission especially for subtype C [15]. The viral α4β7 binding site is located within V2 at positions 179 to 181. Despite only covering 3 of the 857 AA in Env, optimal pairs incident to this domain were particularly prevalent but especially so for subtype C: 4.2% of optimal pairs (S1 Table) relative to 0.35% of Env and a ‘prevalence ratio’ of 11.9 (4.2%/0.35%). Of interest is that several α4β7 binding site optimal pairs were also incident to the external fusion or membrane internal domains of Env gp41. While structural connections between V2 and the gp41 fusion domains can be envisioned from the crystal and cryo-EM structures of Env gp160 [20, 22], these are not anticipated across the membrane, and suggest a functional connection, possibly in cell-cell viral transfer. The remainder of V2 was also prevalent in these optimal pairs for subtype C with a prevalence ratio of 4.6. The α4β7 binding site was less likely to be represented in this analysis for subtype B with a prevalence ratio of 3.3, but still higher than the other domains. However these results are in contrast to findings that blocking binding to α4β7 had little impact on infection by founder subtype C virus [43], or that HIV Env is in general a poor ligand for this integrin [44]. Nevertheless monoclonal antibodies that targeted α4β7 integrin-expressing CD4+ T cells protected rhesus macaques from SIV infection, possibly suggesting a large role played by early virus in binding to these cells [45], although it is unclear whether this effect was through altered Env binding and entry, or through altered gut mucosal trafficking. The subtype C α4β7 pairs were also incident to the domains: signal (1), V2 (1), C3 (2) and gp41 (2). Of note is that the pairs between α4β7 and C3, 181–346 (Val Gln) and 181–364 (Ile Ala), are incident to N-linked glycosylation sites (339–341 or 344–346, 362–364). Glycosylation sites in C3 have been observed to affect α4β7 reactivity [15].

AA within gp41 or the α2 helix in C3 (positions 335 to 353), tended to be central to the networks of all optimal pairs for each subtype (Fig 3B), exhibiting the greatest node degrees within those networks. In particular the C3 α2 helix positions formed a strongly connected subnetwork for subtype C that also included V2 and α4β7 positions, suggesting these Founder characteristics evolved under immune pressure in a related manner. This is consistent with observations of the staged development of antibodies targeting the α2 helix of C3 followed by those targeting V1-V2 for subtype C [46].

We identified all optimal solutions relative to the problem of minimizing the number of pairs with covariance over all sequences, there being many more for subtype C than for subtype B (Table 3) possibly related to the higher covariance network exhibited by subtype C (Fig 2) and perhaps reflecting a more diverse pathway of evolution for this genotype. Every subtype B optimal solution contained a core set of 9 covarying pairs some of which have been identified through their interactions with bNAbs. Glycans at positions 234 and 276 are essential for the reactivity of bNAb 8ANC195 [18], where these would be impacted by linked mutations at nearby sites in the optimal pairs 232–236 and 279–315, that lie at the gp120-gp41 ectodomain structural interface [20, 22, 47]. Mutations at positions 89, 90, 227, 232 and 243 also diminished neutralization by bNAb 35O22 [48]. Ten of the 39 founder sequences contained the 232–236 Thr Ser combination, while none of the chronic sequences displayed this feature, revealing potential differences in N-linked glycosylation that may optimise the viral entry functions involving the gp120-gp41 interface, but expose a neutralisation sensitive epitope.

As may be expected by the much larger number of optimal solutions for subtype C, there were no pairs that were common to all solutions. This would suggest that there may be no single Env target that will be effective in blocking transmission to an individual in a community where subtype C is prevalent. Investigation of glycans for an elite subtype C neutralizer showed N611D and N637K as well as E647A had the greatest effect on neutralization which may be relevant to the 644–781 pair appearing in 74/75 optimal solutions [49]. That the covarying AA lie beneath the membrane and are shielded from antibody access, raises the possibility of either structural alterations across the membrane, or effects on cell-cell transmission involving intracellular molecular interactions. One of the glycosylation sites that varied across subtype C virus commenced at position 133 [50], part of the optimal pair in 73 of 75 optimal solutions. Positions on the α2 helix of gp120: 236, 305, 332, 335, 336, 337, 343, 350 and 393 are strongly selected in subtype C for resistance to NAb [51], which may be relevant to the appearance in multiple optimal solutions of 10–346 Tyr Gly (72/75) and 307–337 Ile Gln (73/74).

Testing whether these optimal subtype B Founder pairs were also represented in early Env sequences from 5 pre-seroconverters showed little overlap: only 2 pairs connecting C3 to gp41, and three within gp41 with 2 of these involving gp41 sequence inside the membrane (Table 4). However one strength of this analysis is that the pairs that did overlap with the seroconverter sequences were not contained in any of the 78 Chronic sequences. Only positions 44 to 752 of Env were sequenced in the cloning process which therefore omitted comparison with some of the optimal pairs. These pre-seroconverter clones may also be variants from the virus that established the infections approximately 20 days earlier. Effects on cell signalling or cell-cell transmission may have promoted the identification of AA located beneath the membrane. Nevertheless identified pairs incident to glycosylation sites were most likely to play a role with these occurring in the pre-seroconverters at or near positions 354 in C3 and 624 in gp41. A comparison of the seroconverter sequences with the 12 optimal pairs determined for subtype B with no covariance restriction revealed 3 features that were exhibited by at least one clone in each individual: 553–624 (Ser Asn), 724–747 (Arg Arg) (these two were identified in the above analysis) and 46–293 (Arg Glu).

The optimization and comparison calculations determining differences between Founders and Chronics can only be performed on a single alignment. Given the high degree of variability and the presence of indels in these sequences there will be a number of possible alignments with differences occurring mainly at, or adjacent to, gaps in the alignment. Partly because of this we omitted positions where more than 10% of sequences contained gaps or were uncertain. Including these positions leads to optimal separating networks more highly incident to indels (data not shown) and as such less likely to be robust targets for vaccine stimulated antibodies. Omitting positions that were predominantly gaps meant that we were less able to investigate differences in length of some Env regions between Founders from Chronics as previously observed by us and others [9, 52]. However the number of positions omitted were not markedly different between Founders vs Chronics—for subtype B there were 26 positions in the Founder sequences and not in Chronics consisting entirely of gaps and 18 positions in Chronics but not in Founders (20 vs 17 respectively for subtype C). Dropping positions also excluded 10 of the 31 HXB2 glycosylation sites so that our analysis will not be able to use changes in their number as a factor, but which we had previously analysed [9].

A further limitation of our approach is that these sequences were not matched for factors apart from subtype and geographical region. Matching Founders and Chronics by gender, transmission mode, etc., would assist in removing extraneous components that confound separation due to vaccine-relevant factors, but it would reduce the power of the analysis and the generality of the results. Our general approach of searching among covarying positions, also assumes that mutation away from the transmitted virus occurs along a few pathways in response to immune pressure and that then induce covariance. This is generally true for the mutational pathways within HIV resulting from the development of drug resistance to a particular antiretroviral agent. Founder envelope sequences have also been observed to co-evolve with broadly neutralizing antibody [12, 46], which may indicate that initially susceptible features in Env can be deduced from the covariance induced by the resulting mutational pathways. However all covariances within Env need not be related to changes due to immune pressure so that the differences we determine between Founders and Chronics can encompass aspects that will not be relevant to vaccine targets.

The covarying pairs link regions within Env. Linkage across domains, which is not surprising given the complex trimer structure which brings variable and conserved regions into close proximity [22], is known to impact on infectivity and function [53, 54]. Some of the linkages determined here may be due to the similar geographical regions from which the sequences originated. Although chronic sequences were frequency-matched by geographical location, this cannot be completely ruled out as contributing to some aspects of the optimal networks determined here. The optimal pairs are likely components of more diverse networks that may more precisely describe the multiple binding sites of any bNAbs or the mutational pathways that the virus follows to evade them. The effects measured on gp41 AA beneath the membrane could also arise from emergence of T cell responses. How these larger networks can be extracted from this analysis is a more complex problem. It would also be beneficial to investigate these in silico results in an in vitro setting.

In summary we have used operations research methods to determine the most prominent features of Env that differ between founder and chronic sequences for subtype B and subtype C. Our results suggest that the gut-homing α4β7 integrin plays a role in establishing infection and may indicate key AAs desirable for inclusion in vaccine strains. Unlike subtype B where 9 AA pairs were in all optimal solutions, no single AA pair was present in all subtype C optimal solutions (Table 3). This may highlight difficulties in targeting transmission with a single vaccine strain, especially for subtype C.

Materials and methods

Founder and chronic DNA sequences, as well as the HXB2 reference envelope sequence [55], were converted to AA sequences (nt2aa, MATLAB 2012b, The MathWorks Inc., Natick MA, USA), and then aligned using a progressive multiple alignment method (multialign). Pairwise distances were calculated with the Jukes-Cantor method, with the phylogenetic tree generated using the Unweighted Pair Group Method Average (seqlinkage). AA positions in the aligned sequences were numbered relative to HXB2 according to the convention of Korber et al. [31].

Pre-seroconverters

Env was sequenced from an additional set of subtype B virus for 5 newly infected males (estimated 20 to 23 days from transmission, Fiebig Stage II (n = 4), and III (n = 1)) [56], where transmission was through men having sex with men. Due to the cloning process only positions 44 to 752 of Env were sequenced from these transmission HIV strains.

These pre-seroconverter individuals were enrolled in a naturally history cohort study, the Primary HIV and Early Disease Research: Australian Cohort (PHAEDRA), that was established by the National Centre in HIV Epidemiology and Clinical Research to monitor immunological and virological characteristics of individuals with acute and early HIV-1 infection. Research ethics approval (number 02244) was given by St Vincent's Hospital, Sydney, Research Ethics Committee. All participants signed an informed consent form before study entry.

Networks

The construction of a covariance network was previously described by Murray et al. [30], such that the nodes of the network are given by each AA position, and each pair of covarying positions above the cut-off value provides an edge. Subtype B and C sequences were investigated separately, along with three different sequence groups: i) all sequences within the subtype (All), ii) the founder sequences and iii) the chronic sequences. As such, for a given subtype and sequence group all covarying positions are contained in the set Psubtype,group. A network was then constructed for each subtype and sequence group combination. We calculated optimal networks based on each of these sets, extracting the “best” separating pairs (a separating pair is defined as a covarying pair of positions and an AA combination present in at least one of the founder sequences but in none of the chronic sequences), according to certain criteria and where each founder sequence contained at least one of these separating pairs. The optimal separation of Founders from Chronics was performed using two criteria: i) separating with the fewest pairs, and ii) using a weighting that simultaneously minimized the number of pairs and maximized a measure of total covariance. The latter objective was achieved by applying the weight wk = (Ŝ-Sk)2 to the cost of pair k, where Ŝ is the integer ceiling of the maximum of all covariance values.

For these calculations we considered separating pairs for each AA combination so that optimal solutions will extract single AA combinations at each pair of positions in the optimal network [30].

As an example, the problem PB,All with optimal criterion (i) determines the smallest set of AA combinations expressed by Founder sequences at covarying pairs where these were calculated on all subtype B sequences. In this instance, an optimal solution consists of the set of 13 pairs and AA described by {2–92 (Arg Lys), 24–293 (Ile Lys), 166–836 (Arg Thr), 178–535 (Asn Val), 232–236 (Thr Ser), 240–340 (Lys Lys), 279–315 (Asp Lys), 291–792 (Ala Ile), 322–347 (Asp Thr), 535–620 (Lys Asp), 612–836 (Ala Lys), 724–837 (Arg Phe), 750–836 (Asp Ile)}, where each Founder sequence contains at least one (and possibly more) of these features, while no Chronic sequence exhibits any of these features. In this way this set of AA pairs separates all Founders from all Chronics. The optimality aspect guarantees that there are no other combinations that separate the groups with fewer combinations, although there may be other sets with the same number of pairs and hence are also optimal.

The binary integer programming method used to extract this smallest set from all separating pairs (formulated as described in Murray et al. [30]), is a standard optimization procedure, where a variable xi is assigned to the ith separating pair and given the value 1 if it is included in the optimal set and value 0 otherwise. For this problem optimality is determined by allocating the fewest 1 values among the xi.

Multiple optimal solutions for each subtype and group combination

These problems determining optimal separating pairs do not necessarily have unique solutions. Multiple optimal solutions were identified using an iterative approach. Initially, the problem was solved to identify an optimal solution. A constraint was then added to the problem to exclude the current solution. The problem was then resolved to identify another optimal solution. For example, after calculating the optimal solution above, the separating pair 2–92 (Arg Lys) was excluded and the problem resolved. The solution to this restricted problem also contained 13 pairs, and so comprised an additional optimal solution. On the other hand, excluding 750–836 (Asp Ile) resulted in a solution requiring 14 pairs to separate the two groups and so any optimal solution must include this combination.

Optimal solutions with no lower bound on covariance

The optimal solutions above were calculated over all pairs with a covariance value of at least 0.5. This implied that the pairs might exhibit some functional relationship between them otherwise the AAs at these positions would more likely combine in a random manner. Different cut-off values would allow more or fewer covarying pairs which would impact the optimal solutions. To determine how much this could change we calculated optimal networks, only using criterion (i) that minimized the number of pairs, over all pairs regardless of covariance. The size of these problems, 177,906 subtype B pairs and 157,641 subtype C pairs, made calculating the optimal solution with criterion (ii), impractical for this NP Hard problem. As long as we are not interested in generating connected networks where the AA match at the connecting position [32], we can reduce the number of separating pairs by only including those that are not dominated by others—a separating pair is said to be dominant if it is expressed by a set of Founders IG and where there is no other separating pair which is expressed by Founders IG’ such that IG ⊆ IG’., For example if the separating pair 180Ser-230Asn is exhibited by Founders 1, 2 and 3 while 130Ala-485Tyr is exhibited by Founders 1, 2, 3 and 5, then the first pair is excluded from the calculations of separating Founders from Chronics with the fewest number of pairs. The pair 130Ala-485Tyr dominates 180Ser-230Asn. After this calculation there are only 1,358 dominant pairs for subtype B Founders and 928 dominant pairs for subtype C. The binary integer programming method above can then be applied to this problem with no covariance restriction.

All calculations were performed with Matlab version R2012b (The MathWorks Inc., Natick MA, USA). The binary integer programming problems were solved using the bintprog routine and the CPLEX toolbox (IBM, Armonk NY, USA).

Supporting information

S1 Table. Pairs of AA in optimal networks that separate Founders from Chronics.

Each item lists the pairs in the optimal network when calculations are performed over covariance calculations determined on sequences in All, Founders or Chronics (Sep. Set). These are features exhibited by some founder sequences but by no chronic sequence. The number of sequences that exhibit that feature for that AA pair are listed as (n). Optimality was determined either through choosing the fewest number of pairs (Prob = No) or through maximizing a measure of total covariance (Prob = Yes).

https://doi.org/10.1371/journal.pone.0171572.s001

(DOCX)

Acknowledgments

We thank PJ Klasse for helpful discussions. The authors would also like to thank the patients and clinicians participating in the PHAEDRA study.

Author Contributions

  1. Conceptualization: JMM.
  2. Data curation: TM DP RJC KS ADK.
  3. Formal analysis: JMM SM.
  4. Investigation: TM DP RJC KS ADK.
  5. Methodology: JMM.
  6. Project administration: JMM DP.
  7. Resources: DP ADK.
  8. Software: JMM SM.
  9. Validation: JMM SM.
  10. Visualization: JMM.
  11. Writing – original draft: JMM SM DP.
  12. Writing – review & editing: JMM SM RJC DP.

References

  1. 1. Day TA, Kublin JG. Lessons Learned from HIV Vaccine Clinical Efficacy Trials. Current HIV research. 2013;11(6):441–9. pmid:24033299
  2. 2. Rerks-Ngarm S, Pitisuttithum P, Nitayaphan S, Kaewkungwal J, Chiu J, Paris R, et al. Vaccination with ALVAC and AIDSVAX to Prevent HIV-1 Infection in Thailand. N Engl J Med. 2009;361(23):2209–20. pmid:19843557
  3. 3. Geretti AM. HIV-1 subtypes: epidemiology and significance for HIV management. Curr Opin Infect Dis. 2006;19(1):1–7. pmid:16374210
  4. 4. Shankarappa R, Margolick JB, Gange SJ, Rodrigo AG, Upchurch D, Farzadegan H, et al. Consistent Viral Evolutionary Changes Associated with the Progression of Human Immunodeficiency Virus Type 1 Infection. J Virol. 1999;73(12):10489–502. pmid:10559367
  5. 5. Abrahams M-R, Anderson JA, Giorgi EE, Seoighe C, Mlisana K, Ping L-H, et al. Quantitating the Multiplicity of Infection with Human Immunodeficiency Virus Type 1 Subtype C Reveals a Non-Poisson Distribution of Transmitted Variants. J Virol. 2009;83(8):3556–67. pmid:19193811
  6. 6. Asmal M, Hellmann I, Liu W, Keele BF, Perelson AS, Bhattacharya T, et al. A Signature in HIV-1 Envelope Leader Peptide Associated with Transition from Acute to Chronic Infection Impacts Envelope Processing and Infectivity. PLoS ONE. 2011;6(8):e23673. pmid:21876761
  7. 7. Gnanakaran S, Bhattacharya T, Daniels M, Keele BF, Hraber PT, Lapedes AS, et al. Recurrent Signature Patterns in HIV-1 B Clade Envelope Glycoproteins Associated with either Early or Chronic Infections. PLoS Pathog. 2011;7(9):e1002209. pmid:21980282
  8. 8. Keele BF, Giorgi EE, Salazar-Gonzalez JF, Decker JM, Pham KT, Salazar MG, et al. Identification and characterization of transmitted and early founder virus envelopes in primary HIV-1 infection. Proceedings of the National Academy of Sciences. 2008;105(21):7552–7.
  9. 9. Mota T, Murray J, Center R, Purcell DF, McCaw J. Application of a case-control study design to investigate genotypic signatures of HIV-1 transmission. Retrovirology. 2012;9(1):54.
  10. 10. Sagar M, Laeyendecker O, Lee S, Gamiel J, Wawer MJ, Gray RH, et al. Selection of HIV variants with signature genotypic characteristics during heterosexual transmission. J Infect Dis. 2009;199(4):580–9. pmid:19143562
  11. 11. 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. The Journal of Experimental Medicine. 2009;206(6):1273–89. pmid:19487424
  12. 12. Liao H-X, Lynch R, Zhou T, Gao F, Alam SM, Boyd SD, et al. Co-evolution of a broadly neutralizing HIV-1 antibody and founder virus. Nature. 2013;advance online publication. http://www.nature.com/nature/journal/vaop/ncurrent/abs/nature12053.html#supplementary-information.
  13. 13. Cocchi F, DeVico AL, Garzino-Demo A, Cara A, Gallo RC, Lusso P. The V3 domain of the HIV-1 gp120 envelope glycoprotein is critical for chemokine-mediated blockade of infection. Nat Med. 1996;2(11):1244–7. Epub 1996/11/01. pmid:8898753
  14. 14. Arthos J, Cicala C, Martinelli E, Macleod K, Van Ryk D, Wei D, et al. HIV-1 envelope protein binds to and signals through integrin α4β7, the gut mucosal homing receptor for peripheral T cells. Nat Immunol. 2008;9(3):301–9. http://www.nature.com/ni/journal/v9/n3/suppinfo/ni1566_S1.html. pmid:18264102
  15. 15. Nawaz F, Cicala C, Van Ryk D, Block KE, Jelicic K, McNally JP, et al. The Genotype of Early-Transmitting HIV gp120s Promotes α4β7 –Reactivity, Revealing α4β7+CD4+ T cells As Key Targets in Mucosal Transmission. PLoS Pathog. 2011;7(2):e1001301. pmid:21383973
  16. 16. Wang X, Xu H, Gill AF, Pahar B, Kempf D, Rasmussen T, et al. Monitoring α4β7 integrin expression on circulating CD4+ T cells as a surrogate marker for tracking intestinal CD4+ T-cell loss in SIV infection. Mucosal immunology. 2009;2(6):518–26. pmid:19710637
  17. 17. McLellan JS, Pancera M, Carrico C, Gorman J, Julien J-P, Khayat R, et al. Structure of HIV-1 gp120 V1/V2 domain with broadly neutralizing antibody PG9. Nature. 2011;480(7377):336–43. http://www.nature.com/nature/journal/v480/n7377/abs/nature10696.html#supplementary-information. pmid:22113616
  18. 18. Scharf L, Scheid Johannes F, Lee Jeong H, West Anthony P Jr, Chen C, Gao H, et al. Antibody 8ANC195 Reveals a Site of Broad Vulnerability on the HIV-1 Envelope Spike. Cell Reports. 2014;7(3):785–95. http://dx.doi.org/10.1016/j.celrep.2014.04.001. pmid:24767986
  19. 19. Keele BF. Identifying and characterizing recently transmitted viruses. Current Opinion in HIV and AIDS. 2010;5(4):327. pmid:20543609
  20. 20. Julien J-P, Cupo A, Sok D, Stanfield RL, Lyumkis D, Deller MC, et al. Crystal Structure of a Soluble Cleaved HIV-1 Envelope Trimer. Science. 2013.
  21. 21. Lyumkis D, Julien J-P, de Val N, Cupo A, Potter CS, Klasse P-J, et al. Cryo-EM Structure of a Fully Glycosylated Soluble Cleaved HIV-1 Envelope Trimer. Science. 2013.
  22. 22. Pancera M, Zhou T, Druz A, Georgiev IS, Soto C, Gorman J, et al. Structure and immune recognition of trimeric pre-fusion HIV-1 Env. Nature. 2014;514(7523):455–61. http://www.nature.com/nature/journal/v514/n7523/abs/nature13808.html#supplementary-information. pmid:25296255
  23. 23. Schneidewind A, Brockman MA, Yang R, Adam RI, Li B, Le Gall S, et al. Escape from the Dominant HLA-B27-Restricted Cytotoxic T-Lymphocyte Response in Gag Is Associated with a Dramatic Reduction in Human Immunodeficiency Virus Type 1 Replication. J Virol. 2007;81(22):12382–93. pmid:17804494
  24. 24. Kelleher AD, Long C, Holmes EC, Allen RL, Wilson J, Conlon C, et al. Clustered Mutations in HIV-1 Gag Are Consistently Required for Escape from Hla-B27–Restricted Cytotoxic T Lymphocyte Responses. The Journal of Experimental Medicine. 2001;193(3):375–86. pmid:11157057
  25. 25. Aurora R, Donlin MJ, Cannon NA, Tavis JE. Genome-wide hepatitis C virus amino acid covariance networks can predict response to antiviral therapy in humans. The Journal of Clinical Investigation. 2009;119(1):225–36. pmid:19104147
  26. 26. Altschuh D, Lesk AM, Bloomer AC, Klug A. Correlation of co-ordinated amino acid substitutions with function in viruses related to tobacco mosaic virus. Journal of Molecular Biology. 1987;193(4):693–707. pmid:3612789
  27. 27. Eyal E, Frenkel-Morgenstern M, Sobolev V, Pietrokovski S. A pair-to-pair amino acids substitution matrix and its applications for protein structure prediction. Proteins: Structure, Function, and Bioinformatics. 2007;67(1):142–53.
  28. 28. Gobel U, Sander C, Schneider R, Valencia A. Correlated mutations and residue contacts in proteins. Proteins: Structure, Function, and Bioinformatics. 1994;18(4):309–17.
  29. 29. Lee B-C, Park K, Kim D. Analysis of the residue?residue coevolution network and the functionally important residues in proteins. Proteins: Structure, Function, and Bioinformatics. 2008;72(3):863–72.
  30. 30. Murray JM, Moenne-Loccoz R, Velay A, Habersetzer F, Doffoël M, Gut J-P, et al. Genotype 1 Hepatitis C Virus Envelope Features That Determine Antiviral Response Assessed through Optimal Covariance Networks. PLoS ONE. 2013;8(6):e67254. pmid:23840641
  31. 31. Korber BT, Foley BT, Kuiken CL, Pillai SK, Sodroski JG. Numbering Positions in HIV Relative to HXB2CG. In: Korber B, KC L., Foley B, Hahn B, McCutchan F, Mellors JW, et al., editors. Human Retroviruses and AIDS 1998. Los Alamos: Theoretical Biology and Biophysics Group, Los Alamos National Laboratory; 1998. p. 102–11.
    • 32. Maher SJ, Murray JM. The unrooted set covering connected subgraph problem differentiating between HIV envelope sequences. Eur J Oper Res. 2016;248(2):668–80. http://dx.doi.org/10.1016/j.ejor.2015.07.011.
    • 33. Arthos J, Cicala C, Martinelli E, Macleod K, Van Ryk D, Wei D, et al. HIV-1 envelope protein binds to and signals through integrin alpha4beta7, the gut mucosal homing receptor for peripheral T cells. Nature immunology. 2008;9(3):301–9. Epub 2008/02/12. pmid:18264102
    • 34. Zhou T, Georgiev I, Wu X, Yang ZY, Dai K, Finzi A, et al. Structural basis for broad and potent neutralization of HIV-1 by antibody VRC01. Science. 2010;329(5993):811–7. Epub 2010/07/10. pmid:20616231
    • 35. Kwong PD, Wyatt R, Robinson J, Sweet RW, Sodroski J, Hendrickson WA. Structure of an HIV gp120 envelope glycoprotein in complex with the CD4 receptor and a neutralizing human antibody. Nature. 1998;393(6686):648–59. Epub 1998/06/26. pmid:9641677
    • 36. Gorny MK, Xu JY, Karwowska S, Buchbinder A, Zolla-Pazner S. Repertoire of neutralizing human monoclonal antibodies specific for the V3 domain of HIV-1 gp120. Journal of immunology. 1993;150(2):635–43. Epub 1993/01/15.
    • 37. Kuhlmann AS, Steckbeck JD, Sturgeon TJ, Craigo JK, Montelaro RC. Unique functional properties of conserved arginine residues in the lentivirus lytic peptide domains of the C-terminal tail of HIV-1 gp41. The Journal of biological chemistry. 2014;289(11):7630–40. Epub 2014/02/06. pmid:24497632
    • 38. Korazim O, Sackett K, Shai Y. Functional and structural characterization of HIV-1 gp41 ectodomain regions in phospholipid membranes suggests that the fusion-active conformation is extended. J Mol Biol. 2006;364(5):1103–17. Epub 2006/10/19. pmid:17045292
    • 39. Costin JM, Rausch JM, Garry RF, Wimley WC. Viroporin potential of the lentivirus lytic peptide (LLP) domains of the HIV-1 gp41 protein. Virology journal. 2007;4:123. Epub 2007/11/22. pmid:18028545
    • 40. Kennedy RC, Henkel RD, Pauletti D, Allan JS, Lee TH, Essex M, et al. Antiserum to a synthetic peptide recognizes the HTLV-III envelope glycoprotein. Science. 1986;231(4745):1556–9. Epub 1986/03/28. pmid:3006246
    • 41. Sanders RW, Venturi M, Schiffner L, Kalyanaraman R, Katinger H, Lloyd KO, et al. The mannose-dependent epitope for neutralizing antibody 2G12 on human immunodeficiency virus type 1 glycoprotein gp120. J Virol. 2002;76(14):7293–305. pmid:12072528
    • 42. Cicala C, Martinelli E, McNally JP, Goode DJ, Gopaul R, Hiatt J, et al. The integrin α4β7 forms a complex with cell-surface CD4 and defines a T-cell subset that is highly susceptible to infection by HIV-1. Proceedings of the National Academy of Sciences. 2009;106(49):20877–82.
    • 43. Parrish NF, Wilen CB, Banks LB, Iyer SS, Pfaff JM, Salazar-Gonzalez JF, et al. Transmitted/Founder and Chronic Subtype C HIV-1 Use CD4 and CCR5 Receptors with Equal Efficiency and Are Not Inhibited by Blocking the Integrin α4β7. PLoS Pathog. 2012;8(5):e1002686. pmid:22693444
    • 44. Perez LG, Chen H, Liao H-X, Montefiori DC. Envelope Glycoprotein Binding to the Integrin α4β7 Is Not a General Property of Most HIV-1 Strains. J Virol. 2014;88(18):10767–77. pmid:25008916
    • 45. Byrareddy SN, Kallam B, Arthos J, Cicala C, Nawaz F, Hiatt J, et al. Targeting [alpha]4[beta]7 integrin reduces mucosal transmission of simian immunodeficiency virus and protects gut-associated lymphoid tissue from infection. Nat Med. 2014;20(12):1397–400. http://www.nature.com/nm/journal/v20/n12/abs/nm.3715.html#supplementary-information. pmid:25419708
    • 46. Moore PL, Ranchobe N, Lambson BE, Gray ES, Cave E, Abrahams M-R, et al. Limited Neutralizing Antibody Specificities Drive Neutralization Escape in Early HIV-1 Subtype C Infection. PLoS Pathog. 2009;5(9):e1000598. pmid:19763271
    • 47. Derking R, Ozorowski G, Sliepen K, Yasmeen A, Cupo A, Torres JL, et al. Comprehensive Antigenic Map of a Cleaved Soluble HIV-1 Envelope Trimer. PLoS Pathog. 2015;11(3):e1004767. pmid:25807248
    • 48. Huang J, Kang BH, Pancera M, Lee JH, Tong T, Feng Y, et al. Broad and potent HIV-1 neutralization by a human antibody that binds the gp41-gp120 interface. Nature. 2014;515(7525):138–42. http://www.nature.com/nature/journal/v515/n7525/abs/nature13601.html#supplementary-information. pmid:25186731
    • 49. Falkowska E, Le Khoa M, Ramos A, Doores Katie J, Lee Jeong H, Blattner C, et al. Broadly Neutralizing HIV Antibodies Define a Glycan-Dependent Epitope on the Prefusion Conformation of gp41 on Cleaved Envelope Trimers. Immunity. 2014;40(5):657–68. http://dx.doi.org/10.1016/j.immuni.2014.04.009. pmid:24768347
    • 50. Go EP, Hewawasam G, Liao H-X, Chen H, Ping L-H, Anderson JA, et al. Characterization of Glycosylation Profiles of HIV-1 Transmitted/Founder Envelopes by Mass Spectrometry. J Virol. 2011;85(16):8270–84. pmid:21653661
    • 51. Rong R, Gnanakaran S, Decker JM, Bibollet-Ruche F, Taylor J, Sfakianos JN, et al. Unique Mutational Patterns in the Envelope α2 Amphipathic Helix and Acquisition of Length in gp120 Hypervariable Domains Are Associated with Resistance to Autologous Neutralization of Subtype C Human Immunodeficiency Virus Type 1. J Virol. 2007;81(11):5658–68. pmid:17360739
    • 52. Derdeyn CA, Decker JM, Bibollet-Ruche F, Mokili JL, Muldoon M, Denham SA, et al. Envelope-Constrained Neutralization-Sensitive HIV-1 After Heterosexual Transmission. Science. 2004;303(5666):2019–22. pmid:15044802
    • 53. Moore JP, Willey RL, Lewis GK, Robinson J, Sodroski J. Immunological evidence for interactions between the first, second, and fifth conserved domains of the gp120 surface glycoprotein of human immunodeficiency virus type 1. J Virol. 1994;68(11):6836–47. pmid:7933065
    • 54. Wang WK, Essex M, Lee TH. Single amino acid substitution in constant region 1 or 4 of gp120 causes the phenotype of a human immunodeficiency virus type 1 variant with mutations in hypervariable regions 1 and 2 to revert. J Virol. 1996;70(1):607–11. pmid:8523579
    • 55. Leonard CK, Spellman MW, Riddle L, Harris RJ, Thomas JN, Gregory TJ. Assignment of intrachain disulfide bonds and characterization of potential glycosylation sites of the type 1 recombinant human immunodeficiency virus envelope glycoprotein (gp120) expressed in Chinese hamster ovary cells. J Biol Chem. 1990;265(18):10373–82. Epub 1990/06/25. pmid:2355006
    • 56. Fiebig EW, Wright DJ, Rawal BD, Garrett PE, Schumacher RT, Peddada L, et al. Dynamics of HIV viremia and antibody seroconversion in plasma donors: implications for diagnosis and staging of primary HIV infection. AIDS. 2003;17(13):1871–9. Epub 2003/09/10. pmid:12960819