Humoral Immune Pressure Selects for HIV-1 CXC-chemokine Receptor 4-using Variants

==Highlights==
* CXCR4- as compared to coexisting CCR5-using HIV are less neutralization sensitive.
* CXCR4-using virus can emerge by passaging a CCR5 HIV in the presence of plasma.
* Humoral immune response may lead to the emergence of CXCR4-utilizing variants.

The emergence of HIV-1 variants that use the CXCR4 as opposed to strictly utilizing the CCR5 receptor (termed R5) has been associated with an accelerated progression to AIDS. Potentially, the emergence of CXCR4-using viruses hastens disease or individuals with profound immunosuppression have a greater likelihood of producing CXCR4-utilizing variants. We show that variants that exclusively use the CXCR4 (X4) are neutralization escape variants. Furthermore, CXCR4-using virus can emerge by passaging a R5 variant in the presence of autologous contemporaneous plasma. This suggests that in some individuals, humoral immune pressure against HIV-1 selects for CXCR4-utilizing variants, which accelerates disease progression.

==Abstract==

Although both C-C chemokine receptor 5 (CCR5)- and CXC chemokine receptor 4 (CXCR4)-using HIV-1 strains cause AIDS, the emergence of CXCR4-utilizing variants is associated with an accelerated decline in CD4 + T cells. It remains uncertain if CXCR4-using viruses hasten disease or if these variants only emerge after profound immunological damage. We show that exclusively CXCR4- as compared to cocirculating CCR5-utilizing variants are less sensitive to neutralization by both contemporaneous autologous plasma and plasma pools from individuals that harbor only CCR5-using HIV-1. The CXCR4-utilizing variants, however, do not have a global antigenic change because they remain equivalently susceptible to antibodies that do not target coreceptor binding domains. Studies with envelope V3 loop directed antibodies and chimeric envelopes suggest that the neutralization susceptibility differences are potentially influenced by the V3 loop. In vitro passage of a neutralization sensitive CCR5-using virus in the presence of autologous plasma and activated CD4 + T cells led to the emergence of a CXCR4-utilizing virus in 1 of 3 cases. These results suggest that in some but not necessarily all HIV-1 infected individuals humoral immune pressure against the autologous virus selects for CXCR4-using variants, which potentially accelerates disease progression. Our observations have implications for using antibodies for HIV-1 immune therapy.

==Graphical abstract==


[[Image:draft_García_292415893-1-s2.0-S2352396416301852-fx1.jpg|center|500px|Image 1]]

==Keywords==

HIV ; Envelope ; Neutralization ; Co-receptor ; Evolution ; Resistance ; Tropism

==1. Introduction==

In nearly all cases HIV-1 infection leads to a progressive loss of CD4 + T cells, and both host immunologic factors and viral characteristics influence the rate of progression towards AIDS. Although HIV-1 can enter target cells by either using the C-C chemokine receptor 5 (CCR5) or CXC-chemokine receptor 4 (CXCR4) co-receptor after CD4 receptor engagement ([[#bb0080|Deng et al., 1996]]  and [[#bb0115|Feng et al., 1996]] ), presence of syncytium-inducing (SI) CXCR4-using viruses is associated with accelerated disease progression ([[#bb0190|Koot et al., 1993]] , [[#bb0325|Richman and Bozzette, 1994]]  and [[#bb0065|Connor et al., 1997]] ). In general, viruses that exclusively utilize the CCR5 receptor (termed R5 or non-syncytium-inducing (NSI)) establish initial infection, and these viruses persist throughout the course of disease ([[#bb0065|Connor et al., 1997]] , [[#bb0390|van't Wout et al., 1994]]  and [[#bb0360|Schuitemaker et al., 1992]] ). During the chronic phase of infection, some individuals have a dual-mixed (DM) population consisting of R5 viruses and CXCR4-using variants, which can either only use CXCR4 (termed X4) or be dual-tropic that can engage both co-receptors (termed dual or R5X4) ([[#bb0360|Schuitemaker et al., 1992]]  and [[#bb0065|Connor et al., 1997]] ). Previous studies have not clarified if CXCR4-using viruses accelerate disease progression or if faster CD4 + T cell decline facilitates the emergence of DM virus populations. Deciphering mechanisms for the emergence of CXCR4-utilizing variants will provide insights into this viral factor important for progression towards AIDS.      

It is generally agreed that R5 variants give rise to R5X4 viruses, which can further evolve to X4 HIV-1 ([[#bb0280|Pastore et al., 2004]] ). Specific modifications in the envelope glycoprotein (env) variable loop 3 (V3), which binds the co-receptor, are required for the transition from R5 to X4 virus ([[#bb0160|Huang et al., 2007]] , [[#bb0055|Chesebro et al., 1996]]  and [[#bb0060|Cocchi et al., 1996]] ). The selective pressures which induce the transition from CCR5 to CXCR4 usage during the natural course of infection remain poorly understood. It has been suggested that CXCR4-using viruses potentially emerge because of a reduced number of susceptible CCR5 + target cells or replication differences in various immune cells among the variants that use different co-receptors ([[#bb0315|Regoes and Bonhoeffer, 2005]]  and [[#bb0320|Ribeiro et al., 2006]] ). Indeed, individuals heterozygous for a 32 base pair deletion in CCR5 gene (CCRΔ32), who have lower levels of the CCR5 receptor on their cells, have an increased likelihood of harboring CXCR4-using viruses compared to people with both wild-type alleles ([[#bb0015|Brumme et al., 2005]]  and [[#bb0150|Henrich et al., 2014]] ). Other non-mutually exclusive hypotheses argue that CXCR4-utilizing variants arise because of random mutations, and these new viruses persist because declining adaptive and innate immune responses fail to clear the less fit dual-tropic and/or X4 HIV-1. Previous studies have also suggested that neutralizing antibody (nAb) pressure do not drive co-receptor switching ([[#bb0380|Trkola et al., 1998]] , [[#bb0250|Montefiori et al., 1998]] , [[#bb0040|Cecilia et al., 1998]]  and [[#bb0195|Lacasse et al., 1998]] ). Other investigations have argued that CXCR4-using viruses are highly neutralization susceptible, and the variants that can use the CXCR4 receptor only emerge after collapse of the humoral immune response ([[#bb0020|Bunnik et al., 2007]]  and [[#bb0155|Ho et al., 2007]] ). Within weeks to months after HIV-1 acquisition, infected individuals develop nAbs against their autologous viruses, but the virus changes its env to escape this humoral response ([[#bb0330|Richman et al., 2003]]  and [[#bb0400|Wei et al., 2003]] ). Neutralization resistance is frequently associated with modifications in and around the env variable portions, such as the env V1–V2 and V3 loop ([[#bb0345|Sagar et al., 2006]] , [[#bb0340|Rong et al., 2009]]  and [[#bb0265|Moore et al., 2009]] ). In a cohort of HIV-1 subtype C (HIV-1C) infected individuals, we have previously shown that the env quasispecies diversity was significantly higher among individuals with DM as compared to strictly R5 viruses ([[#bb0215|Lin et al., 2012]] ). In addition, we observed that all HIV-1C X4 as compared to R5 variants had either a 2 amino acid insertion prior to or basic amino acid substitutions in the generally invariant V3 loop crown. Together, these observations led us to hypothesize that difference in host antibody responses drive the greater env diversity observed in individuals with DM as compared to R5 virus populations, and the V3 loop genotypic signatures associated with X4 variants potentially modify neutralization susceptibility. We also postulated that the emergence of specific host antibodies in some but not all individuals force env sequence and structural changes. These env modifications influence co-receptor switching, and the continual presence of the specific host antibodies allows for the persistence of CXCR4-using viruses. In this study, we examined the role of nAbs on HIV-1C co-receptor switching.      

==2. Materials and Methods==

===2.1. Subjects===

HIV-1 envs were isolated from plasma collected from women in the Mashi study previously shown to harbor DM virus ([[#bb0375|Thior et al., 2006]] ).      

===2.2. Ethics Statement===

The use of samples was approved by the institutional review boards at the Harvard School of Public Health, Partners Healthcare Systems, Boston University, and by the Botswana Ministry of Health. Investigators who analyzed the data had no contact with the women enrolled in the Mashi study, and thus all analysis was done in an anonymous fashion.

===2.3. Viruses===

The env expression fragments were generated as previously described using bulk nested RT-PCRs ([[#bb0205|Lin et al., 2010]] ). Briefly, RNA was isolated from 140 μl of individual plasma samples using QIAamp Viral RNA mini kit (Qiagen). For each sample, 1st round RT-PCR was performed in triplicate to limit re-sampling ([[#bb0220|Liu et al., 1996]] ). In order to minimize possible recombination associated with bulk PCR cloning ([[#bb0350|Salazar-Gonzalez et al., 2008]]  and [[#bb0180|Keele et al., 2008]] ), the RT-PCR reactions were performed with long extension times, limited cycle numbers, and use of proofreading polymerase ([[#bb0200|Lahr and Katz, 2009]] , [[#bb0175|Judo et al., 1998]] , [[#bb0245|Meyerhans et al., 1990]]  and [[#bb0110|Fang et al., 1998]] ). Each reaction was done with 4 μl of the extracted viral RNA, SuperScript III RT/Platinum Taq HiFi mix (ThermoFisher Scientific), and an extension time of 7 min for each of the 20 cycles. These 1st round reactions were combined and cloned using TOPO TA cloning (Invitrogen). A CMV promoter was attached to each individually isolated env clone using overlap PCR. All primers and amplification conditions have been described previously ([[#bb0205|Lin et al., 2010]] ). All env sequences have been submitted to Genbank previously (accession numbers KF770248–KF770455, KX060565–KX060578) ([[#bb0215|Lin et al., 2012]] ). Pseudotyped viruses were generated by co-transfecting HEK293T cells with the CMV promoter PCR product attached to the amplified env and an env deficient HIV-1 backbone (pNL4-3R-E-) plasmid expressing luciferase using the Fugene protocol (Roche Molecular Biochemicals).      

===2.4. Generation of Replication Competent Viruses===

Replication competent recombinant viruses were generated with minor modifications from a previously described protocol ([[#bb0045|Chatziandreou et al., 2012]] ). In the majority of cases, we found that full-length HIV-1C envs did not yield replication competent recombinant viruses when incorporated into a NL4-3 (HIV-1B) or Q23-17 (HIV-1A) background ([[#bb0100|Etemad et al., 2014]]  and [[#bb0295|Pena-Cruz et al., 2013]] ). Because the same envs yielded infectious pseudotypes, we assumed that incompatible ''cis''  interactions between the env tail and matrix protein accounted for the lack of infectious virus ( [[#bb0125|Freed and Martin, 1995]] ). To overcome this issue, we generated a plasmid pCMV-Q23-PBS → LTRΔEnvEcto in which the env ecto-domain was replaced with a uracil biosynthesis (''URA3'' ) gene. Briefly, the env ecto-domain in pCMV-Q23-17-PBS → LTR was replaced with URA3 PCR product, which was amplified with the primers AGAAAGAGCAGAAGACAGTGGCAATGA''TTAATTAA'' ACCACCTTTTCAATTCATC and AAGCCTCCTACTATCATTAT CTGATATAATTAAATTGAAGC from pRS316. The underlined portions correspond to primer sequences used to amplify env ectodomains (HXB2 env amino acid numbering 1–692) from subject samples, and the sequences in italics contain a ''Pac'' I restriction site. Different env ecto-domains were shuttled into Q23 plasmid by transforming yeast with pCMV-Q23-PBS → LTRΔEnvEcto linearized with ''Pac'' I and amplified env ecto-domain PCR products from subject samples. Replication competent viruses were generated by co-transfecting HEK293T cells with a plasmid containing the subject env ecto-domain within pCMV-Q23-17-PBS → LTR and another plasmid with Q23–17 sequences from 5′ LTR to early portion of gag, pCMV-Q23-17-LTR → Gag4 ([[#bb0045|Chatziandreou et al., 2012]] ). All virus stocks were collected 48 h post transfection and titered on TZM-bl cells.      

===2.5. Reagents===

The pNL4-3R-E- plasmid, U87 and TZM-bl cell lines, and VRC01 and 10E8 monoclonal antibody (mAb) were obtained through the NIH AIDS Reagent Program. The 293T cells were obtained from American Type Culture Collection. The PGT128 and 447-52D mAbs were provided by International AIDS Vaccine Initiative and Dr. Susan Zolla-Pazner respectively. The CD4-Ig was produced by gene synthesis and cloning into pcDNA3.1 (Invitrogen).

===2.6. Co-receptor Phenotype===

Co-receptor usage was determined by assessing entry capacity in the presence and absence of a co-receptor antagonist on U87 cells expressing CD4 and either CXCR4 or CCR5 as previously described ([[#bb0205|Lin et al., 2010]] ). To examine replication, U87 cells were incubated with virus and polybrene (8 μg/ml) for approximately 16 h. Post-incubation, cells were washed to remove virus. Culture supernatants were examined for p24<sup>gag</sup>  antigen 5 days post-infection. In all cases NL4-3 (X4) and YU-2 (R5) were used as controls.      

===2.7. HIV-1 gp120 Expression and Purification===

The codon-optimized gene encoding the R5 DM268M full-length gp120 was synthesized with a C-terminal Avi-6xHis tag and cloned into the pcDNA3.1(−) expression vector (Life Technologies). The X4 V3 (DM268Y), differed from the R5 V3 (DM268M) by 5 amino acids; these modifications and the CD4 knock-out mutation D368R were introduced by sequential site-directed mutagenesis using the TagMaster Site-Directed Mutagenesis Kit (GM Biosciences) ([[#bb0215|Lin et al., 2012]] ). The gp120 proteins were expressed by transient transfection using the expi293 system (Life Technologies). Six days post transfection, culture supernatants were harvested, filtered, and buffer-exchanged into PBS, and then purified using HisTALON Gravity Column Purification Kit (Clontech Laboratories).      

===2.8. Neutralization Assay===

Infection of TZM-bl cells in the absence and presence of 2-fold serial dilutions was used to estimate the 50% inhibitory concentration (IC<sub>50</sub> ) as previously described ([[#bb0095|Etemad et al., 2009]] ). All reported IC<sub>50</sub> s are mean estimates from a minimum of 3 independent assays. For virus/plasma combinations that failed to achieve more than or equal to 50% inhibition at the highest tested dilution, they were assigned values of 25, which is the midpoint between 0 and highest tested dilution 50. For virus/mAb combinations, a value equal to the highest tested concentration was assigned when 50% inhibition could not be achieved. For the neutralization competition assays, plasma at various dilutions was incubated with media alone or 500 μg/ml of gp120 for 1 h at 37 °C. Virus was added, and wells were incubated for an additional 30 min. Then TZM-bls were added and infections were allowed to proceed for 48 h. Relative light units generated in the absence (positive wells) as compared to the presence of plasma were used to estimate % inhibition at specific plasma dilutions. For estimating inhibition in the presence DM268M-R5V3 or DM268M-X4V3 gp120, the positive wells also contained the respective gp120. In all cases, plasmas showed no neutralization capacity against NL4-3 pseudovirions with a vesicular stomatitis virus (VSV) envelope G protein (data not shown). Thus, observed neutralization was specific against the HIV-1 envelope.      

===2.9. ELISA===

Plates were coated with gp120 antigen in PBS (pH 7.4) at 2 μg/ml overnight at 4 °C, blocked with 1% BSA in PBS for 1 h at 37 °C, and then incubated with serially diluted plasma or antibodies in blocking buffer for 1 h at 37 °C. Horseradish peroxidase (HRP)-conjugated goat anti-human IgG Fc antibody (Jackson ImmunoResearch Laboratories Inc.) was added in blocking buffer at 1:10,000 for 1 h at 37 °C. Plates were washed between each step with 0.1% Tween 20 in PBS and developed with 3,3′,5,5′-tetramethylbenzidine (TMB) (Life Technologies) for 20 min at room temperature. 1 M sulfuric acid was added to terminate the reaction and the plates were read at 450 nm.

===2.10. Virus Passage===

Viruses were passaged on CD4 + T cells in the presence of increasing concentrations of heat inactivated autologous plasma. Peripheral blood mononuclear cells (PBMCs) were isolated from HIV-1 negative blood donors' buffy packs using Ficoll-Hypaque density centrifugation method. CD4 + T cells were isolated from the PBMCs using magnetic bead isolation (Stem Cell Technologies). Around 2 × 10<sup>6</sup>  phytohemagglutinin (PHA) and interleukin-2 (IL-2) activated CD4 + T cells were exposed to infectious virus in the presence of 20 μg/ml diethylaminoethyl (DEAE) dextran. Plasma was added after infections were established, generally two days after infection. Cultures were monitored for viral replication twice a week using a p24<sup>gag</sup>  antigen assay (Perkin Elmer). The cultures were fed twice weekly with activated CD4 + T cells, and plasma concentration was doubled if virus showed continued growth.      

===2.11. Statistical Analysis===

In the majority of cases, individual env variants IC<sub>50</sub> s was a mean of 3 independent assays. Comparisons within an individual were done using Wilcoxon rank sum test. For comparisons across multiple individuals, the median IC<sub>50</sub>  for strains with specific co-receptor usage was determined for each subject. Inter-subject comparison were done using Kruskal Wallis test, and pair-wise inter-subject comparisons were done using the Wilcoxon rank sum test. The matched pairs Wilcoxon rank sum test was used when all individuals had variants with the co-receptor usage being compared. All statistical analysis were done using Stata version 8.0 (Stata Corporation). All p-values are based on two sided tests.      

==3. Results==

===3.1. Subjects===

Previously, we determined the co-receptor usage in 148 HIV-1C infected women in Botswana, and we found that 22 of them harbored DM populations while 126 had strictly R5 virus ([[#bb0210|Lin et al., 2011]] ). Sequence characteristics and co-receptor use of individual envs were examined from 10 of the DM samples that showed relatively high CXCR4 usage by the phenotypic assay ([[#bb0215|Lin et al., 2012]] ). Examination of individuals that harbored X4 variants were prioritized because X4 as compared to co-existing R5 and R5X4 envs were demonstrated to harbor unique genotypic features ([[#bb0215|Lin et al., 2012]] ). In addition, individual envs were not examined from DM samples that demonstrated low amount of CXCR4 usage in the co-receptor assay because it would likely require more extensive cloning to isolate CXCR4-using variants. Only 5 from the original 148 samples fulfilled these criteria. We chose to examine an additional 2 samples with only co-circulating R5 and R5X4 variants. We assessed neutralization susceptibility of a median of 11 co-circulating envs (range: 9–18) from these 7 different women ([[#t0005|Table 1]] ). The co-existing env populations consisted of X4 and R5 variants in 3 (DM8, DM178, and DM268), R5X4 and R5 variants in 2 (DM146 and DM269), R5X4 and X4 variants in 1 (DM173), and all three co-receptor phenotypes in 1 (DM159) of the individuals. An additional 6 DM8 clones failed to yield other X4 envs. Intra-individual env variants with the different co-receptor usage were phylogenetically related, which ruled out contamination and super-infection (Fig. S1). Even though all 88 primary envs examined in this study were generated using bulk-cloning PCR as described previously ([[#bb0205|Lin et al., 2010]]  and [[#bb0215|Lin et al., 2012]] ), all variants had unique sequences. Thus, we avoided re-sampling the same variant multiple times, which is a bias often associated with the bulk PCR cloning technique used to isolate the envs in this study ([[#bb0350|Salazar-Gonzalez et al., 2008]]  and [[#bb0180|Keele et al., 2008]] ).

<span id='t0005'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Table 1.

Subject demographics and patient characteristics.

|-

! rowspan="2" | Subject
! rowspan="2" | Plasma HIV (copies/ml)                                                    
! rowspan="2" | CD4 count (cells/μl)                                                    
! rowspan="2" | Previous ART[[#tf0005|<sup>a</sup>]]
! colspan="3" | Env variants
|-

! R5
! Dual
! X4
|-

! DM8
| 522,000
| 43
| N
| 10
| 0
| 1
|-

! DM146
| 653,000
| 172
| N
| 10
| 8
| 0
|-

! DM159
| 154,000
| 116
| N
| 5
| 2
| 3
|-

! DM173
| 1290
| 316
| N
| 0
| 6
| 3
|-

! DM178
| 230,000
| 21
| Y
| 8
| 0
| 6
|-

! DM268
| 333,000
| 375
| Y
| 5
| 0
| 9
|-

! DM269
| 10,600
| 319
| Y
| 9
| 3
| 0
|}

a. Women were sampled either prior to or a minimum of 6 months after the discontinuation of antiretroviral therapy (ART) because of virological failure.

===3.2. X4 Are Less Neutralization Susceptible to Contemporaneous Plasma===

To assess neutralization sensitivity differences among 88 co-circulating variants with different observed co-receptor usage, we first examined susceptibility to contemporaneous autologous plasma. In 5 individuals, X4 strains were less sensitive to contemporaneous autologous plasma as compared to co-existing R5 (DM8, DM178, DM268, DM159) or R5X4 (DM173) viruses ([[#f0005|Fig. 1]] a–e). Indeed, in 3 individuals (DM8, DM178, and DM268), the X4 envs were not neutralized above 50% at the highest tested plasma dilution (1:50) ([[#f0005|Fig. 1]] a–c). These virus-plasma combinations were assigned a value of 25, which is the mid-point value between 0 and 50, for statistical comparisons. In 3 individuals (DM178, DM268, and DM173) differences were statistically significant (p < 0.05, Wilcoxon rank sum test), while in the other 2 subjects (DM8 and DM159), limited number of envs with different co-receptor usage potentially limited the power to detect significant differences. In 2 individuals (DM146 and DM269) co-circulating R5 and dual-tropic envs showed similar neutralization sensitivity (p > 0.05, Wilcoxon rank sum test) ([[#f0005|Fig. 1]] f–g). Aggregate comparisons were done by estimating a median IC<sub>50</sub>  for each individuals envs that used different co-receptor. In aggregate, neutralization IC<sub>50</sub> s to autologous contemporaneous plasma was significantly different among R5, R5X4, and X4 envs (p = 0.03, Kruskal Wallis test). The median neutralization susceptibility to contemporaneous autologous plasma of the X4 viruses (median IC<sub>50</sub> : 25.0, range: 25.0–140.7) was around 5 fold lower compared to R5 (median IC<sub>50</sub> : 136.0, range: 94.2–222.7, p = 0.03, Wilcoxon rank sum test) and R5X4 viruses (median IC<sub>50</sub> : 156.5, range: 135.2–160.8, p = 0.03, Wilcoxon rank sum test) ([[#f0005|Fig. 1]] h).

<span id='f0005'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_292415893-1-s2.0-S2352396416301852-gr1.jpg|center|px|Fig. 1]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 1.

The X4 as compared to co-circulating R5 and R5X4 envs are more neutralization resistant to autologous contemporaneous plasma. (a–g) Figures show neutralization IC<sub>50</sub> s for R5 (circles), dual-tropic (triangles), and X4 (boxes) viruses against autologous contemporaneous plasma. Symbols represent one of the 88 independent envs from the 7 subjects labeled above each graph. Y-axis shows the IC<sub>50</sub>  against autologous contemporaneous plasma as reciprocal dilution. Figure h shows the aggregate comparisons in all subjects, and each symbol shows the median IC<sub>50</sub> s within a subject for envs with the defined co-receptor phenotype. In every case, the long line represents the median of the dot plots, and the whiskers show the interquartile range. Star indicates statistically significant difference (p < 0.05). The dashed line shows the highest tested plasma dilution (1:50). Symbols below the dotted lines are envs that were not neutralized by at least 50% at the highest tested plasma dilution. For statistical comparisons, these envs were assigned a midpoint value between 0 and 50.                  

</span>
|}

In preliminary data, we also found an HIV-1B infected individual also harbored X4 as compared to CCR5-using variants that were significantly less sensitive to neutralization by autologous contemporaneous plasma (Fig. S2). This suggests our observations are not specific to HIV-1C infections only. Thus, in every individual with isolated X4 viruses (n = 6), the X4 as compared to co-existing R5 and/or R5X4 variants had greater neutralization resistance to autologous contemporaneous plasma ([[#f0005|Figs. 1]]  and S2). Previous studies have noted similar IC<sub>50</sub>  fold changes among longitudinally collected viruses to sequential autologous plasma ([[#bb0330|Richman et al., 2003]] ). In addition, the observed neutralization response against contemporaneous viruses in the individuals examined in this study is not necessarily unique because recent work shows that plasma and isolated antibodies neutralize some autologous contemporaneous viruses to a similar extent, which subsequently provides pressure for envelope evolution ([[#bb0255|Moody et al., 2015]] ). In aggregate, this suggests that the × X4 as compared to co-existing CCR5-using viruses have both statistical and biological meaningful differences in neutralization susceptibility. Thus, X4 variants are potentially neutralization escape variants to the antibodies present in the contemporaneous autologous plasma.      

===3.3. X4 Are More Resistant to Heterologous Pooled Plasma===

We next hypothesized that X4 envs are escape variants from antibodies often present in individuals that harbor only CCR5-using viruses. To examine this possibility, we compared sensitivity of X4 and R5 strains to heterologous pooled R5 plasma, which was generated by pooling plasma from 5 women documented to harbor R5 viruses only ([[#bb0215|Lin et al., 2012]] ). We examined 47 envs from the 4 individuals (DM8, DM159, DM178, and DM268) in which the median IC<sub>50</sub>  against autologous contemporaneous plasma was lower for the X4 as compared to co-circulating R5 strains. We did not examine variants from DM146, DM173, and DM269 because in these individuals we did not isolate co-existing R5 and X4 strains. In all 4 individuals, the X4 viruses demonstrated decreased sensitivity to this pooled R5 plasma compared to the co-existing CCR5 utilizing viruses ([[#f0010|Fig. 2]] a–d). Significant differences (p < 0.05, Wilcoxon rank sum test) were observed in 2 cases (DM178 and DM268) but not in the other 2 individuals (DM8 and DM159) with a limited number of envs with different co-receptor usage. In aggregate, X4 (median IC<sub>50</sub> : 161.7, range: 122.8–179.6) as compared to R5 (median IC<sub>50</sub> : 425.3, range: 174.5–684.8) viruses showed a statistical trend of around 3 fold greater neutralization resistance to the pooled R5 plasma (p = 0.1, paired Wilcoxon rank sum test) ([[#f0010|Fig. 2]] e). This suggests that X4 as compared to co-existing R5 viruses are less neutralization susceptible to antibodies commonly circulating in individuals that harbor only R5 variants.

<span id='f0010'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_292415893-1-s2.0-S2352396416301852-gr2.jpg|center|px|Fig. 2]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 2.

The X4 as compared to co-circulating R5 are more neutralization resistant to heterologous pooled R5 plasma. (a–d) Graphs show neutralization IC<sub>50</sub> s (reciprocal plasma dilution) for R5 (circles) and X4 (boxes) viruses against heterologous pooled R5 plasma for 4 subjects. Each symbol represents one of the 47 independent envs from the 4 individuals shown above the graph. Figure E shows the aggregate comparisons in the 4 subjects, and each symbol shows the median IC<sub>50</sub> s for the R5 or X4 envs within an individual. In every case, the long line represents the median of the dot plots with the whiskers showing the interquartile range. Star indicates statistically significant difference (p < 0.05). The dashed line shows the highest tested plasma dilution (1:50).                  

</span>
|}

===3.4. X4 Variants Remain Sensitive to Antibodies Targeting env Non-co-receptor Binding Sites===

It is possible that mutations in the X4 viruses compared to the R5 variants leads to a “global” antigenic change, similar to some previously described strains ([[#bb0075|del Prete et al., 2009]]  and [[#bb0335|Roederer et al., 2013]] ). If the X4 as compared to the R5 variants have an overall antigenic difference, they may have different neutralization susceptibility to antibodies targeting env regions not implicated in co-receptor binding. To explore this possibility, we assessed sensitivity to two broadly neutralizing antibodies (bnAbs), VRC01, which targets the CD4 binding site (CD4bs) ([[#bb0410|Wu et al., 2010]]  and [[#bb0415|Zhou et al., 2010]] ), and 10E8, which binds the membrane proximal external region (MPER) ([[#bb0165|Huang et al., 2012]] ). In individual subjects, the X4 as compared to R5 viruses demonstrated similar sensitivity to neutralization (p > 0.05, Wilcoxon rank sum test) to VRC01 ([[#f0015|Fig. 3]] a–d) and 10E8 ([[#f0015|Fig. 3]] e–h). Although in some individuals the highest tested bnAb concentration failed to yield 50% inhibition, this lack of neutralization was observed in both the X4 and co-existing R5 variants. In these cases, virus-antibody combination was assigned a value equivalent to the highest tested concentration for statistical comparisons. In aggregate, the median VRC01 and 10E8 IC<sub>50</sub>  of an individuals X4 and co-existing R5 variants were not significantly different (p > 0.05, paired Wilcoxon rank sum test). Thus, the co-existing X4 and R5 variants demonstrated equivalent sensitivity to antibodies that target env regions not involved in co-receptor binding, which suggests that changes that confer co-receptor switching do not lead to overall structural change in CD4 binding site or domains close to the membrane. Comparing sensitivity to other bnAbs that target unique sites, such as the V1–V2 loop or the interface between the trans-membrane domain and surface unit, will need to be examined to assess structural changes in other env regions.

<span id='f0015'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_292415893-1-s2.0-S2352396416301852-gr3.jpg|center|px|Fig. 3]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 3.

The X4 variants do not have a “global” antigenic change. Neutralization IC<sub>50</sub> s for R5 (circles), dual-tropic (triangles), and X4 (boxes) viruses to VRC01 (a–d) and 10E8 (e–h). Y-axis shows the bnAb concentration in μg/ml. Dotted line shows the highest tested antibody concentration. For statistical comparisons, virus antibody combination that failed to yield > 50% neutralization at the highest tested antibody concentration were assigned an IC<sub>50</sub>  value equivalent to the highest tested antibody concentration. In every case, the long line represents the median of the dot plots with the whiskers showing the interquartile range.                  

</span>
|}

===3.5. Infected Hosts Develop Antibodies Against X4 Variants===

A study using a simian-human immunodeficiency virus with a subtype B env (SHIV-1B) has suggested that CXCR4-using variants are exquisitely neutralization sensitive, and these viruses potentially arise in setting where there are limited nAb responses in the infected host ([[#bb0155|Ho et al., 2007]] ). If CXCR4 viruses emerge after immunological collapse then infected hosts will potentially not be able to develop antibody responses against X4 variants over time. We assessed sensitivity to an autologous plasma sample collected around 6 months after the time point at which the DM8 X4 and R5 variants were isolated to assess this notion. The DM8 X4 compared to the co-circulating R5 variants demonstrated similar sensitivity to plasma collected 6 month after env isolation (Fig. S3a). Importantly, all R5 and the one X4 env were more susceptible to the later compared to the contemporaneous plasma sample (p = 0.008, paired Wilcoxon rank sum test) (Fig. S3b). Autologous plasma samples from later time-points were not available from other individuals to further examine if infected hosts could generate nAbs against previously circulating envs after the emergence of DM virus populations. These observations suggest that subject DM8 with a DM virus population continues to develop humoral responses over time against previously co-circulating R5 and X4 variants.      

===3.6. Neutralization Susceptibility Differences Are Influenced by env V3 Loop Changes===

Our results suggest that X4 variants are relatively insensitive to specific antibodies present in autologous plasma and heterologous antibodies in individuals with only R5 circulating variants. We hypothesized that these antibodies potentially target the env V3 loop because the V3 loop is the primary determinant for co-receptor usage ([[#bb0025|Cann et al., 1992]] ). We compared sensitivity of co-existing X4 and R5 variants to the bnAb, PGT128, because V3 loop stem sequence partially influences susceptibility to this bnAb ([[#bb0395|Walker et al., 2011]]  and [[#bb0290|Pejchal et al., 2011]] ). In DM268, X4 viruses (median IC<sub>50</sub> : 0.14 μg/ml, range: 0.04–0.22) were > 100 fold less neutralization susceptible to PGT128 compared to co-circulating R5 (median IC<sub>50</sub> : 0.001 μg/ml, range: 0.001–0.02) variants (p = 0.003, Wilcoxon rank sum test) ([[#f0020|Fig. 4]] a). In contrast, neutralization sensitivity to PGT128 for the one DM8 X4 variant (IC<sub>50</sub> : 0.41 μg/ml) was within the range for the co-circulating R5 variants (median IC<sub>50</sub> : 1.8 μg/ml, range: 0.37–3.67) ([[#f0020|Fig. 4]] b). None of the DM159 and DM178 X4 and R5 variants were inhibited > 50% at the highest tested PGT128 concentration (5 μg/ml), and thus an IC<sub>50</sub> s could not be calculated for comparisons.

<span id='f0020'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_292415893-1-s2.0-S2352396416301852-gr4.jpg|center|px|Fig. 4]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 4.

Neutralization susceptibility differences among X4 and R5 viruses are influenced by env V3 loop. Figures a–b show mean IC<sub>50</sub> s (μg/ml) to PGT128 among co-circulating variants with different co-receptor usage in DM268 (a) and DM8 (b). Figures c–g show IC<sub>50</sub>  against pooled R5 plasma (reciprocal dilution) (c–e) and PGT128 (μg/ml) (f–g) for the original R5 (white), original X4 (black), env with a R5 V3 loop in a X4 background (gray) or X4 V3 loop in a R5 background (checkers). In each graph, the subject is denoted above. Label below the x-axis denotes the env clone followed by the observed co-receptor usage in parenthesis. Star denotes significant p-value for the comparison (p < 0.05, Wilcoxon rank sum test). DM268 chimeric env MYM was not infectious ([[#bb0215|Lin et al., 2012]] ).                  

</span>
|}

To further examine the role of V3 loop in the observed neutralization differences among co-circulating R5 and X4 variants, we generated V3 loop chimeric envs. In 3 subjects (DM8, DM159, and DM268) with documented differences in neutralization sensitivity among co-existing X4 and R5 variants, we exchanged V3 loops between a R5 and X4 variant. The chimeric DM178 envs with V3 loop exchanges did not yield infectious virus. Swapping the V3 loops resulted in co-receptor changes in every case other than in DM8 ([[#f0020|Fig. 4]] c–e) ([[#bb0215|Lin et al., 2012]] ). The chimeric envs with X4 V3s were less sensitive to pooled R5 heterologous plasma compared to the envs containing the R5 V3 loops in every case. The differences, however, were relatively small especially compared to the susceptibility differences among the original envs. Furthermore, in both DM268 and DM8, introduction of an X4 V3 loop within the R5 background yielded a virus with greater resistance to PGT128 compared to the original env or the recombinant with a R5 V3 loop ([[#f0020|Fig. 4]] f–g). In aggregate, these results suggest that the V3 loop potentially influences the observed neutralization differences among co-circulating R5 and X4 variants.      

ELISAs were used to further examine V3 directed antibody binding variation to co-existing DM268 X4 and R5 variants. We only assessed DM268 envs because of the observed sensitivity difference to bnAb, PGT128. We generated gp120 proteins based on the predicted amino acids in a DM268 R5 variant (clone M). One gp120 contained an R5 V3, DM268-M-R5V3, and the other harbored an X4 V3 from clone Y, DM268-M-X4V3. Both contained a D368R mutation to avoid CD4 interference ([[#bb0275|Olshevsky et al., 1990]] ). As expected, the CD4 binding site bnAb, VRC01, showed no binding to either gp120 because of the D368R mutation ([[#f0025|Fig. 5]] a). On the other hand, both gp120s showed relatively equivalent weak binding to CD4-Ig ([[#f0025|Fig. 5]] b), which is similar to what has been shown with other gp120s containing the D368R mutation ([[#bb0170|Jia et al., 2016]] ). The DM268-M-R5V3 gp120 as compared to DM268-M-X4V3, however, showed markedly greater binding to PGT128 ([[#f0025|Fig. 5]] c), V3 loop crown antibody, 447-52D ([[#f0025|Fig. 5]] d), and contemporaneous autologous serum ([[#f0025|Fig. 5]] e) ([[#bb0135|Gorny et al., 1992]] ). Although, 447-52D demonstrated ELISA binding, this mAb was unable to neutralize any DM268M variant at the highest tested concentration (5 μg/ml) (data not shown). These ELISA results suggest that the DM268 plasma harbors anti V3 loop antibodies that preferentially bind R5 but not X4 envs.

<span id='f0025'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_292415893-1-s2.0-S2352396416301852-gr5.jpg|center|px|Fig. 5]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 5.

Co-circulating R5 and X4 viruses have different reactivity to anti-env V3 loop antibodies and autologous plasma. Figure a thru e show ELISA binding profiles of the DM268_D368R gp120 with either a R5 V3 loop (circle) or X4 V3 loop (square) to VRC01 (a), CD4-Ig (b), PGT128 (c), 447-52D (d), and autologous contemporary plasma (e). The y-axis shows the OD and the x-axis shows the antibody concentration or plasma dilution.

</span>
|}

Neutralization competition experiments were carried out to assess if the presence anti V3 loop antibodies in the autologous plasma is responsible for the differential neutralization of DM268 R5 and X4 variants. If anti V3 loop antibodies that preferentially bind R5 but not X4 envs were primarily responsible for neutralization susceptibility differences among DM268 X4 and R5 variants, it would be expected that gp120s with R5 but not X4 V3 would prevent autologous contemporaneous plasma from neutralizing the virus. In the presence of DM268-M-R5V3 gp120 (reciprocal plasma dilution IC<sub>50</sub> : 132.1) neutralization of the DM268M virus required greater amount of plasma compared to in the presence of DM268-M-X4V3 gp120 (reciprocal plasma dilution IC<sub>50</sub> : 168.3) or the absence of gp120 (reciprocal plasma dilution IC<sub>50</sub> : 200.0) ([[#f0030|Fig. 6]] ). These estimated IC<sub>50</sub>  differences, however, were not statistically significantly different. This suggests that even though the plasma contains anti V3 loop antibodies that preferentially bind R5 but not X4 gp120s, these antibodies had a small impact in neutralizing a DM268 R5 strain.

<span id='f0030'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_292415893-1-s2.0-S2352396416301852-gr6.jpg|center|358px|Fig. 6]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 6.

gp120s with different env V3 loops fails to prevent neutralization. Percent inhibition observed in the absence or presence of DM268_D368R gp120 with R5 V3 or X4 V3. The y-axis shows the percent inhibition and the x-axis shows the reciprocal plasma dilution.

</span>
|}

===3.7. Emergence of CXCR4-using Virus With Serial Passage in Plasma===

To understand if CXCR4-using variants emerge from R5 viruses as a consequence of nAb pressure, we independently passaged a neutralization sensitive R5 variant from 3 different subjects (DM8, DM178, and DM268) in activated CD4 + T cells in the absence or presence of heat inactivated contemporaneous autologous plasma ([[#f0035|Fig. 7]] a–c). While all the previous experiments used pseudovirions with full-length envs, replication competent recombinant viruses were used in these passages. These replication competent recombinant viruses incorporated an individuals neutralization sensitive R5 env ecto-domain rather than the full-length env. Incorporation of full-length envs within a heterologous HIV-1 backbone did not yield infectious virus. In all 3 cases, virus passaged in the presence as compared to the absence of autologous plasma yielded virus supernatants that were neutralization resistant ([[#f0035|Fig. 7]] d–f). Passage of DM268M R5 variant in the presence but not the absence of serially increasing amounts of contemporaneous autologous plasma led to the emergence of a CXCR4-using virus ([[#f0040|Fig. 8]] a–b). The day 75 DM268 virus supernatant passaged in the presence as opposed to the absence of autologous plasma replicated in U87–CXCR4 cells that lacked the CCR5 receptor ([[#f0040|Fig. 8]] a). Replication decreased in the presence of CXCR4 inhibitor, AMD3100, suggesting that cell entry was mediated by engaging the CXCR4 co-receptor. CXCR4 usage was highly inefficient because day 75 plasma passaged DM268 virus supernatant yielded nearly 10 fold less HIV-1 antigen (p24<sup>gag</sup> ) in the U87–CXCR4 as compared to U87–CCR5 cells ([[#f0040|Fig. 8]] b). Similar results were observed with day 68 DM268 virus supernatant passaged in the presence of autologous plasma (data not shown). In contrast, all other virus supernatants passaged either in the absence or presence of plasma infected U87–CCR5 cells but not U87–CXCR4 cells.

<span id='f0035'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_292415893-1-s2.0-S2352396416301852-gr7.jpg|center|px|Fig. 7]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 7.

Passage in the presence of autologous contemporaneous plasma leads to neutralization resistant virus. Figure a thru c show passage of DM8 clone S (a), DM178 clone A (b), and DM268 clone M (c) in the absence (circle) or presence (square) of autologous contemporaneous plasma. In each graph, the x-axis shows the days post-infection. The left y-axis shows the estimated p24 concentration in the culture and the right y-axis shows contemporaneous autologous plasma dilution in the culture (triangles). Figures d thru f show neutralization susceptibility of passaged virus stock to contemporaneous autologous plasma. The y-axis shows the percent inhibition and the x-axis shows the reciprocal of the autologous plasma dilution. Any virus/plasma combination that did not show a positive inhibition value was assigned a % inhibition of 10% to improve visualization on the graph. Virus stocks passaged in the presence (square) or absence of (circle) of autologous plasma were from day 34 (DM8) (d), day 92 (DM178) (e), and day 75 (DM268) (f) post-infection.

</span>
|}

<span id='f0040'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_García_292415893-1-s2.0-S2352396416301852-gr8.jpg|center|px|Fig. 8]]


|-
| <span style="text-align: center; font-size: 75%;">

Fig. 8.

A CXCR4-using virus emerged with passage in the presence as opposed to the absence of autologous plasma in DM268. Replication in U87-CD4-CXCR4 (a) and U87-CD4-CCR5 (b) cells for day 34 (DM8), day 92 (DM178), and day 75 (DM268) virus stocks. The table below indicates whether the virus stock was from passage in the presence or absence of autologous plasma. The table also shows replication observed in the presence (white bars) or absence (black bars) of specific CXCR4 (AMD3100) and CCR5 (TAK779) inhibitors. Graphs also show replication of control viruses NL43 (X4) and YU-2(R5). Each bar shows standard error from 2 independent measurements. Dotted line shows the limit of detection (25 pg/ml). Cultures in which the p24 antigen level was below the limit of detection were assigned a value midpoint between 0 and 25 to improve visualization on the graph.

</span>
|}

We hypothesized that the day 75 DM268 virus stock passaged in the presence of plasma (DM268 + plasma-d75) contained a possible mixture of R5, R5X4, and X4 viruses. To efficiently isolate the CCR5- versus CXCR4-using viruses in DM268 + plasma-d75, we examined viruses in the U87–CXCR4 and U87–CCR5 cell culture supernatants. We reasoned that the U87–CXCR4 and the U87–CCR5 cell culture supernatants should only contain CXCR4- and CCR5-utilizing viruses respectively. Interestingly, the culture supernatants from the U87–CCR5 and U87–CXCR4 cells exposed to DM268 + plasma-d75 showed identical env bulk sequences. This suggests that passage in the presence of autologous plasma led to the emergence of a neutralization resistant R5X4 variant with relatively inefficient CXCR4 usage. Surprisingly, this R5X4 variant had an amino acid change at position 88 (HXB2 numbering), which introduced a predicted asparagine (N) linked glycosylation site (PNGS) (Fig. S4). In contrast, passage in the absence of plasma generated virus stocks that contained a mixture of amino acids at 2 positions but was otherwise invariant from the starting input virus. Interestingly, all other R5 and X4 DM268 variants (n = 17) previously isolated from this individual harbored the PNGS at position 88. Site directed mutagenesis could not be used to assess if this modification would change co-receptor usage in other strains because all other isolated R5 envs contained the position 88 PNGS. The effect of the position 88 change on co-receptor usage, however, is context dependent because none of the other previously isolated DM268 env variants with a position 88 PNGS were able to use both the CCR5 and CXCR4 co-receptor for cell entry. These observations imply that passage in the presence of autologous plasma led to the emergence of a PNGS at position 88, which subsequently conferred CXCR4 usage.

==4. Discussion==

This study demonstrates that X4 variants are more neutralization resistant compared to R5 and R5X4 viruses and provides a potential biological mechanism for HIV-1 co-receptor switching. In multiple subjects, we show that co-existing X4 as compared to R5 variants are more resistant to antibodies present in autologous contemporaneous and heterologous pooled plasma from individuals with only R5 strains. This decreased neutralization sensitivity did not result from changes that led to a global antigenic change. Thus, this argues that X4 viruses are more resistant to certain antibodies, likely against the env V3 loop, generated by the infected individual. Furthermore, we show that passage of neutralization sensitive CCR5-using variants in the presence of increasing concentrations of autologous contemporaneous plasma was sufficient to force the evolution of a CXCR4-using virus from an R5 strain in 1 of 3 subjects. In aggregate, our studies suggest that humoral immune pressure can drive the emergence of CXCR4-utilizing viruses in some individuals.

Studies have suggested that dual-tropic strains evolve from R5 variants, and X4 viruses arise from the R5X4 strains ([[#bb0280|Pastore et al., 2004]] ). It is unlikely that the emergence of a dual-tropic virus during the passage experiments occurred by chance because CXCR4-using variants were not detected in the absence of plasma and changes required for co-receptor switching are generally associated with diminished replication fitness ([[#bb0280|Pastore et al., 2004]] ). The DM268 passage results were surprising because we observed no significant neutralization susceptibility difference to contemporaneous autologous plasma among co-existing R5 and dual-tropic strains ([[#f0005|Fig. 1]] ). Although previous studies have demonstrated that diverse glycosylation changes can render viruses neutralization resistant ([[#bb0345|Sagar et al., 2006]]  and [[#bb0400|Wei et al., 2003]] ), it remains unclear how the addition of PNGS at position 88 leads to low level CXCR4 usage. It has been shown that co-circulating R5 and R5X4 variants often have the same predicted V3 loop sequence ([[#bb0160|Huang et al., 2007]]  and [[#bb0215|Lin et al., 2012]] ), and thus it possible that modifications in other env segments lead to a structural change that can influence co-receptor usage. Interestingly, DM viruses did not emerge in the other 2 cases (DM8 and DM178) even though their X4 as compared to R5 envs were more neutralization resistant to autologous contemporaneous plasma. Our observations suggest that the passaged DM8 and DM178 env variants accommodated changes required for neutralization resistance to autologous plasma without co-receptor switching. Other R5 variants from these individuals, not examined in the passaging experiments, may acquire the ability to use the CXCR4 receptor as they become resistant to autologous antibodies.      

The earliest studies implied that lab-adapted SI, presumably CXCR4-using, strains were more neutralization sensitive compared to NSI–CCR5-utilizing variants ([[#bb0240|Mascola et al., 1996]] , [[#bb0250|Montefiori et al., 1998]] , [[#bb0090|D'Souza et al., 1997]]  and [[#bb0355|Sawyer et al., 1994]] ). Subsequent studies showed that compared to CCR5-utilizing strains primary as opposed to lab-adapted CXCR4-using envs displayed a relatively similar range of neutralization susceptibilities to some mAbs and unrelated heterologous plasma ([[#bb0380|Trkola et al., 1998]] , [[#bb0250|Montefiori et al., 1998]] , [[#bb0040|Cecilia et al., 1998]]  and [[#bb0195|Lacasse et al., 1998]] ). This led to the notion that nAbs are unlikely to be the major selection pressure that drives co-receptor switching. Our results also show that X4 and co-existing R5 variants have equivalent sensitivity to antibodies targeting env domains not involved in co-receptor binding. More recent investigations with SHIV-1B and HIV-1B infected subjects argue that relatively neutralization sensitive CXCR4-utilizing variants emerge later in disease after HIV-1 induced immunological damage has led to a compromised humoral immune response ([[#bb0020|Bunnik et al., 2007]]  and [[#bb0155|Ho et al., 2007]] ). We, however, observed no case in which X4 or R5X4 as compared to co-circulating R5 envs were significantly more neutralization susceptible to autologous contemporaneous plasma, even though we sampled more envs from each subject and examined a greater number of individuals. Furthermore, the previous human study examined PBMC passaged variants ([[#bb0020|Bunnik et al., 2007]] ), which may have led to the different results because passaging virus in PBMCs significantly impacts its neutralization properties in a unique manner for each virus antibody combination ([[#bb0105|Etemad et al., 2015]]  and [[#bb0310|Provine et al., 2012]] ).      

Our results also differ from the previous human study ([[#bb0020|Bunnik et al., 2007]] ) potentially because we examined HIV-1C versus HIV-1B infected individuals. Our preliminary results, however, suggest this is not the case (Fig. S2). While HIV-1C constitutes the majority of worldwide infections, emergence of CXCR4 usage occurs at a lower frequency compared to HIV-1B ([[#bb0210|Lin et al., 2011]]  and [[#bb0005|Ataher et al., 2012]] ). Reasons for the difference in incidence remain uncertain. Although it has been reported that CCR5 and CXCR4-using viruses preferentially infect memory and naïve T cells respectively ([[#bb0070|Davenport et al., 2002]]  and [[#bb0010|Blaak et al., 2000]] ), there is no evidence that this preference differs among HIV-1B versus HIV-1C variants ([[#bb0030|Cashin et al., 2014]]  and [[#bb0320|Ribeiro et al., 2006]] ). While, it is known that the CCR5Δ32 allele frequency is higher among Europeans than Africans ([[#bb0130|Galvani and Novembre, 2005]] ), the density of the CCR5 receptor on CD4 + T cells is not significantly different among people with European ancestry who constitute the majority of HIV-1B infections versus those with African heritage who are predominantly infected with HIV-1C ([[#bb0305|Picton et al., 2012]] ). Thus, the higher incidence of DM virus population in HIV-1B as compared to HIV-1C infected individuals is unlikely because of lower availability of appropriate CCR5 expressing target cells. It has been suggested that HIV-1B and HIV-1C envs have structural differences, especially in the V3 loop ([[#bb0285|Patel et al., 2008]] ). Potentially in response to these env structural differences, HIV-1C infected individuals develop antibodies targeting the C3-V4 region while those with HIV-1B generate responses against the V1-V2 domain and the base of the V3 loop ([[#bb0340|Rong et al., 2009]] , [[#bb0260|Moore et al., 2008]]  and [[#bb0370|Tang et al., 2011]] ). Interestingly, the main determinants for co-receptor switching are located in the V1 to V3 domains ([[#bb0270|O'Brien et al., 1990]] , [[#bb0050|Chesebro et al., 1991]] , [[#bb0365|Shioda et al., 1992]] , [[#bb0140|Groenink et al., 1993]]  and [[#bb0120|Fouchier et al., 1995]] ). Greater prevalence of anti-V1-V2 or V3 loop antibodies may correlate with the greater frequency of DM virus populations in HIV-1B as compared to HIV-1C infected individuals.      

While most HIV-1 infected individuals develop nAbs against their circulating variants and viruses escape this humoral pressure, there are potentially two important factors that explain why CXCR4-using variants do not evolve in everyone. First, CXCR4-utilizing viruses will not emerge if the circulating env variants cannot accommodate changes required for CXCR4 usage regardless of the selection pressure. Both laboratory and human studies have documented that some HIV-1 env variants do not evolve to use the CXCR4 receptor even in setting where they are exposed to CCR5 inhibitors ([[#bb0405|Westby et al., 2007]] , [[#bb0235|Marozsan et al., 2005]]  and [[#bb0385|Tsibris et al., 2008]] ). Second, similar to the observation that some but not all HIV-1 infected individuals develop cross-neutralizing antibodies ([[#bb0085|Doria-Rose et al., 2009]] ), co-receptor switching due to humoral selection may only occur in people that develop certain types of antibodies over the course of infection. Thus, in people that do not generate antibodies with the appropriate specificity neutralization escape may occur without changes in co-receptor usage. Unfortunately, we did not have access to live PBMCs to generate mAbs from these individuals, which would have allowed us to map the epitope for the antibodies that neutralize R5 but not co-existing X4 variants.      

Interestingly, treatment of HIV-1 infected humanized mice with anti-V3 loop bnAb, PGT128, but not anti-CD4bs bnAb, 45-46GW, has been predicted to select for a significantly higher proportion of CXCR4-using variants ([[#bb0300|Pfeifer et al., 2014]] ). Furthermore, CXCR4 as compared to CCR5-using variants are overrepresented among variants resistant to env V1–V2 directed bnAbs, PG9 and PG16. Thus similar to our suggestions, this investigation implies that antibodies directed against the env variable loops potentially select for neutralization resistant CXCR4-utilizing variants ([[#bb0300|Pfeifer et al., 2014]] ). One study has argued that this may also be the case during longitudinal HIV-2 infection ([[#bb0230|Marcelino et al., 2012]] ). Even though our conclusions are based on observations from a small number of individuals, and in some cases only a limited number of envs with different co-receptor phenotype, our studies strongly support the idea that X4 variants are neutralization escape variants.      

In contrast to the previously proposed mechanisms that lead to co-receptor switching, such as random mutations, waning immunity, and/or decreasing number of CCR5 + T cells, our observations suggest that antibody pressure selects for CXCR4-utilizing variants ([[#bb0315|Regoes and Bonhoeffer, 2005]] ). It is likely that in some but not necessarily all individuals that eventually develop CXCR4-using viruses antibodies force env structural changes that directly influence co-receptor switching. Continual presence of these humoral responses eventually leads to the emergence of X4 variants resistant to these and other host generated anti V3 antibodies. This new biological mechanism for co-receptor switching has a number of important implications. It is generally agreed that presence of CXCR4-using variants leads to faster CD4 decline and earlier onset of AIDS ([[#bb0190|Koot et al., 1993]] , [[#bb0325|Richman and Bozzette, 1994]]  and [[#bb0065|Connor et al., 1997]] ), although the causal link is unclear. Our studies would argue that an infected individuals humoral immune response leads to the emergence of CXCR4-utilizing variants, which results in an adverse outcome. Thus, our observations imply that CXCR4-using variants cause and are not a result of decreasing immunity. Importantly, this suggests that in some individuals, the humoral immune response against their autologous virus worsens infection outcome. Although antibodies have been suggested to adversely affect disease progression in other viral infections, notably dengue and influenza, no previous study has implied this for HIV-1 ([[#bb0145|Halstead and O'Rourke, 1977]]  and [[#bb0185|Khurana et al., 2013]] ). In addition, in contrast to dengue or influenza where antibodies enhance entry of the infecting virus into host cells, our studies suggest a mechanism where antibodies force phenotypic changes in replicating variants to expand the cellular tropism. Finally, our results have implications for proposed strategies of using bnAbs as prophylaxis, treatment, or cure ([[#bb0035|Caskey et al., 2015]]  and [[#bb0225|Lynch et al., 2015]] ). Our results and previous studies suggest that specific antibodies may be less potent against CXCR4-utilizing variants, which may select for their emergence in settings where antibodies are unable to fully inhibit virus replication ([[#bb0300|Pfeifer et al., 2014]] ). Thus, the bnAbs selected for immunotherapy will need to be evaluated carefully to ensure that they do not promote the emergence of CXCR4-using strains, which can accelerate disease progression among infected individuals.      

==Funding Sources==

This project was supported by NIH  grants [[#gts0005|AI081545]]  (NL) and [[#gts0005|AI122209]]  (MS) and BU Pilot Grant  (MS). No author received compensation from a pharmaceutical company or other agency for completion of this manuscript.      

==Conflict of Interest Statement==

There are no conflicts of interest associated with this study. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

==Author Contributions==

NL, XW, DRK, and MS designed the study. NL, OG, CB, LR, BE, and HL performed all the experiments. SL, ME, and SM oversaw the original trial and sample collection. MS wrote the paper with input from co-authors.

==Acknowledgements==

We would like to thank all Mashi trial participants, without whom this study would not have been possible.

==Appendix A. Supplementary Data==


[[Image:draft_García_292415893-mmc_pdf.gif|center|link=Lin_et_al_2016-mmc1.pdf]]
<span style="text-align: center; font-size: 75%;">

Model for antibody induced emergence of CXCR4-using virus

</span>

==References==

<ol style='list-style-type: none;margin-left: 0px;'><li><span id='bb0005'></span>
[[#bb0005|Ataher et al., 2012]] Q. Ataher, S. Portsmouth, L.A. Napolitano, S. ENG, A. Greenacre, A. Kambugu, R. Wood, S. Badal-Faesen, R. Tressler; The epidemiology and clinical correlates of HIV-1 co-receptor tropism in non-subtype B infections from India, Uganda and South Africa; J. Int. AIDS Soc., 15 (2012), p. 2</li>
<li><span id='bb0010'></span>
[[#bb0010|Blaak et al., 2000]] H. Blaak, A.B. Van't Wout, M. Brouwer, B. Hooibrink, E. Hovenkamp, H. Schuitemaker; In vivo HIV-1 infection of CD45RA(+)CD4(+) T cells is established primarily by syncytium-inducing variants and correlates with the rate of CD4(+) T cell decline; Proc. Natl. Acad. Sci. U. S. A., 97 (2000), pp. 1269–1274</li>
<li><span id='bb0015'></span>
[[#bb0015|Brumme et al., 2005]] Z.L. Brumme, J. Goodrich, H.B. Mayer, C.J. Brumme, B.M. Henrick, B. Wynhoven, J.J. Asselin, P.K. Cheung, R.S. Hogg, J.S. Montaner, P.R. Harrigan; Molecular and clinical epidemiology of CXCR4-using HIV-1 in a large population of antiretroviral-naive individuals; J. Infect. Dis., 192 (2005), pp. 466–474</li>
<li><span id='bb0020'></span>
[[#bb0020|Bunnik et al., 2007]] E.M. Bunnik, E.D. Quakkelaar, A.C. van Nuenen, B. Boeser-Nunnink, H. Schuitemaker; Increased neutralization sensitivity of recently emerged CXCR4-using human immunodeficiency virus type 1 strains compared to coexisting CCR5-using variants from the same patient; J. Virol., 81 (2007), pp. 525–531</li>
<li><span id='bb0025'></span>
[[#bb0025|Cann et al., 1992]] A.J. Cann, M.J. Churcher, M. Boyd, W. O'Brien, J.Q. Zhao, J. Zack, I.S. Chen; The region of the envelope gene of human immunodeficiency virus type 1 responsible for determination of cell tropism; J. Virol., 66 (1992), pp. 305–309</li>
<li><span id='bb0030'></span>
[[#bb0030|Cashin et al., 2014]] K. Cashin, G. Paukovics, M.R. Jakobsen, L. Ostergaard, M.J. Churchill, P.R. Gorry, J.K. Flynn; Differences in coreceptor specificity contribute to alternative tropism of HIV-1 subtype C for CD4(+) T-cell subsets, including stem cell memory T-cells; Retrovirology, 11 (2014), p. 97</li>
<li><span id='bb0035'></span>
[[#bb0035|Caskey et al., 2015]] M. Caskey, F. Klein, J.C. Lorenzi, M.S. Seaman, A.P. West JR., N. Buckley, G. Kremer, L. Nogueira, M. Braunschweig, J.F. Scheid, J.A. Horwitz, I. Shimeliovich, S. Ben-Avraham, M. Witmer-Pack, M. Platten, C. Lehmann, L.A. Burke, T. Hawthorne, R.J. Gorelick, B.D. Walker, T. Keler, R.M. Gulick, G. Fatkenheuer, S.J. Schlesinger, M.C. Nussenzweig; Viraemia suppressed in HIV-1-infected humans by broadly neutralizing antibody 3BNC117; Nature, 522 (2015), pp. 487–491</li>
<li><span id='bb0040'></span>
[[#bb0040|Cecilia et al., 1998]] D. Cecilia, V.N. Kewalramani, J. O'Leary, B. Volsky, P. Nyambi, S. Burda, S. Xu, D.R. Littman, S. Zolla-Pazner; Neutralization profiles of primary human immunodeficiency virus type 1 isolates in the context of coreceptor usage; J. Virol., 72 (1998), pp. 6988–6996</li>
<li><span id='bb0045'></span>
[[#bb0045|Chatziandreou et al., 2012]] N. Chatziandreou, A.B. Arauz, I. Freitas, P.H. Nyein, G. Fenton, S.H. Mehta, G.D. Kirk, M. Sagar; Sensitivity changes over the course of infection increases the likelihood of resistance against fusion but not CCR5 receptor blockers; AIDS Res. Hum. Retrovir., 28 (2012), pp. 1584–1593</li>
<li><span id='bb0050'></span>
[[#bb0050|Chesebro et al., 1991]] B. Chesebro, J. Nishio, S. Perryman, A. Cann, W. O'Brien, I.S. Chen, K. Wehrly; Identification of human immunodeficiency virus envelope gene sequences influencing viral entry into CD4-positive HeLa cells, T-leukemia cells, and macrophages; J. Virol., 65 (1991), pp. 5782–5789</li>
<li><span id='bb0055'></span>
[[#bb0055|Chesebro et al., 1996]] B. Chesebro, K. Wehrly, J. Nishio, S. Perryman; Mapping of independent V3 envelope determinants of human immunodeficiency virus type 1 macrophage tropism and syncytium formation in lymphocytes; J. Virol., 70 (1996), pp. 9055–9059</li>
<li><span id='bb0060'></span>
[[#bb0060|Cocchi et al., 1996]] F. Cocchi, A.L. Devico, A. Garzino-Demo, A. Cara, R.C. Gallo, P. Lusso; The V3 domain of the HIV-1 gp120 envelope glycoprotein is critical for chemokine-mediated blockade of infection; Nat. Med., 2 (1996), pp. 1244–1247</li>
<li><span id='bb0065'></span>
[[#bb0065|Connor et al., 1997]] R.I. Connor, K.E. Sheridan, D. Ceradini, S. Choe, N.R. Landau; Change in coreceptor use correlates with disease progression in HIV-1-infected individuals; J. Exp. Med., 185 (1997), pp. 621–628</li>
<li><span id='bb0070'></span>
[[#bb0070|Davenport et al., 2002]] M.P. Davenport, J.J. Zaunders, M.D. Hazenberg, H. Schuitemaker, R.P. Van Rij; Cell turnover and cell tropism in HIV-1 infection; Trends Microbiol., 10 (2002), pp. 275–278</li>
<li><span id='bb0075'></span>
[[#bb0075|del Prete et al., 2009]] G.Q. del Prete, B. Haggarty, G.J. Leslie, A.P. Jordan, J. Romano, N. Wang, J. Wang, M.C. Holmes, D.C. Montefiori, J.A. Hoxie; Derivation and characterization of a simian immunodeficiency virus SIVmac239 variant with tropism for CXCR4; J. Virol., 83 (2009), pp. 9911–9922</li>
<li><span id='bb0080'></span>
[[#bb0080|Deng et al., 1996]] H. Deng, R. Liu, W. Ellmeier, S. Choe, D. Unutmaz, M. Burkhart, P. di Marzio, S. Marmon, R.E. Sutton, C.M. Hill, C.B. Davis, S.C. Peiper, T.J. Schall, D.R. Littman, N.R. Landau; Identification of a major co-receptor for primary isolates of HIV-1; Nature, 381 (1996), pp. 661–666</li>
<li><span id='bb0085'></span>
[[#bb0085|Doria-Rose et al., 2009]] N.A. Doria-Rose, R.M. Klein, M.M. Manion, S. O'Dell, A. Phogat, B. Chakrabarti, C.W. Hallahan, S.A. Migueles, J. Wrammert, R. Ahmed, M. Nason, R.T. Wyatt, J.R. Mascola, M. Connors; Frequency and phenotype of human immunodeficiency virus envelope-specific B cells from patients with broadly cross-neutralizing antibodies; J. Virol., 83 (2009), pp. 188–199</li>
<li><span id='bb0090'></span>
[[#bb0090|D'Souza et al., 1997]] M.P. D'Souza, D. Livnat, J.A. Bradac, S.H. Bridges; Evaluation of monoclonal antibodies to human immunodeficiency virus type 1 primary isolates by neutralization assays: performance criteria for selecting candidate antibodies for clinical trials. AIDS Clinical Trials Group Antibody Selection Working Group; J. Infect. Dis., 175 (1997), pp. 1056–1062</li>
<li><span id='bb0095'></span>
[[#bb0095|Etemad et al., 2009]] B. Etemad, A. Fellows, B. Kwambana, A. Kamat, Y. Feng, S. Lee, M. Sagar; Human immunodeficiency virus type 1 V1-to-V5 envelope variants from the chronic phase of infection use CCR5 and fuse more efficiently than those from early after infection; J. Virol., 83 (2009), pp. 9694–9708</li>
<li><span id='bb0100'></span>
[[#bb0100|Etemad et al., 2014]] B. Etemad, O.A. Gonzalez, L. White, O. Laeyendecker, G.D. Kirk, S. Mehta, M. Sagar; Characterization of HIV-1 envelopes in acutely and chronically infected injection drug users; Retrovirology, 11 (2014), p. 106</li>
<li><span id='bb0105'></span>
[[#bb0105|Etemad et al., 2015]] B. Etemad, M. Ghulam-Smith, O. Gonzalez, L.F. White, M. Sagar; Single genome amplification and standard bulk PCR yield HIV-1 envelope products with similar genotypic and phenotypic characteristics; J. Virol. Methods, 214 (2015), pp. 46–53</li>
<li><span id='bb0110'></span>
[[#bb0110|Fang et al., 1998]] G. Fang, G. Zhu, H. Burger, J.S. Keithly, B. Weiser; Minimizing DNA recombination during long RT-PCR; J. Virol. Methods, 76 (1998), pp. 139–148</li>
<li><span id='bb0115'></span>
[[#bb0115|Feng et al., 1996]] Y. Feng, C.C. Broder, P.E. Kennedy, E.A. Berger; HIV-1 entry cofactor: functional cDNA cloning of a seven-transmembrane, G protein-coupled receptor; Science, 272 (1996), pp. 872–877</li>
<li><span id='bb0120'></span>
[[#bb0120|Fouchier et al., 1995]] R.A. Fouchier, S.M. Broersen, M. Brouwer, M. Tersmette, A.B. van't Wout, M. Groenink, H. Schuitemaker; Temporal relationship between elongation of the HIV type 1 glycoprotein 120 V2 domain and the conversion toward a syncytium-inducing phenotype; AIDS Res. Hum. Retrovir., 11 (1995), pp. 1473–1478</li>
<li><span id='bb0125'></span>
[[#bb0125|Freed and Martin, 1995]] E.O. Freed, M.A. Martin; Virion incorporation of envelope glycoproteins with long but not short cytoplasmic tails is blocked by specific, single amino acid substitutions in the human immunodeficiency virus type 1 matrix; J. Virol., 69 (1995), pp. 1984–1989</li>
<li><span id='bb0130'></span>
[[#bb0130|Galvani and Novembre, 2005]] A.P. Galvani, J. Novembre; The evolutionary history of the CCR5-Delta32 HIV-resistance mutation; Microbes Infect., 7 (2005), pp. 302–309</li>
<li><span id='bb0135'></span>
[[#bb0135|Gorny et al., 1992]] M.K. Gorny, A.J. Conley, S. Karwowska, A. Buchbinder, J.Y. Xu, E.A. Emini, S. Koenig, S. Zolla-Pazner; Neutralization of diverse human immunodeficiency virus type 1 variants by an anti-V3 human monoclonal antibody; J. Virol., 66 (1992), pp. 7538–7542</li>
<li><span id='bb0140'></span>
[[#bb0140|Groenink et al., 1993]] M. Groenink, R.A. Fouchier, S. Broersen, C.H. Baker, M. Koot, A.B. van't Wout, H.G. Huisman, F. Miedema, M. Tersmette, H. Schuitemaker; Relation of phenotype evolution of HIV-1 to envelope V2 configuration; Science, 260 (1993), pp. 1513–1516</li>
<li><span id='bb0145'></span>
[[#bb0145|Halstead and O'Rourke, 1977]] S.B. Halstead, E.J. O'Rourke; Dengue viruses and mononuclear phagocytes. I. Infection enhancement by non-neutralizing antibody; J. Exp. Med., 146 (1977), pp. 201–217</li>
<li><span id='bb0150'></span>
[[#bb0150|Henrich et al., 2014]] T.J. Henrich, P.J. McLaren, S.S. Rao, N.H. Lin, E. Hanhauser, F. Giguel, R.M. Gulick, H. Ribaudo, P.I. de Bakker, D.R. Kuritzkes; Genome-wide association study of human immunodeficiency virus (HIV)-1 coreceptor usage in treatment-naive patients from an AIDS clinical trials group study; Open Forum Infect. Dis., 1 (2014), p. ofu018</li>
<li><span id='bb0155'></span>
[[#bb0155|Ho et al., 2007]] S.H. Ho, S. Tasca, L. Shek, A. Li, A. Gettie, J. Blanchard, D. Boden, C. Cheng-Mayer; Coreceptor switch in R5-tropic simian/human immunodeficiency virus-infected macaques; J. Virol., 81 (2007), pp. 8621–8633</li>
<li><span id='bb0160'></span>
[[#bb0160|Huang et al., 2007]] W. Huang, S.H. Eshleman, J. Toma, S. Fransen, E. Stawiski, E.E. Paxinos, J.M. Whitcomb, A.M. Young, D. Donnell, F. Mmiro, P. Musoke, L.A. Guay, J.B. Jackson, N.T. Parkin, C.J. Petropoulos; Coreceptor tropism in human immunodeficiency virus type 1 subtype D: high prevalence of CXCR4 tropism and heterogeneous composition of viral populations; J. Virol., 81 (2007), pp. 7885–7893</li>
<li><span id='bb0165'></span>
[[#bb0165|Huang et al., 2012]] J. Huang, G. Ofek, L. Laub, M.K. Louder, N.A. Doria-Rose, N.S. Longo, H. Imamichi, R.T. Bailer, B. Chakrabarti, S.K. Sharma, S.M. Alam, T. Wang, Y. Yang, B. Zhang, S.A. Migueles, R. Wyatt, B.F. Haynes, P.D. Kwong, J.R. Mascola, M. Connors; Broad and potent neutralization of HIV-1 by a gp41-specific human antibody; Nature, 491 (2012), pp. 406–412</li>
<li><span id='bb0170'></span>
[[#bb0170|Jia et al., 2016]] M. Jia, H. Lu, M. Markowitz, C. Cheng-Mayer, X. Wu; Development of broadly neutralizing antibodies and their mapping by monomeric gp120 in HIV-1 infected humans and SHIVSF162P3N infected macaques; J. Virol., 90 (2016), pp. 4017–4031</li>
<li><span id='bb0175'></span>
[[#bb0175|Judo et al., 1998]] M.S. Judo, A.B. Wedel, C. Wilson; Stimulation and suppression of PCR-mediated recombination; Nucleic Acids Res., 26 (1998), pp. 1819–1825</li>
<li><span id='bb0180'></span>
[[#bb0180|Keele et al., 2008]] B.F. Keele, E.E. Giorgi, J.F. Salazar-Gonzalez, J.M. Decker, K.T. Pham, M.G. Salazar, C. Sun, T. Grayson, S. Wang, H. Li, X. Wei, C. Jiang, J.L. Kirchherr, F. Gao, J.A. Anderson, L.H. Ping, R. Swanstrom, G.D. Tomaras, W.A. Blattner, P.A. Goepfert, J.M. Kilby, M.S. Saag, E.L. Delwart, M.P. Busch, M.S. Cohen, D.C. Montefiori, B.F. Haynes, B. Gaschen, G.S. Athreya, H.Y. Lee, N. Wood, C. Seoighe, A.S. Perelson, T. Bhattacharya, B.T. Korber, B.H. Hahn, G.M. Shaw; Identification and characterization of transmitted and early founder virus envelopes in primary HIV-1 infection; Proc. Natl. Acad. Sci. U. S. A., 105 (2008), pp. 7552–7557</li>
<li><span id='bb0185'></span>
[[#bb0185|Khurana et al., 2013]] S. Khurana, C.L. Loving, J. Manischewitz, L.R. King, P.C. Gauger, J. Henningson, A.L. Vincent, H. Golding; Vaccine-induced anti-HA2 antibodies promote virus fusion and enhance influenza virus respiratory disease; Sci. Transl. Med., 5 (2013), p. 200ra114</li>
<li><span id='bb0190'></span>
[[#bb0190|Koot et al., 1993]] M. Koot, I.P. Keet, A.H. Vos, R.E. de Goede, M.T. Roos, R.A. Coutinho, F. Miedema, P.T. Schellekens, M. Tersmette; Prognostic value of HIV-1 syncytium-inducing phenotype for rate of CD4 + cell depletion and progression to AIDS; Ann. Intern. Med., 118 (1993), pp. 681–688</li>
<li><span id='bb0195'></span>
[[#bb0195|Lacasse et al., 1998]] R.A. Lacasse, K.E. Follis, T. Moudgil, M. Trahey, J.M. Binley, V. Planelles, S. Zolla-Pazner, J.H. Nunberg; Coreceptor utilization by human immunodeficiency virus type 1 is not a primary determinant of neutralization sensitivity; J. Virol., 72 (1998), pp. 2491–2495</li>
<li><span id='bb0200'></span>
[[#bb0200|Lahr and Katz, 2009]] D.J. Lahr, L.A. Katz; Reducing the impact of PCR-mediated recombination in molecular evolution and environmental studies using a new-generation high-fidelity DNA polymerase; Biotechniques, 47 (2009), pp. 857–866</li>
<li><span id='bb0205'></span>
[[#bb0205|Lin et al., 2010]] N.H. Lin, D.M. Negusse, R. Beroukhim, F. Giguel, S. Lockman, M. Essex, D.R. Kuritzkes; The design and validation of a novel phenotypic assay to determine HIV-1 coreceptor usage of clinical isolates; J. Virol. Methods, 169 (2010), pp. 39–46</li>
<li><span id='bb0210'></span>
[[#bb0210|Lin et al., 2011]] N.H. Lin, L.M. Smeaton, F. Giguel, V. Novitsky, S. Moyo, R.M. Mitchell, J. Makhema, M. Essex, S. Lockman, D.R. Kuritzkes; Prevalence and clinical associations of CXCR4-using HIV-1 among treatment-naive subtype C-infected women in Botswana; J. Acquir. Immune Defic. Syndr., 57 (2011), pp. 46–50</li>
<li><span id='bb0215'></span>
[[#bb0215|Lin et al., 2012]] N.H. Lin, C. Becerril, F. Giguel, V. Novitsky, S. Moyo, J. Makhema, M. Essex, S. Lockman, D.R. Kuritzkes, M. Sagar; Env sequence determinants in CXCR4-using human immunodeficiency virus type-1 subtype C; Virology, 433 (2012), pp. 296–307</li>
<li><span id='bb0220'></span>
[[#bb0220|Liu et al., 1996]] S.L. Liu, A.G. Rodrigo, R. Shankarappa, G.H. Learn, L. Hsu, O. Davidov, L.P. Zhao, J.I. Mullins; HIV quasispecies and resampling; Science, 273 (1996), pp. 415–416</li>
<li><span id='bb0225'></span>
[[#bb0225|Lynch et al., 2015]] R.M. Lynch, E. Boritz, E.E. Coates, A. Dezure, P. Madden, P. Costner, M.E. Enama, S. Plummer, L. Holman, C.S. Hendel, I. Gordon, J. Casazza, M. Conan-Cibotti, S.A. Migueles, R. Tressler, R.T. Bailer, A. Mcdermott, S. Narpala, S. O'Dell, G. Wolf, J.D. Lifson, B.A. Freemire, R.J. Gorelick, J.P. Pandey, S. Mohan, N. Chomont, R. Fromentin, T.W. Chun, A.S. Fauci, R.M. Schwartz, R.A. Koup, D.C. Douek, Z. Hu, E. Capparelli, B.S. Graham, J.R. Mascola, J.E. Ledgerwood; Virologic effects of broadly neutralizing antibody VRC01 administration during chronic HIV-1 infection; Sci. Transl. Med., 7 (2015), p. 319ra206</li>
<li><span id='bb0230'></span>
[[#bb0230|Marcelino et al., 2012]] J.M. Marcelino, P. Borrego, C. Nilsson, C. Familia, H. Barroso, F. Maltez, M. Doroana, F. Antunes, A. Quintas, N. Taveira; Resistance to antibody neutralization in HIV-2 infection occurs in late stage disease and is associated with X4 tropism; AIDS, 26 (2012), pp. 2275–2284</li>
<li><span id='bb0235'></span>
[[#bb0235|Marozsan et al., 2005]] A.J. Marozsan, S.E. Kuhmann, T. Morgan, C. Herrera, E. Rivera-Troche, S. Xu, B.M. Baroudy, J. Strizki, J.P. Moore; Generation and properties of a human immunodeficiency virus type 1 isolate resistant to the small molecule CCR5 inhibitor, SCH-417690 (SCH-D); Virology, 338 (2005), pp. 182–199</li>
<li><span id='bb0240'></span>
[[#bb0240|Mascola et al., 1996]] J.R. Mascola, S. Snyder, O. Weislow, S. Belay, R. Belshe, D. Schwartz, M. Clements, R. Dolin, B. Graham, G. Gorse, M. Keefer, M. Mcelrath, M. Walker, K. Wagner, J. Mcneil, F. Mccutchan, D. Burke; Immunization with envelope subunit vaccine products elicits neutralizing antibodies against laboratory-adapted but not primary isolates of human immunodeficiency virus type 1. The National Institute of Allergy and Infectious Diseases AIDS Vaccine Evaluation Group; J. Infect. Dis., 173 (1996), pp. 340–348</li>
<li><span id='bb0245'></span>
[[#bb0245|Meyerhans et al., 1990]] A. Meyerhans, J.P. Vartanian, S. Wain-Hobson; DNA recombination during PCR; Nucleic Acids Res., 18 (1990), pp. 1687–1691</li>
<li><span id='bb0250'></span>
[[#bb0250|Montefiori et al., 1998]] D.C. Montefiori, R.G. Collman, T.R. Fouts, J.Y. Zhou, M. Bilska, J.A. Hoxie, J.P. Moore, D.P. Bolognesi; Evidence that antibody-mediated neutralization of human immunodeficiency virus type 1 by sera from infected individuals is independent of coreceptor usage; J. Virol., 72 (1998), pp. 1886–1893</li>
<li><span id='bb0255'></span>
[[#bb0255|Moody et al., 2015]] M.A. Moody, F. Gao, T.C. Gurley, J.D. Amos, A. Kumar, B. Hora, D.J. Marshall, J.F. Whitesides, S.M. Xia, R. Parks, K.E. Lloyd, K.K. Hwang, X. Lu, M. Bonsignori, A. Finzi, N.A. Vandergrift, S.M. Alam, G. Ferrari, X. Shen, G.D. Tomaras, G. Kamanga, M.S. Cohen, N.E. Sam, S. Kapiga, E.S. Gray, N.L. Tumba, L. Morris, S. Zolla-Pazner, M.K. Gorny, J.R. Mascola, B.H. Hahn, G.M. Shaw, J.G. Sodroski, H.X. Liao, D.C. Montefiori, P.T. Hraber, B.T. Korber, B.F. Haynes; Strain-specific V3 and CD4 binding site autologous HIV-1 neutralizing antibodies select neutralization-resistant viruses; Cell Host Microbe, 18 (2015), pp. 354–362</li>
<li><span id='bb0260'></span>
[[#bb0260|Moore et al., 2008]] P.L. Moore, E.S. Gray, I.A. Choge, N. Ranchobe, K. Mlisana, S.S. Abdool Karim, C. Williamson, L. Morris; The c3-v4 region is a major target of autologous neutralizing antibodies in human immunodeficiency virus type 1 subtype C infection; J. Virol., 82 (2008), pp. 1860–1869</li>
<li><span id='bb0265'></span>
[[#bb0265|Moore et al., 2009]] P.L. Moore, N. Ranchobe, B.E. Lambson, E.S. Gray, E. Cave, M.R. Abrahams, G. Bandawe, K. Mlisana, S.S. Abdool Karim, C. Williamson, L. Morris; Limited neutralizing antibody specificities drive neutralization escape in early HIV-1 subtype C infection; PLoS Pathog., 5 (2009), p. e1000598</li>
<li><span id='bb0270'></span>
[[#bb0270|O'Brien et al., 1990]] W.A. O'Brien, Y. Koyanagi, A. Namazie, J.Q. Zhao, A. Diagne, K. Idler, J.A. Zack, I.S. Chen; HIV-1 tropism for mononuclear phagocytes can be determined by regions of gp120 outside the CD4-binding domain; Nature, 348 (1990), pp. 69–73</li>
<li><span id='bb0275'></span>
[[#bb0275|Olshevsky et al., 1990]] U. Olshevsky, E. Helseth, C. Furman, J. Li, W. Haseltine, J. Sodroski; Identification of individual human immunodeficiency virus type 1 gp120 amino acids important for CD4 receptor binding; J. Virol., 64 (1990), pp. 5701–5707</li>
<li><span id='bb0280'></span>
[[#bb0280|Pastore et al., 2004]] C. Pastore, A. Ramos, D.E. Mosier; Intrinsic obstacles to human immunodeficiency virus type 1 coreceptor switching; J. Virol., 78 (2004), pp. 7565–7574</li>
<li><span id='bb0285'></span>
[[#bb0285|Patel et al., 2008]] M.B. Patel, N.G. Hoffman, R. Swanstrom; Subtype-specific conformational differences within the V3 region of subtype B and subtype C human immunodeficiency virus type 1 Env proteins; J. Virol., 82 (2008), pp. 903–916</li>
<li><span id='bb0290'></span>
[[#bb0290|Pejchal et al., 2011]] R. Pejchal, K.J. Doores, L.M. Walker, R. Khayat, P.S. Huang, S.K. Wang, R.L. Stanfield, J.P. Julien, A. Ramos, M. Crispin, R. Depetris, U. Katpally, A. Marozsan, A. Cupo, S. Maloveste, Y. Liu, R. Mcbride, Y. Ito, R.W. Sanders, C. Ogohara, J.C. Paulson, T. Feizi, C.N. Scanlan, C.H. Wong, J.P. Moore, W.C. Olson, A.B. Ward, P. Poignard, W.R. Schief, D.R. Burton, I.A. Wilson; A potent and broad neutralizing antibody recognizes and penetrates the HIV glycan shield; Science, 334 (2011), pp. 1097–1103</li>
<li><span id='bb0295'></span>
[[#bb0295|Pena-Cruz et al., 2013]] V. Pena-Cruz, B. Etemad, N. Chatziandreou, P.H. Nyein, S. Stock, S.J. reynolds, O. Laeyendecker, D. Serwadda, R.H. Gray, S. Lee, T.C. Quinn, M. Sagar; HIV-1 envelope replication and α4β7 utilization among newly infected subjects and their corresponding heterosexual partners; Retrovirology, 10 (2013), p. 162</li>
<li><span id='bb0300'></span>
[[#bb0300|Pfeifer et al., 2014]] N. Pfeifer, H. Walter, T. Lengauer; Association between HIV-1 coreceptor usage and resistance to broadly neutralizing antibodies; J. Acquir. Immune Defic. Syndr., 67 (2014), pp. 107–112</li>
<li><span id='bb0305'></span>
[[#bb0305|Picton et al., 2012]] A.C. Picton, S. Shalekoff, M. Paximadis, C.T. Tiemessen; Marked differences in CCR5 expression and activation levels in two South African populations; Immunology, 136 (2012), pp. 397–407</li>
<li><span id='bb0310'></span>
[[#bb0310|Provine et al., 2012]] N.M. Provine, V. Cortez, V. Chohan, J. Overbaugh; The neutralization sensitivity of viruses representing human immunodeficiency virus type 1 variants of diverse subtypes from early in infection is dependent on producer cell, as well as characteristics of the specific antibody and envelope variant; Virology, 427 (2012), pp. 25–33</li>
<li><span id='bb0315'></span>
[[#bb0315|Regoes and Bonhoeffer, 2005]] R.R. Regoes, S. Bonhoeffer; The HIV coreceptor switch: a population dynamical perspective; Trends Microbiol., 13 (2005), pp. 269–277</li>
<li><span id='bb0320'></span>
[[#bb0320|Ribeiro et al., 2006]] R.M. Ribeiro, M.D. Hazenberg, A.S. Perelson, M.P. Davenport; Naive and memory cell turnover as drivers of CCR5-to-CXCR4 tropism switch in human immunodeficiency virus type 1: implications for therapy; J. Virol., 80 (2006), pp. 802–809</li>
<li><span id='bb0325'></span>
[[#bb0325|Richman and Bozzette, 1994]] D.D. Richman, S.A. Bozzette; The impact of the syncytium-inducing phenotype of human immunodeficiency virus on disease progression; J. Infect. Dis., 169 (1994), pp. 968–974</li>
<li><span id='bb0330'></span>
[[#bb0330|Richman et al., 2003]] D.D. Richman, T. Wrin, S.J. Little, C.J. Petropoulos; Rapid evolution of the neutralizing antibody response to HIV type 1 infection; Proc. Natl. Acad. Sci. U. S. A., 100 (2003), pp. 4144–4149</li>
<li><span id='bb0335'></span>
[[#bb0335|Roederer et al., 2013]] M. Roederer, B.F. Keele, S.D. Schmidt, R.D. Mason, H.C. Welles, W. Fischer, C. LABRANCHE, K.E. Foulds, M.K. Louder, Z.Y. Yang, J.P. Todd, A.P. Buzby, L.V. Mach, L. Shen, K.E. Seaton, B.M. Ward, R.T. Bailer, R. Gottardo, W. Gu, G. Ferrari, S.M. Alam, T.N. Denny, D.C. Montefiori, G.D. Tomaras, B.T. Korber, M.C. Nason, R.A. Seder, R.A. Koup, N.L. Letvin, S.S. Rao, G.J. Nabel, J.R. Mascola; Immunological and virological mechanisms of vaccine-mediated protection against SIV and HIV; Nature, 505 (2013), pp. 502–508</li>
<li><span id='bb0340'></span>
[[#bb0340|Rong et al., 2009]] R. Rong, B. Li, R.M. Lynch, R.E. Haaland, M.K. Murphy, J. Mulenga, S.A. Allen, A. Pinter, G.M. Shaw, E. Hunter, J.E. Robinson, S. Gnanakaran, C.A. Derdeyn; Escape from autologous neutralizing antibodies in acute/early subtype C HIV-1 infection requires multiple pathways; PLoS Pathog., 5 (2009), p. e1000594</li>
<li><span id='bb0345'></span>
[[#bb0345|Sagar et al., 2006]] M. Sagar, X. Wu, S. Lee, J. Overbaugh; HIV-1 V1-V2 envelope loop sequences expand and add glycosylation sites over the course of infection and these modifications affect antibody neutralization sensitivity; J. Virol., 80 (2006), pp. 9586–9598</li>
<li><span id='bb0350'></span>
[[#bb0350|Salazar-Gonzalez et al., 2008]] J.F. Salazar-Gonzalez, E. Bailes, K.T. Pham, M.G. Salazar, M.B. Guffey, B.F. Keele, C.A. Derdeyn, P. Farmer, E. Hunter, S. Allen, O. Manigart, J. Mulenga, J.A. Anderson, R. Swanstrom, B.F. Haynes, G.S. Athreya, B.T. Korber, P.M. Sharp, G.M. Shaw, B.H. Hahn; Deciphering human immunodeficiency virus type 1 transmission and early envelope diversification by single-genome amplification and sequencing; J. Virol., 82 (2008), pp. 3952–3970</li>
<li><span id='bb0355'></span>
[[#bb0355|Sawyer et al., 1994]] L.S. Sawyer, M.T. Wrin, L. Crawford-Miksza, B. Potts, Y. Wu, P.A. Weber, R.D. Alfonso, C.V. Hanson; Neutralization sensitivity of human immunodeficiency virus type 1 is determined in part by the cell in which the virus is propagated; J. Virol., 68 (1994), pp. 1342–1349</li>
<li><span id='bb0360'></span>
[[#bb0360|Schuitemaker et al., 1992]] H. Schuitemaker, M. Koot, N.A. Kootstra, M.W. Dercksen, R.E. de Goede, R.P. van Steenwijk, J.M. Lange, J.K. Schattenkerk, F. Miedema, M. Tersmette; Biological phenotype of human immunodeficiency virus type 1 clones at different stages of infection: progression of disease is associated with a shift from monocytotropic to T-cell-tropic virus population; J. Virol., 66 (1992), pp. 1354–1360</li>
<li><span id='bb0365'></span>
[[#bb0365|Shioda et al., 1992]] T. Shioda, J.A. Levy, C. Cheng-Mayer; Small amino acid changes in the V3 hypervariable region of gp120 can affect the T-cell-line and macrophage tropism of human immunodeficiency virus type 1; Proc. Natl. Acad. Sci. U. S. A., 89 (1992), pp. 9434–9438</li>
<li><span id='bb0370'></span>
[[#bb0370|Tang et al., 2011]] H. Tang, J.E. Robinson, S. Gnanakaran, M. Li, E.S. Rosenberg, L.G. Perez, B.F. Haynes, H.X. Liao, C.C. Labranche, B.T. Korber, D.C. Montefiori; Epitopes immediately below the base of the V3 loop of gp120 as targets for the initial autologous neutralizing antibody response in two HIV-1 subtype B-infected individuals; J. Virol., 85 (2011), pp. 9286–9299</li>
<li><span id='bb0375'></span>
[[#bb0375|Thior et al., 2006]] I. Thior, S. Lockman, L.M. Smeaton, R.L. Shapiro, C. Wester, S.J. Heymann, P.B. Gilbert, L. Stevens, T. Peter, S. Kim, E. van Widenfelt, C. Moffat, P. Ndase, P. Arimi, P. Kebaabetswe, P. Mazonde, J. Makhema, K. Mcintosh, V. Novitsky, T.H. Lee, R. Marlink, S. Lagakos, M. Essex; Breastfeeding plus infant zidovudine prophylaxis for 6 months vs formula feeding plus infant zidovudine for 1 month to reduce mother-to-child HIV transmission in Botswana: a randomized trial: the Mashi Study; JAMA, 296 (2006), pp. 794–805</li>
<li><span id='bb0380'></span>
[[#bb0380|Trkola et al., 1998]] A. Trkola, T. Ketas, V.N. Kewalramani, F. Endorf, J.M. Binley, H. Katinger, J. Robinson, D.R. Littman, J.P. Moore; Neutralization sensitivity of human immunodeficiency virus type 1 primary isolates to antibodies and CD4-based reagents is independent of coreceptor usage; J. Virol., 72 (1998), pp. 1876–1885</li>
<li><span id='bb0385'></span>
[[#bb0385|Tsibris et al., 2008]] A.M. Tsibris, M. Sagar, R.M. Gulick, Z. Su, M. Hughes, W. Greaves, M. Subramanian, C. Flexner, F. Giguel, K.E. Leopold, E. Coakley, D.R. Kuritzkes; In vivo emergence of vicriviroc resistance in a human immunodeficiency virus type 1 subtype C-infected subject; J. Virol., 82 (2008), pp. 8210–8214</li>
<li><span id='bb0390'></span>
[[#bb0390|van't Wout et al., 1994]] A.B. van't Wout, N.A. Kootstra, G.A. Mulder-Kampinga, N. Albrecht-van Lent, H.J. Scherpbier, J. Veenstra, K. Boer, R.A. Coutinho, F. Miedema, H. Schuitemaker; Macrophage-tropic variants initiate human immunodeficiency virus type 1 infection after sexual, parenteral, and vertical transmission; J. Clin. Invest., 94 (1994), pp. 2060–2067</li>
<li><span id='bb0395'></span>
[[#bb0395|Walker et al., 2011]] L.M. Walker, M. Huber, K.J. Doores, E. Falkowska, R. Pejchal, J.P. Julien, S.K. Wang, A. Ramos, P.Y. Chan-Hui, M. Moyle, J.L. Mitcham, P.W. Hammond, O.A. Olsen, P. Phung, S. Fling, C.H. Wong, S. Phogat, T. Wrin, M.D. Simek, W.C. Koff, I.A. Wilson, D.R. Burton, P. Poignard; Broad neutralization coverage of HIV by multiple highly potent antibodies; Nature, 477 (2011), pp. 466–470</li>
<li><span id='bb0400'></span>
[[#bb0400|Wei et al., 2003]] X. Wei, J.M. Decker, S. Wang, H. Hui, J.C. Kappes, X. Wu, J.F. Salazar-Gonzalez, M.G. Salazar, J.M. Kilby, M.S. Saag, N.L. Komarova, M.A. Nowak, B.H. Hahn, P.D. Kwong, G.M. Shaw; Antibody neutralization and escape by HIV-1; Nature, 422 (2003), pp. 307–312</li>
<li><span id='bb0405'></span>
[[#bb0405|Westby et al., 2007]] M. Westby, C. Smith-Burchnell, J. Mori, M. Lewis, M. Mosley, M. Stockdale, P. Dorr, G. Ciaramella, M. Perros; Reduced maximal inhibition in phenotypic susceptibility assays indicates that viral strains resistant to the CCR5 antagonist maraviroc utilize inhibitor-bound receptor for entry; J. Virol., 81 (2007), pp. 2359–2371</li>
<li><span id='bb0410'></span>
[[#bb0410|Wu et al., 2010]] X. Wu, Z.Y. Yang, Y. Li, C.M. Hogerkorp, W.R. Schief, M.S. Seaman, T. Zhou, S.D. Schmidt, L. Wu, L. Xu, N.S. Longo, K. McKee, S. O'Dell, M.K. Louder, D.L. Wycuff, Y. Feng, M. Nason, N. Doria-Rose, M. Connors, P.D. Kwong, M. Roederer, R.T. Wyatt, G.J. Nabel, J.R. Mascola; Rational design of envelope identifies broadly neutralizing human monoclonal antibodies to HIV-1; Science, 329 (2010), pp. 856–861</li>
<li><span id='bb0415'></span>
[[#bb0415|Zhou et al., 2010]] T. Zhou, I. Georgiev, X. Wu, Z.Y. Yang, K. Dai, A. Finzi, Y.D. Kwon, J.F. Scheid, W. Shi, L. Xu, Y. Yang, J. Zhu, M.C. Nussenzweig, J. Sodroski, L. Shapiro, G.J. Nabel, J.R. Mascola, P.D. Kwong; Structural basis for broad and potent neutralization of HIV-1 by antibody VRC01; Science, 329 (2010), pp. 811–817</li>
</ol>