- Short report
- Open Access
A strongly selected mutation in the HIV-1 genome is independent of T cell responses and neutralizing antibodies
Retrovirologyvolume 14, Article number: 46 (2017)
Mutations rapidly accumulate in the HIV-1 genome after infection. Some of those mutations are selected by host immune responses and often cause viral fitness losses. This study is to investigate whether strongly selected mutations that are not associated with immune responses result in fitness losses.
Strongly selected mutations were identified by analyzing 5′-half HIV-1 genome (gag/pol) sequences from longitudinal samples of subject CH0131. The K43R mutation in the gag gene was first detected at day 91 post screening and was fixed in the viral population at day 273 while the synonymous N323tc mutation was first detected at day 177 and fixed at day 670. No conventional or cryptic T cell responses were detected against either mutation sites by ELISpot analysis. However, when fitness costs of both mutations were measured by introducing each mutation into their cognate transmitted/founder (T/F) viral genome, the K43R mutation caused a significant fitness loss while the N323tc mutation had little impact on viral fitness.
The rapid fixation, the lack of detectable immune responses and the significant fitness cost of the K43R mutation suggests that it was strongly selected by host factors other than T cell responses and neutralizing antibodies.
After infection, human immunodeficiency virus type 1 (HIV-1) rapidly evolves into a quasispecies population due to the high level of mutability of the virus and host selection pressure. Some of the most strongly selected mutations are associated with CD8+ T cell responses [1,2,3,4,5] or neutralizing antibodies (nAbs) [6, 7]. By precisely determining T cell responses using peptides based on the autologous transmitted/founder (T/F) virus proteome sequences in each infected individual and analyzing longitudinal near full-length HIV-1 genome sequences obtained by single genome amplification (SGA), we previously found that some predominant or fixed mutations were associated with T cell responses [5, 8]. However, we also observed that about one-third of rapidly predominating mutations were not associated with cytotoxic T lymphocyte (CTL), especially for those mutations not in the env gene . This indicates that host factors other than immune responses contribute to the selection of those mutations. Investigating the landscape of host selection pressure exerting on the HIV-1 genome will provide a better understanding of the complex virus-host interactions.
Mutations that have a significant fitness cost predominate because they have a replicative advantage over the wild type virus under strong host selection pressure. Mutations selected by CTL responses generally impair viral fitness [9,10,11,12,13]. Fitness costs of immune escape mutations can have important implications for viral pathogenesis, transmission and vaccine development [14,15,16,17,18,19,20,21]. By introducing select mutations into their cognate T/F viral genomes, fitness change can be precisely determined by a competitive fitness assay. Our previous studies showed that both mutations selected by T cell and nAb responses significantly affected the fitness of the T/F viruses [22,23,24,25]. However, whether mutations that are not associated with immune selection can affect viral fitness has not been studied.
Our previous analysis of longitudinal 5′-half genome sequences obtained by SGA from subject CH0131 showed that strong selection pressure from the host drove some mutations to be rapidly fixed in the viral population within 2 years of infection . We showed that mutations at two sites were selected by T cell immune responses  while one was a reversion mutation in the TW10 epitope . In this study, we found that the earliest mutation K43R in the gag gene in CH0131 was strongly selected and rapidly predominated in the viral population within 6 months of infection. No conventional or cryptic T cell responses were detected targeting this mutation site. However, the mutation caused a significant fitness loss. Thus, our results demonstrated that host factors other than adaptive immune responses can also have strong selection pressure on viral evolution and drive genetic diversification of HIV-1 in vivo.
Strongly selected mutations are identified in the 5′-half genome
Analysis of longitudinal 5′-half genome (gag/pol) sequences from subject CH0131 (HLA A*29:01, A*23:01, B*45:01, B*15:03, Cw*02:02, Cw*06:02) identified mutations that were fixed at five sites in the viral population by day 670 post screening (Fig. 1). Mutations (K43R, E161Q/S165N, N242T and N323tc) were at four sites in the gag gene and one mutation (A980V) was at the end of the pol gene. Using autologous overlapping peptides based on the T/F whole proteome sequence, we previously detected CD8+ T cell responses targeting two T cell epitopes VF9156–164 and RY10978–986 . By day 670, mutations E161Q/S165N and A980V in VF9156–164 and RY10978–986 epitopes were fixed in the viral population, suggesting that these mutations were selected by the T cell responses (Fig. 1). However, no T cell responses were detected for the regions containing the other three fixed mutations. The N242T mutation in the TW10240–249 epitope is a reversion mutation [26, 27] since the subject CH0131 did not have the HLA B*57 allele and was infected with the escape mutant virus (Fig. 1) . The N242T mutation alone had no detectable impact on fitness of the cognate T/F virus .
The K43R mutation, located in an HLA-B*35-restricted WF936–44 epitope (WASRELERF) in Gag , was the first detected mutation (day 91). It became predominant (83%) at day 175 and was fixed in the viral population at days 273–670. This indicated that this mutation was strongly selected in the host (Figs. 1 and 2a). The synonymous N323tc mutation in Gag was first detected at day 175 and fixed at day 670 (Figs. 1 and 2b).
The K43R and N323tc mutations are not selected by T cell responses
Since no T cell responses were detected against the WF9 epitope using 18-mer peptides with 10 amino acid overlap in our previous study , we further determined T cell responses using various shorter 9-mer peptides containing the K43 site by ELISpot. None of the peptides were recognized by the T cells (Fig. 3a). To determine if it was selected by cryptic T cell epitopes [29, 30] from alternative expression products in the HIV-1 genome, we examined all potential open reading frames (ORFs) and found that the A → G mutation in the K43 codon (AAA) resulted in amino acid substitutions in three alternative ORFs (Fig. 3b). To determine if these potentially coded proteins were targeted by T cell responses, we tested all possible peptides containing the K43R mutation for their ability to stimulate T cell responses. None of these peptides could be recognized by T cells from days 23, 93 and 121 (Fig. 3c). These results demonstrated that the K43R mutation was not selected by classic or cryptic T cell responses.
T cell responses targeting the site containing the synonymous N323tc mutation were also not detected with autologous Gag peptides . To further investigate if this mutation was selected by cryptic T cell responses, we tested nine peptides from four alternative ORFs (Fig. 4a) to determine if these potentially coded proteins were targeted by T cell responses. ELISpot analysis showed that none of those peptides were recognized by the T cells from days 65 and 602 (Fig. 4b). This suggests that the N323tc mutation was also not selected by classic or cryptic T cell responses.
The K43R mutation causes significant fitness loss
One important parameter to evaluate whether a mutation is selected by a host is to determine if it causes fitness loss. A fixed mutation with significant fitness loss will indicate a strong host selective pressure. To determine if the K43R and N323tc mutations affected viral fitness, we introduced each mutation into their cognate T/F infectious molecular clone (IMC)  and determined their fitness costs. The T/F and the two mutant viruses replicated equally well in primary CD4+ T cells when cultured individually (Fig. 5a). When the K43R mutant was cultured together with the T/F virus, the proportion of the K43R mutant gradually decreased while the T/F virus gradually increased in the culture. Our model analysis showed that the K43R mutant was 12% less fit than the cognate T/F virus (− 12% ± 1%; p = 0.003 by t test) (Fig. 5b). Additional passages showed that the proportion of the K43R mutant continuously decreased after each passage and was almost completely replaced by the T/F virus by passage 3 (Fig. 5c). However, no fitness differences between the N323tc mutant and the T/F virus were observed in the single passage assay (− 0.00% ± 0.01%) (Fig. 5d). After three passages, the proportion of the N323tc mutant only slightly decreased (Fig. 5e). These results demonstrated that the K43R mutation caused a significantly fitness loss while the N323tc mutation had little impact on viral fitness.
The K43R mutation quickly emerged and was rapidly fixed in the viral population after infection in CH0131, suggesting that it was strongly selected in the host. Because no conventional or cryptic T cell responses were detected for the regions containing the K43R mutation and that the Gag protein does not elicit neutralizing antibodies, the strong selection pressure on the K43R mutation could not be from adaptive T and B cell responses. We have shown clearly that strong immunodominant T cell responses exerted the strongest and most rapid selection pressure in our previous study . Using the same approach, we showed that T cell responses to the WF9 epitope were absent or at least below the level of detection in this highly sensitive assay. Therefore, if there were T cell responses targeting the K43R site that were undetectable (and far from immunodominant), it would be difficult to envision how it could exert any significant pressure to select such a rapid escape, with a fitness loss.
T cell responses were also not found to target the N323tc mutation site, but the mutation had little impact on viral fitness. The majority of the N323tc mutation was found in the sequences that had the K43R and N242T mutations (Fig. 1). Thus, the N323tc mutation was not likely selected by the host, but was instead a concurrent mutation together with the N242T or K43R mutation.
Analysis of the gag gene sequences in the Los Alamos HIV sequence database (https://www.hiv.lanl.gov) showed that R43 was present in 97.3% of the distinct viruses while K43 only accounted for 2.2% of the population. Since the K43R mutation was found in the WF9 epitope restricted by the B*35 allele that CH0131 did not carry, it was mostly likely a reversion mutation (to the consensus of subtype C viruses). However, the significant fitness loss of the K43R reversion mutation indicates that it was not simply reverted back to the consensus state due to the lack of allele-specific selection pressure, but instead it was strongly selected in vivo. We have previously found that three other reversion mutations (V247I and T242N in Gag as well as I64T in Tat) also did not cause detectable fitness differences [24, 25]. These results demonstrate that reversion mutations alone do not necessarily lead to more fit viruses. The other possibility is that K43R is a compensatory mutation. It is known that the fitness loss caused by some mutations can be restored by compensatory mutations [10, 12, 16, 22, 31, 32]. This is also not likely since the K43R mutation was the very first mutation that selected and fixed in the viral population before other mutations were detected in the 5′-half genome. Mutations in vpu, env and nef were also detected in the 3′-half genome sequences when the K43R mutation was detected . However, it is unlikely that a mutation in the gag gene can compensate fitness losses due to mutations in those genes. K43R should not be a compensatory mutation for any mutations in LTR because such mutations were not detected in the viral genome. Most importantly, since the role of compensatory mutations is to restore fitness loss caused by other mutations, they alone should not result in fitness losses. Otherwise, they will cause additional fitness damage.
Genetic polymorphisms at many sites in the viral genome have been found to be associated with NK cell responses [33, 34] and anti-Gag antibody responses [35, 36]. However, specific mutations associated with selection pressure from NK cells or anti-Gag Ab responses have not been clearly identified as those mutations selected by T cell responses for in-frame or cryptic epitopes or Ab responses against epitopes in Env [1,2,3,4,5,6,7, 37]. Thus, it is unlikely that the K43R mutation was selected by NK cell responses and anti-Gag antibody responses. It will be important to determine if the K43R mutation or similar mutations are selected by NK cell responses and anti-Gag antibody responses or by T cell responses targeting to the spliced epitopes . While only the K43R mutation was fully documented here, we found other possible examples of selected escapes that were not within an epitope stimulating a measurable T cell response in our previous studies [5, 8]. Therefore, this could be a more widespread phenomenon and this study could raise this as an issue for further study. Similar mutations should be studied further to increase our understanding on how such mutations are strongly selected by the host and if the significant fitness loss caused by the mutation is associated with lower viral loads due to the significantly reduced viral replicative capacity.
The K43R and N323tc mutants were obtained by introducing each mutation into the cognate T/F IMC as previously described .
Ex vivo IFN-γ ELISpot assay
T cell responses to the targeting sites at various days post screening in CH0131 were determined using an in vitro IFN-γ ELISpot assay as we previously described [5, 8]. Peripheral blood mononuclear cells (PBMCs) were obtained through leukophereses from healthy donors under clinical protocols approved by the Duke University Institutional Review Board. ELISpot data are expressed as the mean spot forming units (SFU) per million PBMCs (SFU/106 PBMC) ± SEM. Positive T cell responses were defined as: ≥ 50 SFU/million PBMC and > 4 times above background. All assays were performed in triplicate with a peptide concentration of 2 μg/ml.
PASS fitness assay
PBMCs were obtained through leukophereses from healthy donors under clinical protocols approved by the Duke University Institutional Review Board. CD4+ T cells were purified with the QuadroMACS Separator (Miltenyi Biotec, Auburn, CA) and the CD4+ T cell Isolation Kit II (Miltenyi Biotec, Auburn, CA). CD4+ T cells (5 × 105 cells) were infected with the mixture of two compared viruses (2.5 ng of p24 for each virus) both in the single passage and through multiple passages as described before . The viral cDNA was synthesized with the primer low3 (5′-TTTTTCCTAGGGGCCCTGCAATTT-3′; nucleotide (nt) 1998–2021 in HXB2). The PASS amplifying primers were the forward primer M6F2 (5′-Acry-CTCGACGCAGGACTCGGCTTGCTG-3′; nt 685–708) and the reverse primer M6R2 (5′-TCCTCCCACTCCCTGACATGCTGTCATCATTTC-3′; nt 1822–1854). The bases were determined with sequencing primers juxtaposed to the mutation sites. All assays were performed in triplicate.
Modeling relative fitness and statistical analyses
The relative fitness is determined by measuring the replication slope of each virus in the culture over time using the mathematical model and the statistical analysis was performed as previously described [22, 24].
human immunodeficiency virus type 1
single genome amplification
cytotoxic T lymphocyte
human leukocyte antigen
infectious molecular clone
peripheral blood mononuclear cell
Phillips RE, Rowland-Jones S, Nixon DF, Gotch FM, Edwards JP, Ogunlesi AO, Elvin JG, Rothbard JA, Bangham CR, Rizza CR, et al. Human immunodeficiency virus genetic variation that can escape cytotoxic T cell recognition. Nature. 1991;354:453–9.
Borrow P, Lewicki H, Wei X, Horwitz MS, Peffer N, Meyers H, Nelson JA, Gairin JE, Hahn BH, Oldstone MB, Shaw GM. Antiviral pressure exerted by HIV-1-specific cytotoxic T lymphocytes (CTLs) during primary infection demonstrated by rapid selection of CTL escape virus. Nat Med. 1997;3:205–11.
Goulder PJ, Phillips RE, Colbert RA, McAdam S, Ogg G, Nowak MA, Giangrande P, Luzzi G, Morgan B, Edwards A, et al. Late escape from an immunodominant cytotoxic T-lymphocyte response associated with progression to AIDS. Nat Med. 1997;3:212–7.
Jones NA, Wei X, Flower DR, Wong M, Michor F, Saag MS, Hahn BH, Nowak MA, Shaw GM, Borrow P. Determinants of human immunodeficiency virus type 1 escape from the primary CD8 + cytotoxic T lymphocyte response. J Exp Med. 2004;200:1243–56.
Goonetilleke N, Liu MK, Salazar-Gonzalez JF, Ferrari G, Giorgi E, Ganusov VV, Keele BF, Learn GH, Turnbull EL, Salazar MG. The first T cell response to transmitted/founder virus contributes to the control of acute viremia in HIV-1 infection. J Exp Med. 2009;206:1253–72.
Richman DD, Wrin T, Little SJ, Petropoulos CJ. Rapid evolution of the neutralizing antibody response to HIV type 1 infection. Proc Natl Acad Sci. 2003;100:4144–9.
Wei X, Decker JM, Wang S, Hui H, Kappes JC, Wu X, Salazar-Gonzalez JF, Salazar MG, Kilby JM, Saag MS. Antibody neutralization and escape by HIV-1. Nature. 2003;422:307–12.
Liu MK, Hawkins N, Ritchie AJ, Ganusov VV, Whale V, Brackenridge S, Li H, Pavlicek JW, Cai F, Rose-Abrahams M. Vertical T cell immunodominance and epitope entropy determine HIV-1 escape. J Clin Investig. 2013;123:380.
Maree AF, Keulen W, Boucher CA, De Boer RJ. Estimating relative fitness in viral competition experiments. J Virol. 2000;74:11067–72.
Martinez-Picado J, Prado JG, Fry EE, Pfafferott K, Leslie A, Chetty S, Thobakgale C, Honeyborne I, Crawford H, Matthews P, et al. Fitness cost of escape mutations in p24 Gag in association with control of human immunodeficiency virus type 1. J Virol. 2006;80:3617–23.
Miura T, Brockman MA, Schneidewind A, Lobritz M, Pereyra F, Rathod A, Block BL, Brumme ZL, Brumme CJ, Baker B, et al. HLA-B57/B*5801 human immunodeficiency virus type 1 elite controllers select for rare gag variants associated with reduced viral replication capacity and strong cytotoxic T-lymphocyte [corrected] recognition. J Virol. 2009;83:2743–55.
Troyer RM, McNevin J, Liu Y, Zhang SC, Krizan RW, Abraha A, Tebit DM, Zhao H, Avila S, Lobritz MA, et al. Variable fitness impact of HIV-1 escape mutations to cytotoxic T lymphocyte (CTL) response. PLoS Pathog. 2009;5:e1000365.
Brockman MA, Schneidewind A, Lahaie M, Schmidt A, Miura T, Desouza I, Ryvkin F, Derdeyn CA, Allen S, Hunter E, et al. Escape and compensation from early HLA-B57-mediated cytotoxic T-lymphocyte pressure on human immunodeficiency virus type 1 Gag alter capsid interactions with cyclophilin A. J Virol. 2007;81:12608–18.
Arien KK, Troyer RM, Gali Y, Colebunders RL, Arts EJ, Vanham G. Replicative fitness of historical and recent HIV-1 isolates suggests HIV-1 attenuation over time. AIDS. 2005;19:1555–64.
Troyer RM, Collins KR, Abraha A, Fraundorf E, Moore DM, Krizan RW, Toossi Z, Colebunders RL, Jensen MA, Mullins JI, et al. Changes in human immunodeficiency virus type 1 fitness and genetic diversity during disease progression. J Virol. 2005;79:9006–18.
Chopera DR, Woodman Z, Mlisana K, Mlotshwa M, Martin DP, Seoighe C, Treurnicht F, de Rosa DA, Hide W, Karim SA. Transmission of HIV-1 CTL escape variants provides HLA-mismatched recipients with a survival advantage. PLoS Pathog. 2008;4:e1000033.
Goepfert PA, Lumm W, Farmer P, Matthews P, Prendergast A, Carlson JM, Derdeyn CA, Tang J, Kaslow RA, Bansal A, et al. Transmission of HIV-1 Gag immune escape mutations is associated with reduced viral load in linked recipients. J Exp Med. 2008;205:1009–17.
Abraha A, Nankya IL, Gibson R, Demers K, Tebit DM, Johnston E, Katzenstein D, Siddiqui A, Herrera C, Fischetti L, et al. CCR5- and CXCR4-tropic subtype C human immunodeficiency virus type 1 isolates have a lower level of pathogenic fitness than other dominant group M subtypes: implications for the epidemic. J Virol. 2009;83:5592–605.
Crawford H, Lumm W, Leslie A, Schaefer M, Boeras D, Prado JG, Tang J, Farmer P, Ndung’u T, Lakhi S, et al. Evolution of HLA-B*5703 HIV-1 escape mutations in HLA-B*5703-positive individuals and their transmission recipients. J Exp Med. 2009;206:909–21.
Rodriguez MA, Ding M, Ratner D, Chen Y, Tripathy SP, Kulkarni SS, Chatterjee R, Tarwater PM, Gupta P. High replication fitness and transmission efficiency of HIV-1 subtype C from India: implications for subtype C predominance. Virology. 2009;385:416–24.
Shakirzyanova M, Ren W, Zhuang K, Tasca S, Cheng-Mayer C. Fitness disadvantage of transitional intermediates contributes to dynamic change in the infecting-virus population during coreceptor switch in R5 simian/human immunodeficiency virus-infected macaques. J Virol. 2010;84:12862–71.
Song H, Pavlicek JW, Cai F, Bhattacharya T, Li H, Iyer SS, Bar KJ, Decker JM, Goonetilleke N, Liu MK. Impact of immune escape mutations on HIV-1 fitness in the context of the cognate transmitted/founder genome. Retrovirology. 2012;9:1–14.
Gao F, Bonsignori M, Liao HX, Kumar A, Xia SM, Lu X, Cai F, Hwang KK, Song H, Zhou T, et al. Cooperation of B cell lineages in induction of HIV-1-broadly neutralizing antibodies. Cell. 2014;158:481–91.
Liu D, Zuo T, Hora B, Song H, Kong W, Yu X, Goonetilleke N, Bhattacharya T, Perelson AS, Haynes BF, et al. Preexisting compensatory amino acids compromise fitness costs of a HIV-1 T cell escape mutation. Retrovirology. 2014;11:101.
Song H, Hora B, Bhattacharya T, Goonetilleke N, Liu MK, Wiehe K, Li H, Iyer SS, McMichael AJ, Perelson AS, Gao F. Reversion and T cell escape mutations compensate the fitness loss of a CD8 + T cell escape mutant in their cognate transmitted/founder virus. PLoS ONE. 2014;9:e102734.
Kawashima Y, Pfafferott K, Frater J, Matthews P, Payne R, Addo M, Gatanaga H, Fujiwara M, Hachiya A, Koizumi H, et al. Adaptation of HIV-1 to human leukocyte antigen class I. Nature. 2009;458:641–5.
Schneidewind A, Brumme ZL, Brumme CJ, Power KA, Reyor LL, O’Sullivan K, Gladden A, Hempel U, Kuntzen T, Wang YE, et al. Transmission and long-term stability of compensated CD8 escape mutations. J Virol. 2009;83:3993–7.
Goulder P, Edwards A, Phillips R, McMichael A. Identification of a novel HIA-B* 3501-restricted cytotoxic T lymphocyte epitope using overlapping peptides. Aids. 1997;11:930–2.
Bansal A, Carlson J, Yan J, Akinsiku OT, Schaefer M, Sabbaj S, Bet A, Levy DN, Heath S, Tang J. CD8 T cell response and evolutionary pressure to HIV-1 cryptic epitopes derived from antisense transcription. J Exp Med. 2010;207:51–9.
Bansal A, Mann T, Sterrett S, Peng BJ, Bet A, Carlson JM, Goepfert PA. Enhanced Recognition of HIV-1 Cryptic Epitopes Restricted by HLA Class I Alleles Associated With a Favorable Clinical Outcome. J Acquir Immune Defic Syndr. 2015;70:1–8.
Crawford H, Prado JG, Leslie A, Hue S, Honeyborne I, Reddy S, van der Stok M, Mncube Z, Brander C, Rousseau C, et al. Compensatory mutation partially restores fitness and delays reversion of escape mutation within the immunodominant HLA-B*5703-restricted Gag epitope in chronic human immunodeficiency virus type 1 infection. J Virol. 2007;81:8346–51.
Kelleher AD, Long C, Holmes EC, Allen RL, Wilson J, Conlon C, Workman C, Shaunak S, Olson K, Goulder P, et al. Clustered mutations in HIV-1 gag are consistently required for escape from HLA-B27-restricted cytotoxic T lymphocyte responses. J Exp Med. 2001;193:375–86.
Holzemer A, Thobakgale CF, Jimenez Cruz CA, Garcia-Beltran WF, Carlson JM, van Teijlingen NH, Mann JK, Jaggernath M, Kang SG, Korner C, et al. Selection of an HLA-C*03:04-restricted HIV-1 p24 gag sequence variant is associated with viral escape from KIR2DL3+ natural killer cells: data from an observational cohort in South Africa. PLoS Med. 2015;12:e1001900 (discussion e1001900).
Alter G, Heckerman D, Schneidewind A, Fadda L, Kadie CM, Carlson JM, Oniangue-Ndza C, Martin M, Li B, Khakoo SI, et al. HIV-1 adaptation to NK-cell-mediated immune pressure. Nature. 2011;476:96–100.
French MA, Abudulai LN, Fernandez S. Isotype diversification of IgG antibodies to HIV Gag proteins as a therapeutic vaccination strategy for HIV infection. Vaccines (Basel). 2013;1:328–42.
Tjiam MC, Taylor JP, Morshidi MA, Sariputra L, Burrows S, Martin JN, Deeks SG, Tan DB, Lee S, Fernandez S, French MA. Viremic HIV controllers exhibit high plasmacytoid dendritic cell-reactive opsonophagocytic IgG antibody responses against HIV-1 p24 associated with greater antibody isotype diversification. J Immunol. 2015;194:5320–8.
Chung AW, Isitman G, Navis M, Kramski M, Center RJ, Kent SJ, Stratov I. Immune escape from HIV-specific antibody-dependent cellular cytotoxicity (ADCC) pressure. Proc Natl Acad Sci USA. 2011;108:7505–10.
Liepe J, Marino F, Sidney J, Jeko A, Bunting DE, Sette A, Kloetzel PM, Stumpf MP, Heck AJ, Mishto M. A large fraction of HLA class I ligands are proteasome-generated spliced peptides. Science. 2016;354:354–8.
FG conceived the study and designed experiments. DL, CW, BH, TZ, NG, MKPL, MB, GF performed experiments. DL, CW, BH, TZ, NG, MKPL, MB, GF, AJM analyzed data. TB, ASP developed the algorithm for determination of viral relative fitness and statistical analysis. DL, CW, FG wrote the manuscript. NG, MKPL, GF, AJM, TB, ASP edited the manuscript. All authors read and approved the final manuscript.
We thank Barton Haynes for helpful discussion and support.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Consent for publication
This manuscript does not contain personal data from any individuals.
Ethics approval and consent to participate
Peripheral blood mononuclear cells (PBMC) were obtained through leukophereses from healthy donors under clinical protocols approved by the Duke University Institutional Review Board.
This work was supported by a grant from National Natural Science Foundation of China (31670162), NIH grants to the Center for HIV/AIDS Vaccine Immunology (AI067854), the Center for HIV/AIDS Vaccine Immunology and Immunogen Discovery (AI100645) and NIH Grants AI028433 and OD011095 (A.S.P.).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.