INTRODUCTION
The CD8 T cell response to human/simian immunodeficiency virus (HIV/SIV) is an aggregate of multiple, discrete epitope-specific responses. Patterns of T cell immunodominance are often used to categorize the antiviral efficacy of these individual CD8 T cell responses. At any given time, the most immunodominant CD8 T cell response is the highest frequency epitope-specific CD8 T cell response, while other so-called “subdominant” CD8 T cell responses are present at lower individual frequencies (
1). The earliest immunodominant CD8 T cell responses contribute to the establishment of the viral load set point and define the likelihood of long-term viral control within an individual (
2,
3), while the subdominant CD8 T cell responses contribute to effective control of chronic HIV/SIV replication (
4,
5). Alternatively, the antiviral efficacy of CD8 T cells can be defined by whether they select for immune escape variants rapidly or slowly. Epitopes that do not rapidly accumulate variants have become attractive as vaccine immunogens because they are common in circulating strains of HIV (
6). Nonetheless, it remains unclear to what extent CD8 T cell responses specific for epitopes that rapidly accumulate variants and those specific for epitopes that remain invariant can control early viral replication. Knowledge of which populations of CD8 T cells can best suppress viral replication is needed to design a T cell-based vaccine.
It has been difficult to define whether acute-phase CD8 T cell responses targeting invariant epitopes can establish control of HIV replication. The development of individual epitope-specific CD8 T cell responses is dependent on both the total repertoire of HLA/major histocompatibility complex (MHC) class I alleles in the host and the sequence of the infecting virus. A “pre-escaped” HIV can be transmitted between individuals who share a single HLA class I allele, but acute-phase CD8 T cell responses that develop in the newly infected individual can still select for escape variants in epitopes that did not accumulate variants in the donor (
7–9). Although the efficacy of acute-phase CD8 T cell responses targeting invariant epitopes could be measured when a pre-escaped HIV is transmitted between fully HLA-identical people, this event is exceptionally rare. Thus, it is nearly impossible to evaluate the
in vivo antiviral efficacy of acute-phase CD8 T cell responses that do not rapidly select for variants in HIV-positive individuals.
Nonhuman primate studies with engineered SIV strains containing point mutations in T cell epitopes can be used to model the transmission of pre-escaped HIV. These studies typically focus on responses restricted by a single MHC class I allele, but CD8 T cell responses targeting epitopes restricted by alternate MHC class I alleles can develop and rapidly select for escape variants. Not surprisingly, control of pre-escaped SIV replication has been inconsistent and unpredictable in these studies (
10–12) and the antiviral efficacy of acute-phase CD8 T cell responses targeting invariant epitopes could not be evaluated. Fully MHC class I-matched Mauritian cynomolgus macaques (MCMs) infected with SIV can overcome these limitations by mounting consistent CD8 T cell responses that select for sequence variants in the same epitopes (
13,
14). MCMs that are homozygous for the M3 MHC class I haplotype are present at a frequency of ∼4% (
15). There are now 14 known epitopes restricted by MHC class I molecules encoded by the M3 haplotype (
13,
16–18), of which 5 consistently accumulate sequence variants within the first 12 weeks after SIVmac239 infection in M3 homozygous MCMs. With such exquisite knowledge of host MHC genetics and epitope-specific CD8 T cell immunity, this model can be customized to evaluate the antiviral efficacy of CD8 T cell responses that do not rapidly select for escape variants.
We took advantage of the unique replication kinetics of the live attenuated virus SIVmac239Δnef to understand the contribution of CD8 T cell responses that do not rapidly select for escape variants on elite viral control. In Indian rhesus macaques and Mauritian cynomolgus macaques infected with SIVmac239Δnef, an acute-phase peak of virus replication is subsequently controlled to levels that are not detectable with plasma viral load assays (
19,
20). Viral load decline is coincident with the emergence of CD8 T cell responses (
21,
22) and can be used as an indicator of T cell function.
We created an SIVmac239Δnef derivative, termed m3KOΔnef, in which we ablated four of the five epitopes that accumulate variants within 12 weeks of SIVmac239 infection. The fifth epitope is in an alternate reading frame and was undefined when we created m3KOΔnef, so it remained intact. We eliminated one additional well-characterized epitope in Pol that frequently accumulates variants by 14 weeks after infection (
18). M3 homozygous MCMs in the acute phase of infection with m3KOΔnef are therefore capable of developing acute-phase CD8 T cell responses that target less variable epitopes, making this the first model available to assess the functional potency of these CD8 T cell responses to control early SIV replication. We found that M3 homozygous MCMs could not control m3KOΔnef as well as MCMs that did not have any of the MHC class I alleles present in the M3 MHC haplotype. These results provide compelling evidence that CD8 T cell responses that do not rapidly select for escape variants have limited potency
in vivo and argue that caution should be used when designing vaccines to elicit T cell responses that target conserved epitopes.
MATERIALS AND METHODS
Creation of virus stocks.
The following clonal virus stocks were created for this study: SIVmac239Δnef, barcoded SIVmac239Δnef (BCVΔnef), and m3KOΔnef. For SIVmac239Δnef, we obtained the necessary plasmids containing the 5′ and 3′ viral genomes (p239SpSp5′ and pSIVmac239Δnef deletion mutant) from Ronald Desrosiers through the AIDS Research and Reference Reagent Program, Division of AIDS, NIAID, NIH. The plasmid containing the 5′ SIV genome with 10 synonymous “barcoding” mutations in Gag was used previously (
12). For m3KOΔnef, we first identified consensus variants that accumulated in virus populations replicating in M3
+ MCMs during chronic SIVmac239 infection (
15; data not shown). Twenty-three substitutions were then incorporated into SIVmac239Δnef by custom gene synthesis (GenScript, Piscataway, NJ) to create the 5′ and 3′ plasmids necessary to construct m3KOΔnef. These plasmids were designed to include the SphI restriction site required for coligation, and they contained long terminal repeat sequences matching those in SIVmac239.
All virus stocks were created as described previously (
23,
24). The plasmids containing the 5′ and 3′ halves of the corresponding genomes were digested with SphI, treated with shrimp alkaline phosphatase, precipitated, and then ligated together. Vero cells were transfected with the ligated products and then cocultured with CEMx174 cells for 48 h. Infected CEMx174 cells were grown for about 1 week to produce high-titer viruses. The p27 content of each virus stock was determined by enzyme-linked immunosorbent assay (ELISA; ZeptoMetrix, Buffalo, NY) according to the manufacturer's protocol. The infectious titer of each virus stock was determined by limiting-dilution inoculation of CEMx174 cells. Briefly, cells were inoculated in quadruplicate with serial 4-fold dilutions of virus stock, beginning at a 1:100 dilution. Two weeks later, individual cultures were scored as positive or negative for SIV infection by qualitative p27 antigen capture ELISA (ZeptoMetrix). The infectious titer was then estimated by using the proportionate-distance method of Reed and Muench (
25) and expressed in terms of 50% tissue culture infectious doses (TCID
50) per milliliter of virus stock.
In vitro coculture competition assays.
All virus stocks were normalized for a p27 content of 50 ng/ml. Viruses were mixed in the following combinations: 2.5 ng m3KOΔnef or SIVmac239Δnef to 2.5 ng BCVΔnef (1:1), 0.5 ng m3KOΔnef or SIVmac239Δnef to 4.5 ng BCVΔnef (1:9), and 4.5 ng m3KOΔnef or SIVmac239Δnef to 0.5 ng BCVΔnef (9:1). Each virus mixture was incubated with 106 CEMx174 cells at 37°C for 4 h. After washing, 5 × 105 cells were plated and grown for 1 week. Supernatant was collected at days 3, 5, and 7.
The copy number of each virus was measured in the inoculum and in each supernatant. Viral RNA was isolated with the M48 Biorobot apparatus (Qiagen, Valencia, CA). Quantification was conducted with the SuperScript III Platinum One-Step quantitative PCR kit (Invitrogen, Carlsbad, CA) as previously described (
12). The SIVmac239Δnef and m3KOΔnef viruses were quantified with primers and probes targeting an 84-bp region of
gag (
12,
26). BCVΔnef was quantified in a separate reaction with a distinct set of primers and probes (
12). The appropriate wild-type or barcoded transcript was used for each standard curve. All comparative reactions were run on the same plate with a Roche LightCycler 480 under the following conditions: cDNA synthesis at 50°C for 30 min, an initial denaturation step at 95°C for 2 min, and 45 amplification cycles of 95°C for 10 s and 69°C for 30 s.
The ratio of m3KOΔnef or SIVmac239Δnef to BCVΔnef was determined for each sample and then normalized to the ratio that was present in the inoculum. Replicative differences between viruses were assessed at each time point with unpaired t tests (GraphPad Prism, La Jolla, CA). All data were plotted on a log2 scale.
Animal care and use.
MCMs were purchased from Charles River Laboratories and cared for by the Wisconsin National Primate Research Center according to protocols approved by the University of Wisconsin Graduate School Animal Care and Use Committee. Animals CY0348, CY0379, CY0381, CY0382, CY0383, CY0384, and CY0385 were infected intravenously with 10 ng of p27 m3KOΔnef. Animals CY0332, CY0333, CY0334, CY0335, CY0336, and CY0337 were infected intrarectally with 7,000 TCID
50 of SIVmac239 as part of another study (
13,
16). Animals CY0114, CY0157, CY0160, CY0205, CY0206, CY0213, CY0338, CY0345, and CY0386 were infected intravenously with 10 ng of p27 SIVmac239Δnef (
19).
Plasma viral loads analysis.
Plasma SIV loads were determined essentially as previously described (
15,
26). Viral RNA was isolated from plasma and then reverse transcribed and amplified with the SuperScript III Platinum one-step quantitative reverse transcription (RT)-PCR system (Invitrogen). Samples were quantified on a Roche LightCycler 2.0 and compared to an internal standard curve on each run. The substitutions present in m3KOΔnef did not lie within the primer or probe target sequence.
Differences in the log10 viral load set point were modeled by linear mixed-effects regression with random effects for each animal and fixed effects for each group. The mean log10 viral load set point for a group was estimated by viral load measurements between weeks 14 and 20 after infection. The maximum log10 viral load before 5 weeks was defined as the acute-phase peak viral load, and differences were modeled by linear regression. Statistics and significance measures were calculated in R 2.15.2 with the lmer and gmodels packages.
Genome-wide sequencing of SIV and data analysis.
Genome-wide sequencing of replicating SIV was performed essentially as previously described (
13,
27). Viral RNA was isolated from plasma with the QIAamp MinElute Virus Spin kit (Qiagen). The Superscript III One-Step RT-PCR system with Platinum
Taq High Fidelity (Invitrogen) was used to reverse transcribe the viral RNA in four overlapping amplicons spanning the entire SIV genome. PCR products were purified with the Qiagen MinElute Gel Extraction kit (Qiagen) and quantified with the Quant-IT dsDNA HS Assay kit (Invitrogen).
For pyrosequencing, libraries were prepared with the Nextera DNA Sample Prep kit (Roche Titanium compatible; Epicentre, Madison, WI) and 10-bp multiplex identifier tags were added. Tagged products were purified twice with Agencourt AMPure XP beads (Beckman Coulter Genomics, Danvers, MA) and then quantified with the Quant-IT dsDNA HS Assay kit (Invitrogen) and the Agilent High Sensitivity DNA kit (Agilent Technologies, Santa Clara, CA). Pyrosequencing was performed with a Roche/454 GS Junior instrument with Titanium shotgun chemistry, according to the manufacturer's protocols (454 Life Sciences, Branford, CT).
Nucleotide variation relative to SIVmac239 was determined at each position with Galaxy, an online tool used to analyze next-generation sequence data (
13,
28–30). FASTQ reads were trimmed, and low-quality uncertain sequences were masked with an “N.” Sequence reads were mapped to an SIVmac239 reference (accession no. M33262) with the LASTZ algorithm. The SAM format was converted to a BAM format, and a pileup was created. The total A, C, T, G, and N numbers at each site were calculated. Sites that were greater than 1% variant from the reference were assessed as synonymous or nonsynonymous with the SNPeffect tool (
31). Variants spanning nucleotide positions 1300 to 10200 with >1% variation and greater than three times the percentage of uncertain or indel sequences were plotted with GraphPad Prism.
Amino acid variation relative to SIVmac239 was determined as described previously (
13,
16). FASTQ reads were trimmed, and low-quality sequences were masked with an “N.” FASTQ reads were converted to FASTA reads and translated in all six reading frames. BLAT (BLAST-like alignment tool) was used to align sequence reads with SIVmac239 proteins at a minimum of 50% identity. Custom scripts were used to extract the region of interest, remove reads with ambiguous sequences, and measure the frequency of each variant identified.
IFN-γ ELISPOT assays.
Gamma interferon (IFN-γ) enzyme-linked immunospot (ELISPOT) assays were performed as described previously (
15,
18). First, peripheral blood mononuclear cells (PBMCs) were isolated from EDTA-anticoagulated blood by Ficoll-Paque Plus (GE Healthcare, Piscataway, NJ) density gradient centrifugation. A precoated monkey IFN-γ ELISPOTplus plate (Mabtech, Mariemont, OH) was blocked, and individual peptides were added to each well at a final concentration of 1 μM. Each peptide was tested in duplicate, and concanavalin A (10 μM) was used as a positive control. Assays were performed according to the manufacturer's protocol, and wells were imaged with an AID ELISPOT reader. The average number of spots per peptide was calculated, and we subtracted the average number of spots formed with no peptide. Data were extrapolated to 10
6 PBMCs to calculate the number of spot-forming cells (SFCs) per 10
6 PBMCs.
Generation of peptide-specific T cell lines.
T cell lines were grown as described previously (
17,
18). First, 5 × 10
6 PBMCs were isolated from CY0381 and incubated with the immunogenic peptide (either ARF1
29-43VY15 or ARF1
30-40QL11) for 1 week in medium supplemented with 5 × 10
5 U/ml of interleukin-2 (Prometheus Laboratories, San Diego, CA). Restimulation was conducted every 1 to 2 weeks with autologous irradiated B-lymphoblastoid cells pulsed with the relevant peptide.
Intracellular cytokine staining.
Activation of peptide-specific T cells was measured via the production of IFN-γ and tumor necrosis factor alpha (TNF-α) as described previously (
18). Briefly, peptide-pulsed autologous B-lymphoblastoid cells were used to stimulate cultured T cells at a 1:1 ratio for 4 h in the presence of brefeldin A (10 μg/ml; Sigma-Aldrich). Single MHC allele transferent cell lines, pulsed with peptide and washed extensively, were used similarly to determine MHC restriction. Cells were then stained with anti-CD8–PacBlue (BD Biosciences, San Jose, CA) and anti-CD3–AF700 (BD Biosciences) antibodies, fixed with 2% paraformaldehyde, and left at 4°C overnight. On the following morning, the cells were permeabilized in 1% saponin buffer, stained with anti-IFN-γ–fluorescein isothiocyanate (FITC) (BD Biosciences) and anti-TNF-α–phycoerythrin (PE) (BD Biosciences) antibodies, washed in 1% saponin buffer, and fixed in 2% paraformaldehyde. Flow cytometry was performed with an LSRII instrument (BD Biosciences) and data were analyzed with FlowJo V.9.4.10 software (TreeStar Inc., Ashland, OR).
Tetramer staining.
Fluorochrome-conjugated MHC/peptide tetramers were produced as described previously (
32). The tetramers used in this study were Mafa-A1*063/ARF1
30-40QL11-PE, Mafa-A1*063/Env
338-346RF9-PE, Mafa-B*075/Tat
42-49QA8-PE, and Mafa-B*075/Rev
59-68SP10-PE. Thawed lymph node cells or PBMCs were stained with tetramer at 37°C for 20 min, washed, and then stained with anti-CD3–AF700 (BD Biosciences), anti-CD8–PacBlue (BD Biosciences), and anti-CD69–Texas Red–PE (Beckman Coulter) antibodies. Flow cytometry was performed with an LSRII instrument (BD Biosciences), and data were analyzed with FlowJo V.9.4.10 software (TreeStar Inc.). Tetramer-binding CD8 T cell frequencies were compared between groups with paired or unpaired two-tailed
t tests (GraphPad Prism) after log transformation of the data.
Measurement of immune activation.
Immune activation assays were based on previous reports (
33). Whole blood or cells isolated from bronchoalveolar lavage (BAL) fluid were first incubated with surface stains. Blood samples were then treated with FACSlyse, and BAL fluid samples were treated with 2% paraformaldehyde. Cells were then permeabilized in 0.1% saponin buffer, stained with antibodies specific for intracellular markers, washed, and analyzed by flow cytometry. The antibodies used were anti-CD38–FITC (Stem Cell Technologies, Inc.), anti-CD4–PeCy7 (BD Biosciences), anti-CD3–PacBlue (BD Biosciences), anti-Ki-67–AF647 (BD Biosciences), anti-CD20–AF700 (BioLegend), anti-CD8–APCCy7 (BioLegend, San Diego, CA), and anti-CD69–Texas Red–PE (Beckman Coulter, Indianapolis, IN).
DISCUSSION
In the present study, we took advantage of MHC-identical MCMs to study CD8 T cell-mediated control of live attenuated SIV replication, an MHC-independent model of elite control. We created a derivative of SIVmac239Δnef incorporating commonly detected nucleotide substitutions identified in virus populations replicating in M3 homozygous MCMs. In m3KOΔnef, we ablated five epitopes that all accumulate variants by 14 weeks after SIVmac239 infection of M3 homozygous MCMs. Even though CD8 T cell responses targeting these five epitopes were absent, we detected reproducible CD8 T cell responses targeting three additional epitopes that remained intact in m3KOΔnef, two of which typically do not select for escape variants until the chronic infection phase in M3 homozygous MCMs infected with SIVmac239. Replication of m3KOΔnef was not controlled in M3 homozygous animals, suggesting that the set of CD8 T cell responses that predominated in M3 homozygous MCMs during the acute phase of infection with m3KOΔnef and do not tend to rapidly select for escape variants were unable to control viral replication.
The key observation in this study is that control of live attenuated SIV replication is strongly associated with acute-phase CD8 T cell responses that rapidly select for sequence variants. We found that CD8 T cell responses that do not typically select for escape variants during the first few weeks of SIV infection were insufficient to control the replication of m3KOΔnef. In contrast, m3KOΔnef replication was controlled in animals with entirely different MHC class I alleles, enabling the presentation of alternate epitopes. These results suggest that the CD8 T cell responses that typically emerge during the chronic infection phase and do not rapidly select for escape variants are collectively less potent than the acute-phase CD8 T cell responses that do rapidly select for escape variants.
An alternate explanation for these results is that acute-phase CD8 T cell responses target some epitopes that do rapidly accumulate variants and some that do not. By specifically eliminating the responses that do select for escape variants, we inadvertently reduced the effective breadth of the acute-phase CD8 T cell response, thus reducing viral control. We do not favor this interpretation because we have not observed the entire array of CD8 T cell responses in peripheral blood during the acute phase of infection of M3 homozygous MCMs with SIVmac239 (
13,
16); several responses are detectable only later in infection. It is possible that these responses are exerting an antiviral effect that is not detectable by our immunological assays, but this seems unlikely.
We also detected accelerated variation kinetics in the ARF130-40QL11 and Env338-346RF9 epitopes when CD8 T cell responses targeting these epitopes were present at a high frequency at or near the time of peak viral replication. It is reasonable to suggest that the rapid selection of these epitope variants could be attributed to higher viral replication at the time the CD8 T cell responses were present rather than as a function of enhanced epitope-specific CD8 T cell potency. Indeed, this observation suggests that the correlation between CD8 T cell potency and the selection of escape variants is not absolute. Thus, CD8 T cells that select for escape variants are not always the most potent; however, CD8 T cells that are potent do tend to select for sequence variants. In other words, rapid variation is not sufficient by itself to implicate a specific CD8 T cell response in protective immunity, but this approach can be used as an initial filter to identify potent CD8 T cell responses. Further studies are needed to evaluate whether CD8 T cell responses targeting different epitopes throughout infection are functionally different and, consequently, have distinct abilities to destroy virus-infected cells.
There are two important caveats to this study. First, the animal groups used were necessarily small. While we dedicated considerable resources to assembling a cohort of M3 homozygous MCMs to challenge with m3KOΔnef, we were unable to assemble a similarly sized cohort of M3 homozygous MCMs to challenge with SIVmac239Δnef. We now know that M3 homozygous animals generally exhibit high viral loads when challenged with SIVmac239 (
15,
16), so we cannot formally exclude the possibility that MCMs with this genotype are predisposed to higher viral loads when infected with any SIV strain. Nonetheless, the nearly uniform complete control of SIVmac239Δnef replication in both Indian rhesus and Mauritian cynomolgus macaques (
19,
20) (
Fig. 2), irrespective of MHC background, makes this unlikely. In addition to the small group sizes, our m3KOΔnef controls were not MHC identical to one another. The only MCM haplotypes that do not share any MHC class I alleles with the M3 haplotype are rare (
34), so animals that are fully MHC identical for these alternate haplotypes are exceptionally uncommon. Even though the animal numbers were small, we believe that the use of fully MHC-identical, M3 homozygous MCMs infected with m3KOΔnef provides the best cohort used to date to specifically study the
in vivo antiviral efficacy of CD8 T cell responses that do not rapidly select for escape variants.
Second, we recognize that our results are confounded somewhat by the unexpectedly strong response to the cryptic ARF130-40QL11 epitope. On the basis of our classification of early T cell responses that rapidly select for sequence variants, we would ablate this epitope in the pre-escaped virus if we repeated this study. We did not include ARF130-40QL11 epitope variants in m3KOΔnef because we originally identified common variants that were nonsynonymous in the nine main SIV genes; we did not consider cryptic epitopes. Still, the presence of acute-phase CD8 T cells targeting ARF130-40QL11 was insufficient to enable control of m3KOΔnef in M3 homozygous MCMs. It is provocative to consider that control of m3KOΔnef may have been even weaker in the absence of responses targeting ARF130-40QL11.
The limited potency of acute-phase CD8 T cell responses targeting epitopes that do not rapidly accumulate variants calls into question the value of designing HIV vaccines to elicit CD8 T cell responses targeting conserved elements in the viral proteome. Conserved-element vaccines are attractive because the immunogenic sequences targeted by vaccine-specific immune responses should be present in most circulating virus populations (
6,
40–42). Unfortunately, T cell responses targeting conserved or invariant epitopes are also typically subdominant in HIV infection (
43). Our data suggest that these responses are not particularly effective for long-term control of viral replication, but we do not know whether CD8 T cell responses targeting epitopes that are even less variable than the ones studied here could effectively control viral replication. With the MCM model, we can repeat the strategy used to generate m3KOΔnef and serially ablate additional epitopes restricted by M3 MHC alleles. Given enough generations, we expect to elicit acute-phase CD8 T cell responses targeting epitopes that are practically invariant or potentially eliminate the entire set of CD8 T cell responses (
27). Ultimately, challenging M3 homozygous MCMs with a “T cell epitope-deficient virus” could help quantify the requirement for CD8 T cell responses with reference to MHC-independent viral control.
Overall, the observations from this SIV model can be applied to our understanding of the complex interplay between host immune responses and viral evolution at the population level. As HIV continues to adapt to common HLA alleles in the population, the most potent epitopes may not be immunogenic in transmitted HIV strains, such that alternate CD8 T cell responses targeting less variable epitopes will predominate during the acute phase of infection. This might ultimately lead to less efficient viral control and accelerated disease progression, an alarming forecast with significant implications for the future of the HIV epidemic.
ACKNOWLEDGMENTS
This work was supported by a CHAVI/HVTN Early Stage Investigator award and by University of Wisconsin—Madison startup funds to S.L.O. The research was conducted in part at a facility constructed with support from Research Facilities Improvement Program grants RR15459-01 and RR020141-01. This project was funded by the National Institutes of Health (grant R01AI084787); the National Cancer Institute, National Institutes of Health (contract HHSN26120080001E); and the Intramural Research Program of the Center for Cancer Research, National Cancer Institute, National Institutes of Health. We thank members of the Wisconsin National Primate Research Center, a facility supported by grants P51RR000167 and P51OD011106. D.A.P. is a Wellcome Trust (United Kingdom) Senior Investigator. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government.
We thank Lindsey Traeger and Matt Scarlotta for performing some of the sequencing experiments; Andy VanPay for his careful reading of the manuscript; Esper Kallas at the University of Sao Paulo for his valuable discussions; and Justin Greene, David O'Connor, and members of his laboratory at the University of Wisconsin—Madison for providing samples.