Antigen-responsive CD4+ T cell clones contribute to the HIV-1 latent reservoir

1After integration into the host genome, HIV-1 transcription usually leads to new virus production and cell death. However, HIV-1 can also become latent in a small number of infected CD4+ T cells, and these cells constitute a latent reservoir that is the principle barrier to HIV-1 cure. The latent reservoir has a long half-life of ∼44 mo and persists in memory CD4+ T cells, including some that are HIV-1, CMV, and influenza specific.

2A significant fraction of the circulating latent reservoir is composed of expanded clones of CD4+ T cells containing replication-competent proviruses (≥50%;;;;;;; ). Although the origin of the clones and the mechanisms that govern their expansion is not known, longitudinal analysis indicates that they are dynamic and change in size over time in individuals who maintain viral suppression on antiretroviral therapy (ART;;; ).This dynamic may partially account for the longevity of the reservoir. Thus, understanding the basis for latently infected T cell clonal expansion is important for learning how to control and potentially eliminate the reservoir.

3HIV-1 proviral DNA is enriched in HIV-1–, CMV-, and influenza-responsive T cells obtained from ART-suppressed individuals, but whether or how this might be related to clonal expansion of T cells harboring latent viruses that remain replication competent has not been examined. Here, we report that CD4+ T cells containing clones of replication-competent viruses respond to antigenic stimulation with peptides derived from viruses that cause chronic or recurrent infections.

4To test the hypothesis that expanded clones harboring latent proviruses respond to foreign antigens, we exposed CD4+ T cells from ART-suppressed individuals (;; NCT03571204; Table S1 and Table S2) to overlapping peptide pools from common viral and bacterial antigens including HIV group specific antigen (HIV-gag), CMV phosphoprotein 65 (CMV-pp65), or pooled peptides from CMV, EBV, influenza, and tetanus toxin (CEFT). Some of these antigens have been shown to induce HIV-RNA transcription in vivo after vaccination. Staphylococcal enterotoxin B (SEB) and a self-protein, myelin oligodendrocyte glycoprotein (MOG), served as positive and negative controls, respectively, for T cell activation. After overnight culture with HIV-gag, CMV-pp65, CEFT, or SEB, activated CD4+ T cells from eight donors were purified by cell sorting based on expression of two or more activation-induced markers (AIMs; CD69 and PD-L1 or 4-1BB;;;; Fig. 1, a and b; and Fig. S1). Total live CD4+ T cells were sorted from parallel cultures stimulated with MOG to serve as unfractionated controls that were subjected to the same processing conditions. As expected, there was little detectable response to the MOG self-antigen peptide pool, and all donors showed high-level responses to SEB. In addition, responses to HIV-gag, CMV-pp65, and CEFT varied in magnitude among individuals (Figs. 1 c and S1).

5To determine whether CD4+ T cells harboring HIV-1 proviruses are enriched among antigen- or SEB-responsive cells compared with the MOG control, integrated proviruses were enumerated and further characterized as intact or defective by combining quantitative PCR (qPCR) and next-generation sequencing (quadraplex qPCR [Q4PCR];; Fig. 2 a). The combination of PCR and sequencing methods enables initial selection of likely-intact viruses by PCR, and these are subsequently confirmed as intact or defective by near-full-length (NFL) proviral amplification and sequencing. As expected, we observed variation in the ratio of intact to defective proviruses between individuals (; Fig. 2 a).

6The overall frequency of intact and defective proviruses contained within antigen-responsive cells varied substantially among individuals. For example, proviruses were nearly absent from CMV-pp65–specific CD4+ T cells in participant 603 and absent from HIV-gag–specific CD4+ T cells in participants 605 or and 5106. In contrast, 1 in 1,000 HIV-gag–responsive cells harbored a provirus in 603 (Fig. 2 and Table S3). Individuals 603 and 9247 showed proviral enrichment in CD4+ T cells responding to gag of 2- and 4-fold, respectively, whereas in 605 and 5106, proviruses were 3- and 10-fold enriched in pp65-responsive CD4+ T cells (Fig. 2 a). Finally, antigen-responsive cells were relatively depleted of integrated proviruses in participants 9241 and 5104 (Fig. 2 a). We conclude that HIV-1 proviruses can persist in CD4+ T cells that respond to CMV-pp65 and HIV-gag antigen with strong variation among infected individuals.

7We analyzed all HIV-1 sequences across all groups to determine the fraction of viral sequences contributing to expanded clones of CD4+ T cells. Sequences found more than once were classified as derived from clones, while unique sequences were those identified only once. As expected, clones of viral sequences were identified in all participants (Fig. 2 b). The clonal distribution of HIV-1 sequences in antigen-responsive (HIV-1-gag, CMV-pp65, and CEFT) cells was significantly different from the MOG control in five of the six individuals from whom we were able to obtain >10 sequences (Fig. 2 c). The clonality was not due to T cell division in vitro, as there was no measurable T cell division in 18 h under our culture conditions (Fig. S2). We conclude that cells responding to HIV-gag and CMV-pp65 peptides contain clonal proviruses.

8To determine how unique and clonal proviruses contributed to the overall enrichment observed among total HIV-1 sequences, we analyzed the sequences obtained from defective proviruses. Unique sequences were enriched in CMV-pp65–responsive cells from participants 605, 9255, 5104, and 5106 compared with total live CD4+ T cells in the MOG control (Fig. 3 a). There was also decreased representation of unique sequences in response to some antigens in 9247 and 9241 (Fig. 3 a). Otherwise, no distinct pattern emerged from the analysis of defective single viruses, and there was no significant antigen-dependent enrichment among unique sequences between the other individuals.

9Conversely, compared with total live CD4+ T cells in the MOG control, seven of the eight individuals tested showed clones of identical defective proviruses in HIV-gag–, CMV-pp65–, or CEFT-responsive CD4+ T cells (Fig. 3 b). The overall enrichment of defective clonal sequences among antigen-responsive cells was often attributable to one or more specific clones. For example, in individuals 603 and 9247, the relative enrichment of defective proviruses corresponds to two clones isolated from HIV-gag–responsive CD4+ T cells, whereas two defective clones in 605 and 5104 and four defective clones in B207 account for the HIV-1 enrichment in CMV-pp65–responsive CD4+ T cells (Fig. 3 b). In addition, participant 9241 showed an enrichment of two proviral clones in CEFT-responsive cells (Fig. 3 b). Because defective viruses cannot infect multiple individual cells, identical defective viruses can arise due to cellular proliferation. We conclude that expanded CD4+ T cell clones harboring defective HIV-1 proviruses are responsive to peptides derived from HIV-gag, CMV-pp65, and CEFT antigens.

10SEB is a superantigen that stimulates T cells expressing only a subset of T cell receptors. Consistent with this limited specificity, we found that four of the eight individuals tested (9255, B207, 9247, and 5104) had expanded or novel clones of defective proviruses in SEB-reactive CD4+ T cells, compared with the MOG control (Fig. 3 b).

11Similar analysis was also performed using intact proviruses, which, as expected, were present in much smaller numbers than defective proviruses (Fig. 4;;;;; ). Unique proviruses were enriched in HIV-gag–specific CD4+ T cells in participant B207, while in 605, they were found in CMV-pp65–responsive cells (Fig. 4 a).

12Antigen- or SEB-responsive CD4+ T cell clones harboring intact proviruses were identified in three and four of the eight participants, respectively (Fig. 4 b). In B207 and 9247, we found clones of CD4+ T cells responding to CMV-pp65 and HIV-gag, respectively. In 9247, this activity was due to a single intact clone of identical proviruses. Although intact proviral clones were found in total CD4+ T cells in the MOG control in participants 9255, 9241, 5104, and 5106, they were nearly entirely absent from antigen-responsive cells (Fig. 4 b). However, the absence of such sequences might be due to the limited number of cells and panel of antigens assayed. We conclude that cells harboring intact proviruses can respond to peptides from HIV-1 and CMV antigens.

13To determine whether the intact sequences isolated from antigen-responsive cells are replication competent, we compared the sequences obtained from antigen- or SEB-responsive cells to those obtained from single-cell viral outgrowth cultures (quantitative and qualitative viral outgrowth assay [Q2VOA]) from five of the individuals for whom such data were available. We found identical matches between outgrowth cultures and proviral sequences isolated from either CMV-pp65– or HIV-gag–responsive cells in B207 and 9247 (Fig. 5). Additionally, in participant B207, the intact viruses obtained from the largest CMV-pp65–responsive clone were identical to the clone isolated by latency capture. Each member of this clone of CD4+ T cells expressed identical T cell receptor β and α chains (TRBV7-8 and TRAV-21; ). Identical matches for proviruses harbored in SEB-reactive CD4+ T cells were also found in 9255 and B207. Finally, closely related sequences were found in the remaining individuals (Fig. 5). We conclude that clones of CD4+ T cells responsive to HIV-gag, CMV-pp65, and SEB harbor replication-competent viruses.

14Understanding the mechanism of latent reservoir persistence is critical to developing methods for HIV-1 eradication or functional cure. We have used a limited number of antigens to test the idea that CD4+ T cells harboring identical latent proviruses can undergo clonal expansion in response to antigenic stimulation. Our results suggest a mechanism for viral persistence whereby clones of infected cells harboring latent proviruses are stimulated to divide during immune responses.

15A number of ideas have been proposed to explain HIV-1 persistence. These include low-level ongoing HIV-1 replication in sanctuary sites, where drug concentrations are insufficient to block virus replication; and infected T cell longevity. Cell division was not initially considered as a mechanism of proviral persistence because the signals that induce T cell division, such as NF-κB, also tend to reactivate latent proviruses leading to cell death. However, several lines of evidence support the idea that CD4+ T cells harboring latent proviruses can undergo clonal expansion. Initial indirect support for clonal expansion came from longitudinal human studies showing that ART suppressed individuals produce closely related viruses at low levels over time. Direct evidence for the existence of expanded clones was obtained from independent cultures of latent CD4+ T cells in limiting dilution viral outgrowth assays, suggesting that ≥50% of the circulating reservoir consists of expanded clones. In addition, proviral integration site sequencing revealed collections of cells that share identical proviral integration sites. Furthermore, in vitro studies showed that cell division and latency are not mutually exclusive. Finally, isolation of productively infected cells and paired TCR and full-length viral sequencing and paired full-length viral sequencing and integration site analysis provided definitive evidence for the existence of expanded clones of CD4+ T cells harboring identical intact latent proviruses.

16Despite compelling evidence for CD4+ T cell proliferation as a mechanism for HIV-1 persistence, how latently infected cells divide without succumbing to cytopathic viral reactivation remains unknown. Possible explanations include expression of genes that favor survival by suppressing viral transcription, viral integration in transcriptionally inactive regions of the genome, and expression of antiapoptotic proteins such as BIRC5.

17Although there is general agreement that a large fraction of the HIV-1 reservoir comprises clones of CD4+ T cells that wax and wane, far less is understood about what triggers clonal expansion and maintains reservoir longevity. It has been proposed that proviral integration in the vicinity of genes that control cell division, such as cancer-associated genes, promotes cell growth. However, HIV-1 is preferentially integrated into highly transcribed genes that include many cancer-associated genes. Thus, it has been difficult to definitively determine whether or how integration in the vicinity of cancer-related genes contributes to HIV-1 persistence. Moreover, unlike transforming retroviruses that integrate into cancer genes and cause unrestricted cell growth, HIV-1 is not known to cause T cell cancers by integration.

18Our data may provide mechanistic insights into the finding that immunization with tetanus toxoid or influenza can in some, but not all cases, produce increases in HIV-1 RNA in plasma. Individuals with vaccine-responsive T cell clones that harbor intact latent viruses would be expected to produce variable amounts of virus after vaccination. The data are also consistent with the report that viral load blips are more frequent during winter months, when there is a higher incidence of seasonal viral infections.

19Despite the use of a very limited number of test antigens, our work suggests that expanded clones of CD4+ T cells containing replication-competent proviruses respond to pathogens. Their intermittent exposure to these and other antigens found in the virome and microbiome may account for the reported waxing and waning of individual clones of latently infected CD4+ T cells and their persistence over time. Finally, barrier tissues that are chronically exposed to antigen, such as the skin or gut, may represent reservoir sites that contribute to persistence of the reservoir by supporting ongoing expansion and contraction of CD4+ T cell clones harboring latent proviruses.

20All study participants were recruited at the Rockefeller University Hospital, New York, NY, and the University Hospital Cologne, Cologne, Germany (;; and https://clinicaltrials.gov/ct2/show/NCT03571204). All participants provided written informed consent before participation in the studies. The studies were conducted in accordance with Good Clinical Practice. The protocols were approved by the institutional review boards at the Rockefeller University and the University of Cologne. Peripheral blood mononuclear cells (PBMCs) were isolated from leukapheresis by Ficoll separation and frozen in aliquots. All participants were confirmed to be aviremic at the time of sample collection. The samples used in the study described here were collected before any investigational therapeutic intervention.

21Total CD4+ T cells from baseline leukapheresis were isolated from cryopreserved PBMCs by manual magnetic labeling and negative selection using the CD4+ T Cell Isolation Kit (Miltenyi Biotec).

22PBMCs were thawed and washed, and CD8+ T cells were depleted using human CD8 MicroBeads (Miltenyi Biotec, 130-045-201). CD8+ T cell–depleted PBMCs were cultured in 24-well plates at a concentration of 10 × 106 cells/ml in RPMI 1640 supplemented with Hepes, penicillin/streptomycin, and 10% human serum (Sigma-Aldrich). Cells were rested for 3 h and then stimulated with 0.5 µg/ml of HIV-1 consensus B Gag pool, CMV peptide pool, and CEFT peptide pool (PM-HIV-CONB, PM-PP65-2, PM-CEFT-MHC-II, all from JPT Peptide Technologies) for 18 h at 37°C, 5% CO2. A MOG peptide pool (0.5 µg/ml, PM-MOG, JPT Peptide Technologies) served as negative control. SEB-stimulated cells (0.5 µg/ml, Toxin Technology, BT202) served as positive control. Cells were stained for viability with Aquavivid (Life Technologies) at 1/200 in PBS for 20 min, 4°C; incubated with human FcR Block (Miltenyi Biotec, 130-059-901) at 1/70 in PBS-FBS 1% for 10 min, 4°C; and stained with antibodies to surface markers for 30 min, 4°C with Brilliant Violet Buffer (BD Biosciences, 566349) at 1/4 in PBS-FBS 1% (see Table S4 for antibody staining panel).

23AIM+ cells were gated on live, single cells, CD14−CD19−CD8−CD4+, and were sorted based on expression of activation markers CD69 and PD-L1 and/or 4-1BB, as previously described. This combination of markers was chosen because of its specificity for activated memory CD4+ T cells with limited background binding by other T cell subsets. MOG- and SEB-treated cells served as negative and positive controls, respectively, to set gates for sorting. MOG-stimulated cells were sorted based on CD14−CD19−CD8−CD4+ expression. CD69 single-positive cells and AIM− cells (CD69−, 4-1BB−, PD-L1−) were also sorted as controls. Cells were sorted on a FACS Aria II flow cytometer (BD Biosciences) into tubes containing RPMI 1640 supplemented with 10% FBS, Hepes and penicillin/streptomycin. Cells were pelleted, and cell pellets were flash frozen on dry ice and subsequently processed for DNA extraction as described below.

24Flow cytometric data were analyzed using FlowJo version 10.6.0 for Mac. R programming language was used to apply Fisher’s exact test to evaluate whether there was a statistically significant change in the overall distribution of the clones between the Q4PCR-derived HIV-1 sequences from total CD4+ T cells (MOG control) and all antigen-reactive AIM+ CD4+ T cells.

25Genomic DNA from sorted CD4+ T cells was isolated using phenol-chloroform. Briefly, CD4+ T cells were lysed in proteinase K buffer (100 mM Tris, pH 8, 0.2% SDS, 200 mM NaCl, and 5 mM EDTA) and 20 mg/ml proteinase K at 56°C for 12 h, followed by genomic DNA extraction with phenol/chloroform/isoamyl alcohol (Invitrogen), nonfluorescent pellet paint (Millipore, 70748), and ethanol precipitation. DNA concentration was measured by Qubit High Sensitivity Kit (Thermo Fisher Scientific).

26To get a measurement of cells containing gag+ proviruses, genomic DNA was serially diluted and assayed by gag qPCR (concentrations ranged from 600 to 12.5 CD4+ T cells per well). Based on the Poisson distribution, when <30% of qPCR reactions are positive, each positive PCR reaction has a >80% probability of containing a single copy of HIV-1 DNA. DNA was assayed in a minimum of 16 reactions per concentration in a 384-well plate format using the Applied Biosystem QuantStudio 6 Flex Real-Time PCR System. HIV-1–specific primers and a probe targeting a conserved region in gag were used in a limiting dilution qPCR reaction (forward primer, 5′-ATG​TTT​TCA​GCA​TTA​TCA​GAA​GGA-3′; internal probe, 5′-/6-FAM/CCACCCCAC/ZEN/AAGATTTAAACACCATGCTAA-3′/IABkFQ; reverse primer, 5′-TGC​TTG​ATG​TCC​CCC​CAC​T-3′; Integrated DNA Technologies; ).

27Each qPCR reaction was performed in a 10-µl total reaction volume containing 5 µl of TaqMan Universal PCR Master Mix containing Rox (4304437; Applied Biosystems), 1 µl of diluted genomic DNA, nuclease-free water, and 337.5 nM of forward and reverse primers with 93.75 nm of gag internal probe. gag qPCR conditions were 94°C for 10 min, 50 cycles of 94°C for 15 s, and 60°C for 60 s.

28We used a two-step nested PCR approach to amplify NFL HIV-1 genomes. All reactions were performed in a 20-µl reaction volume using Platinum Taq High Fidelity polymerase (Thermo Fisher Scientific). The outer PCR reaction (NFL1) was performed on genomic DNA at a single-copy dilution (previously determined by gag limiting dilution) using outer PCR primers BLOuterF (5′-AAA​TCT​CTA​GCA​GTG​GCG​CCC​GAA​CAG-3′) and BLOuterR (5′-TGA​GGG​ATC​TCT​AGT​TAC​CAG​AGT​C-3′). Touchdown cycling conditions were 94°C for 2 min and then 94°C for 30 s, 64°C for 30 s, and 68°C for 10 min for three cycles; 94°C for 30 s, 61°C for 30 s, and 68°C for 10 min for three cycles; 94°C for 30 s, 58°C for 30 s, and 68°C for 10 min for three cycles; 94°C for 30 s, 55°C for 30 s, and 68°C for 10 min for 41 cycles; and then 68°C for 10 min.

29Q4PCR was performed as previously described. Briefly, undiluted 1-µl aliquots of the NFL1 PCR product were subjected to a Q4PCR reaction using a combination of four primer/probe sets that target conserved regions in the HIV-1 genome. Each primer/probe set consists of a forward and reverse primer pair as well as a fluorescently labeled internal hydrolysis probe as follows: PS forward, 5′-TCT​CTC​GAC​GCA​GGA​CTC-3′; reverse, 5′-TCT​AGC​CTC​CGC​TAG​TCA​AA-3′; probe, 5′-/Cy5/TTTGGCGTA/TAO/CTCACCAGTCGCC-3′/IAbRQSp (Integrated DNA Technologies; ); env forward, 5′-AGT​GGT​GCA​GAG​AGA​AAA​AAG​AGC-3′; reverse, 5′-GTC​TGG​CCT​GTA​CCG​TCA​GC-3′; probe, 5′-/VIC/CCTTGGGTTCTTGGGA-3′/MGB (Thermo Fisher Scientific; ); gag forward, 5′-ATG​TTT​TCA​GCA​TTA​TCA​GAA​GGA-3′; reverse, 5′-TGC​TTG​ATG​TCC​CCC​CAC​T-3′; probe, 5′-/6-FAM/CCACCCCAC/ZEN/AAGATTTAAACACCATGCTAA-3′/IABkFQ (Integrated DNA Technologies; ); and pol forward, 5′-GCA​CTT​TAA​ATT​TTC​CCA​TTA​GTC​CTA-3′; reverse, 5′-CAA​ATT​TCT​ACT​AAT​GCT​TTT​ATT​TTT​TC-3′; probe, 5′-/NED/AAGCCAGGAATGGATGGCC-3′/MGB (Thermo Fisher Scientific; ).

30Each Q4PCR reaction was performed in a 10-µl total reaction volume containing 5 µl TaqMan Universal PCR Master Mix containing Rox (4304437; Applied Biosystems), 1 µl diluted genomic DNA, nuclease-free water, and the following primer and probe concentrations: PS, 675 nM of forward and reverse primers with 187.5 nM of PS internal probe; env, 90 nM of forward and reverse primers with 25 nM of env internal probe; gag, 337.5 nM of forward and reverse primers with 93.75 nM of gag internal probe; and pol, 675 nM of forward and reverse primers with 187.5 nM of pol internal probe. qPCR conditions were 94°C for 10 min, 40 cycles of 94°C for 15 s, and 60°C for 60 s. All qPCR reactions were performed in a 384-well plate format using the Applied Biosystem QuantStudio 6 Flex Real-Time PCR System.

31We used QuantStudio Real-Time PCR Software version 1.3 (Thermo Fisher Scientific) for data analysis. The same baseline correction (start cycle, 3; end cycle, 10) and normalized reporter signal (ΔRn) threshold (ΔRn threshold = 0.025) was set manually for all targets/probes. Fluorescent signal above the threshold was used to determine the threshold cycle. Samples with a threshold cycle value between 10 and 40 of any probe or probe combination were identified. Samples showing reactivity with two or more of the four qPCR probes were selected for further processing.

32The nested PCR reaction (NFL2) was performed on undiluted 1-µl aliquots of the NFL1 PCR product. Reactions were performed in a 20-µl reaction volume using Platinum Taq High Fidelity polymerase (Thermo Fisher Scientific) and PCR primers 275F (5′-ACA​GGG​ACC​TGA​AAG​CGA​AAG-3′) and 280R (5′-CTA​GTT​ACC​AGA​GTC​ACA​CAA​CAG​ACG-3′; ) at a concentration of 800 nM. Touchdown cycling conditions were 94°C for 2 min and then 94°C for 30 s, 64°C for 30 s, and 68°C for 10 min for three cycles; 94°C for 30 s, 61°C for 30 s, and 68°C for 10 min for three cycles; 94°C for 30 s, 58°C for 30 s, and 68°C for 10 min for three cycles; 94°C for 30 s, 55°C for 30 s, and 68°C for 10 min for 41 cycles; and then 68°C for 10 min.

33All nested PCR products from NFL2 were subjected to library preparation. The Qubit 3.0 Fluorometer and Qubit dsDNA BR Assay Kit (Thermo Fisher Scientific) were used to measure DNA concentrations. Samples were diluted to a concentration of 10–20 ng/µl. Tagmentation reactions were performed using 1 µl of diluted DNA, 0.25 µl Nextera TDE1 Tagment DNA enzyme (15027865), and 1.25 µl TD Tagment DNA buffer (15027866; Illumina). Tagmented DNA was ligated to unique i5/i7 barcoded primer combinations using the Illumina Nextera XT Index Kit v2 and KAPA HiFi HotStart ReadyMix (2X; KAPA Biosystems) and then purified using AmPure Beads XP (Agencourt). 384 purified samples were pooled into one library and then subjected to paired-end sequencing using Illumina MiSeq Nano 300 V2 cycle kits (Illumina) at a concentration of 12 pM.

34HIV-1 sequence assembly was performed by our in-house pipeline (Defective and Intact HIV Genome Assembler), which is capable of reconstructing thousands of HIV genomes within hours via the assembly of raw sequencing reads into annotated HIV genomes. The steps executed by the pipeline are described briefly as follows. First, we removed PCR amplification and performed error correction using clumpify.sh from BBtools package v38.72 (https://sourceforge.net/projects/bbmap/). A quality control check was performed with Trim Galore package v0.6.4 (https://github.com/FelixKrueger/TrimGalore) to trim Illumina adapters and low-quality bases. We also used bbduk.sh from BBtools package to remove possible contaminant reads using HIV genome sequences, obtained from Los Alamos HIV database, as a positive control. We used a k-mer–based assembler, SPAdes v3.13.1, to reconstruct the HIV-1 sequences. The longest assembled contig was aligned via BLAST to a database of HIV genome sequences, obtained from Los Alamos, to set the correct orientation. Finally, the HIV genome sequence was annotated by aligning against Hxb2 using BLAST. Sequences with double peaks, i.e., regions indicating the presence of two or more viruses in the sample (cutoff consensus identity for any residue <70%), or samples with a limited number of reads (empty wells ≤500 sequencing reads) were omitted from downstream analyses. In the end, sequences were classified as intact or defective. Defective sequences were subdivided into more specific classifications according to their sequence structure: Major Splice Donor (MSD) Mutation, Non-functional, or Missing Internal Genes.

35Clones were defined by aligning sequences of each classification (Intact, MSD Mutation, Non-functional, and Missing Internal Genes) to HXB2 and calculating their Hamming distance. Sequences having a maximum of three differences between the first nucleotide of gag and the last nucleotide of nef (reference: HXB2) were considered members of the same clone if found more than once across sequences from baseline or all sorted populations (see Cell sorting) retrieved from each participant.

36All data are available via National Center for Biotechnology Information GenBank under accession nos. MN090896–MN090943 and MT189273–MT191225.

37Fig. S1 shows the AIM assay flow cytometry plots from all participants. Fig. S2 shows the CFSE proliferation plots after 18 h of stimulation in the AIM assay. Table S1 shows the participants’ demographics. Table S2 shows the HLA typing of the study participants. Table S3 shows the gag limiting dilution calculations used for the Q4PCR assay. Table S4 shows the list of antibodies used in flow cytometry for the AIM assay.