Correction: CD8 T Cell Response Maturation Defined by Anentropic Specificity and Repertoire Depth Correlates with SIVΔnef-induced Protection

The live attenuated simian immunodeficiency virus (LASIV) vaccine SIVΔnef is one of the most effective vaccines in inducing protection against wild-type lentiviral challenge, yet little is known about the mechanisms underlying its remarkable protective efficacy. Here, we exploit deep sequencing technology and comprehensive CD8 T cell epitope mapping to deconstruct the CD8 T cell response, to identify the regions of immune pressure and viral escape, and to delineate the effect of epitope escape on the evolution of the CD8 T cell response in SIVΔnef-vaccinated animals. We demonstrate that the initial CD8 T cell response in the acute phase of SIVΔnef infection is mounted predominantly against more variable epitopes, followed by widespread sequence evolution and viral escape. Furthermore, we show that epitope escape expands the CD8 T cell repertoire that targets highly conserved epitopes, defined as anentropic specificity, and generates de novo responses to the escaped epitope variants during the vaccination period. These results correlate SIVΔnef-induced protection with expanded anentropic specificity and increased response depth. Importantly, these findings render SIVΔnef, long the gold standard in HIV/SIV vaccine research, as a proof-of-concept vaccine that highlights the significance of the twin principles of anentropic specificity and repertoire depth in successful vaccine design.


Introduction
Two decades ago, a report described a cohort of rhesus macaques infected with the live attenuated simian immunodeficiency virus (LASIV), SIVΔnef, and subsequently protected from pathogenic wild-type SIV challenge [1]. Since then, numerous studies have demonstrated the efficacy of SIVΔnef-induced protection, ranging from complete protection with sterilizing immunity to partial protection with two or more logs reduction in peak and set-point viremia [2][3][4][5][6]. Remarkably, presumably in part due to the replication of SIVΔnef in mucosal sites [7], robust protection has been also been documented following mucosal challenges [3][4][5]. SIVΔnef has also induced significant protection against heterologous challenge, albeit less effectively than against homologous challenge [4,5]. However, studies describing disease progression in SIVΔnef-vaccinated infant macaques and a subset of adult non-human primates precluded live attenuated HIV from being developed as a vaccine in human subjects [8][9][10].
Due to concerns over safety, research on SIVΔnef and related LASIV vaccines has shifted from safety-and-efficacy determination to mechanism-of-action delineation. As the most effective lentiviral vaccine, SIVΔnef has been extensively studied in order to shed light on the correlates of vaccine-mediated protection. However, no immunological correlate or mode of action has consistently been identified as being responsible for protection against pathogenic challenge. SIVΔnef generates a diverse SIV-specific antibody response [11,12] and macaques vaccinated with the related attenuated virus SIVΔ3 and lacking the Mamu-AÃ01 MHC I allele show effective control of pathogenic viral challenge despite CD8 T cell depletion [13], implying that humoral immunity may play a significant role, at least in some genotypic backgrounds. Innate immunity has also been implicated as a correlate of SIVΔnef-induced protection [3]. SIVΔnef infection induces potent CD8 T cell responses similar in magnitude to wild-type SIV infection [14], and multiple studies have implicated SIV-specific CD8 T cells in the protection induced by SIVΔnef [4][5][6]15]. Most recently, Fukazawa et al. [6] correlated the magnitude of lymph node SIV-specific T cell responses with protection elicited by a range of different LASIV strains, including SIVΔnef. Finally, studies in a related vaccine challenge model involving vaccination with an attenuated SHIV strain have implicated SIV-specific CD8 T cell responses in the female reproductive tract as playing a crucial role in protection [16].
In this work, we sought to scrutinize the interplay of the attenuated virus and the immune response at a high resolution. We chose to focus on the CD8 T cell response as the most likely immune correlate for protection given that previous studies have demonstrated: potent antilentiviral activity of CD8 T cells in vitro [17][18][19][20][21]; increased viremia during chronic SIV infection after CD8 T cell depletion [22]; the ability of a live attenuated SHIV vaccine containing HIV env to protect vaccinated macaques against SIVmac239 challenge [16,23]; SIVΔnefmediated protection of animals challenged with wild-type SIVmac239 containing highly heterologous env sequences [24]; and a correlation between increased LASIV-induced protection and higher magnitude of SIV-specific CD4 and CD8 T cells in the lymph nodes of vaccinated animals [6].
A distinguishing feature of the SIVΔnef and SIVΔ3 vaccines is the gradual increase in protection against wild-type SIV challenge during the first 15-20 weeks of vaccination [2,25]. Given the low-level replication rate of SIVΔnef after initial control of viremia by 8-12 weeks after vaccination, it was unknown whether SIVΔnef undergoes sequence evolution after acute infection. The only previously documented sequence evolution of SIVΔnef was in the Tat SL8 epitope, which escaped as early as 3 weeks post-infection [26]. We hypothesized that, despite low levels of viral replication, SIVΔnef undergoes sequence evolution during the vaccination period, which would induce a shift in the specificity of the CD8 T cell response between week 5, when the animals are not protected, and week 20, where increased protection plateaus [12]. We made use of high-throughput deep sequencing to quantify SIVΔnef sequence variants and their frequencies in vaccinated animals, including those with undetectable plasma viral loads. Concurrently, we fine-mapped the CD8 T cell responses at the epitope level at week 5 and 20 post-vaccination. Studying this cohort of animals longitudinally after SIVΔnef vaccination, we demonstrate concomitant coevolution of the virus and the specificity of the CD8 T cell response, resulting in a CD8 T cell repertoire with increased anentropic specificity, wider cumulative breadth, rearranged immunodominance and enhanced depth.

Ethics statement
The animals included in this study were all female Indian-origin rhesus macaques (Macaca mulatta), housed in a biocontainment facility at the New England Primate Research Center (NEPRC). These experiments and procedures were approved by the Harvard Medical Area Standing Committee on Animals. At approval of the study, the assigned protocol number by the Institutional Animal Care and Use Committee (IACUC) was 04383. The Harvard Medical School animal management program is accredited by the Association for the Assessment and Accreditation of Laboratory Animal Care, International (AAALAC), and meets National Institutes of Health standards as set forth in the 8 th edition of the Guide for the Care and Use of Laboratory Animals [27]. The institution also accepts as mandatory the PHS Policy on Humane Care and Use of Laboratory Animals by Awardee Institutions and NIH Principles for the Utilization and Care of Vertebrate Animals Used in Testing, Research, and Training. There is on file with the Office of Laboratory Animal Welfare (OLAW) an approved Assurance of Compliance (A3431-01).
All animals were housed indoors in an SOP-driven, AAALAC-accredited facility. Husbandry and care met the guidance of the Animal Welfare Regulations, OLAW reporting and the standards set forth in The Guide for the Care and Use of Laboratory Animals. All research animals were enrolled in the NEPRC behavioral management program, including an IACUC-approved plan for Environmental Enrichment for research primates. This program included regular behavioral assessments, and provision of species appropriate manipulanda, and foraging opportunities. This protocol had an IACUC-approved exemption from social housing based on scientific justification. Primary enclosures consisted of stainless steel primate caging provided by a commercial vendor. Animal body weights and cage dimensions were regularly monitored. Overall dimensions of primary enclosures (floor area and height) met the specifications of The Guide for the Care and Use of Laboratory Animals, and the Animal Welfare Regulations (AWR's). Further, all primary enclosures were sanitized every 14 days at a minimum, in compliance with AWRs. Secondary enclosures (room level) met specifications of The Guide with respect to temperature, humidity, lighting and noise level. The animals were provided ad lib access to municipal source water, offered commercial monkey chow twice daily, and offered fresh produce a minimum of three times weekly. Light cycle was controlled at 12/12 hours daily. The animals were subject to twice daily documented observations by trained animal care and veterinary staff, and enrolled in the facility's environmental enrichment, and preventative health care programs. Euthanasia took place at defined experimental endpoints using protocols consistent with the American Veterinary Medical Association (AVMA) guidelines. Animals were first sedated with intramuscular ketamine hydrochloride at 20 mg/kg body followed by sodium pentobarbital (!100 mg/kg) intravenously to achieve euthanasia.

Peptide pools
To determine SIV protein-specific T cell responses, 9 whole protein peptide pools were constructed containing 15-mer peptides overlapping by 11 amino acids at 2.5 μg/ml for each. Peptides spanning the SIVmac239 proteome and corresponding to the peptide sequences available from the AIDS Research and Reference Reagent Program, Division of AIDS, NIAID, NIH, were synthesized by the Massachusetts General Hospital core peptide facility. The sequence of individual peptides can be found at www.aidsreagent.org.
For epitope mapping, a three-dimensional peptide matrix was constructed. The matrix was composed of 87 peptide pools, each containing around 27 to 30 peptides. The 87 peptide pools encompassed the entire SIV proteome with a coverage of 3-fold, such that each of the 823 overlapping peptides spanning the SIVmac239 proteome was represented in 3 different pools. The peptides used for epitope mapping were 15-mer peptides overlapping by 11 amino acids and were obtained from the AIDS Research and Reference Reagent Program, Division of AIDS, NIAID, NIH. The sequence of individual peptides can be found at www.aidsreagent.org.

Calculation of entropy
The entropy for each epitope within SIVmac239 was determined by calculating the entropy of each 9-amino-acid window within SIVmac239, which was determined by aligning 10 proteomes of SIVsmm, the ancestral virus of both the clonal SIVmac239 and the quasi-species SIV-mac251. The aligned proteome sequences were then uploaded into the Shannon Entropy-One calculator from the Los Alamos National Database (http://www.hiv.lanl.gov/content/sequence/ ENTROPY/entropy_one.html). The program calculates an entropy score for each 9 amino acid sequence in the SIVsmm proteome.

CD4 T cell separation
In order to detect only IFN-γ-producing CD8 T cells, ELISPOT assays were carried out with PBMCs or lymph node mononuclear cells (LNMC) after negative selection of CD4 T-lymphocyte populations fractionated by magnetic bead separation (CD4 Dynabeads; Dynal, Oslo, Norway) as previously described [20]. Negatively selected CD4 T cells were >90% CD8 T cells. CD4 T cell-depleted PBMC and LNMC were suspended in R-10 medium and used the same day in ELISPOT assays. Isolated lymph node CD4 T cells were pelleted and used for RNA extraction and cell-associated viral RNA quantitation and sequencing.

Determination of plasma viral loads
For quantitation of plasma viral loads in SIVΔnef-infected animals, highly specific, real-time RT-PCR assays were performed as described previously [28]. The assay specific for SIV-mac239Δnef was developed by designing a reverse primer that uniquely recognizes the sequence generated by the deletion of nef coding sequences in SIVmac239. The nominal threshold for this assay was 30 viral RNA copy equivalents/ml plasma.

Sequencing of SIVΔnef
To sequence lymph node cell-associated viral RNA, cell-associated viral loads were first determined, incorporating qPCR assay primers and probes for SIV gag, rhesus CCR5 and SIV-mac239Δnef as described [28]. Briefly, pelleted cells were rapidly disrupted in 100 μl lysis/ digestion solution (3M GuHCl, 50 mM TrisCl, pH 7.6, 1 mM CaCl 2 and 1 mg/ml Proteinase K), utilizing a Branson 450 sonifier equipped with a high-intensity cup horn and set at 60% power amplitude (Branson Ultrasonics, Danbury, CT). After continued digestion at 42°C for 1 hour, 400 μl of GuSCN/carrier solution (~5.7 M GuSCN, 50 mM TrisCl, pH 7.6, 1 mM EDTA, 600 μg/ml glycogen) was added to completely dissociate RNA-protein complexes. Total nucleic acids were precipitated by addition of 500 μl isopropanol, collected by centrifugation, and washed with 70% ethanol. This precipitate was air-dried, suspended in 60 μl 1X TurboDNAse buffer (ThermoFisher/Ambion), and split into 2 equal 30 μl aliquots for separate RNA and DNA preparation. The samples for DNA determinations were denatured by heating to 100°C for 5 min and quenched on ice prior to qPCR. Associated with determinations of SIVmac239Δnef DNA, cell equivalents based on qPCR for rhesus CCR5 at single haploid copy per genome (2 copy equivalents per cell) were determined on aliquots diluted 1:10 with 5 mM TrisCl, pH 9.0. Copy equivalents of SIVgag and SIVmac239Δnef targets were determined as described [28]. For the aliquots reserved for RNA, 20 μl of a cocktail containing 2 μl (4 Units) TurboDNAseI in 1X buffer (ThermoFisher/Ambion) were added to each, and the samples incubated at 42°C for 30 min to digest DNA. After digestion, RNA was recovered by addition of 200 μl of GuSCN solution without glycogen carrier, followed by precipitation with 250 μl isopropanol and a 70% ethanol wash, as noted above. Genome equivalents based on SIV gag or SIVmac239Δnef target sequences were determined as described. The threshold limits of detection for both SIV DNA and SIV RNA were 30 total copies of SIV sequence normalized to calculated cell equivalents.
Plasma samples were spun at 14,000 rpm at 4°C for 1 hour in a microcentrifuge. Viral RNA was then isolated using the Qiamp MinElute Virus spin kit according to manufacturer's instructions. Plasma viral RNA and cell-associated viral RNA (isolated and quantitated as described above) were reverse transcribed and amplified using 4 sets of primers to produce 4 overlapping amplicons spanning the entire length of the viral coding sequence using the Super-Script III High Fidelity One-Step RT-PCR kit (Invitrogen, Life Technologies, Carlsbad, CA). PCR products were either purified with the Qiagen MinElute Gel Extraction (Qiagen) kit or Agencourt AMPure XP beads (Beckman Coulter). Samples were quantified with the Quant-IT dsDNA HS Assay kit (Invitrogen) The four resultant amplicons were combined and then libraries were created and tagged using the Nextera DNA Sample Prep Kit (Illumina). Individually tagged libraries were quantified with the Quant-IT dsDNA HS Assay kit (Invitrogen) and the Agilent High Sensitivity DNA kit (Agilent Technologies). Libraries were then pooled together, denatured with NaOH, and run on an Illumina MiSeq using either a 300 or 500 cycle MiSeq Kit (Illumina, San Diego, CA, USA).

IFN-γ ELISPOT assays
IFN-γ-producing T-cell responses were enumerated using an enzyme-linked immunospot (ELISPOT) assay for detection of macaque IFN-γ (Mabtech, Mariemont, OH). CD4-depleted PBMC and LN mononuclear cells were stimulated at 2x10 5 cells per well with peptide pools (15-mer peptides overlapping by 11 amino acids at 2.5 μg/ml each). Cells were incubated overnight in multiscreen plates (Millipore) coated with an IFN-γ capture antibody, and spots representing IFN-γ-producing T cells were detected in an enzyme-linked, colorimetric assay for bound IFN-γ. Spots were counted using an automated ELISPOT plate reader (Zellnet Consulting, New York, NY).
For epitope mapping, peptide-pool-stimulated wells, containing a number of spots 3 times higher than background and greater than 30 spots per million mononuclear cells, were selected for second round deconvolution. To deconvolute peptide-specific CD8 T cell responses, peripheral and lymph node mononuclear cells, depleted of CD4 T cells, were stimulated with single peptides that were present in at least 2 positive pools of the 87 peptide pools. Wells with responses to single peptides 3 times higher than background and greater than 60 spots per million were considered significant.

Typing of MHC I alleles
Comprehensive MHC I typing using pyrosequencing was conducted as described [29]. Total cellular RNAs were converted to cDNA using a Superscript III First-Strand Synthesis System (Invitrogen). Primary cDNA-PCR amplicons spanning 190 bp of exon two of macaque class I sequences were generated with high-fidelity Phusion polymerase (New England Biolabs). Each PCR primer contained one of 12 distinct 10 bp MID tags along with adaptor sequences for 454 pyrosequencing. After purification, primary amplicons were normalized to equimolar concentrations and groups of 12 animals were pooled for GS FLX analysis. The emulsion PCR and pyrosequencing steps were performed with Genome Sequencer FLX instruments (Roche/454 Life Sciences) using GS FLX protocols according to the manufacturer's specifications (454 Life Sciences) at the 454 Sequencing Center (Branford, CT) and the University of Illinois at Urbana-Champaign High-Throughput Sequencing Center.

Statistical analysis
All statistical analyses were done using GraphPad Prism software (GraphPad Software v6.0b, Inc., La Jolla, CA, USA). Non-parametric Wilcoxon and Mann-Whitney tests were used for statistical analysis where the sample size was less than 6. Otherwise, parametric t tests were conducted; p values less than 0.05 were assumed to be significant in all analyses.

Widespread sequence evolution occurs after SIVΔnef infection
To study the effect of epitope escape on the evolution of the CD8 T cell response following vaccination with SIVΔnef, we first assessed the extent of viral sequence variation in 12 animals after SIVΔnef infection using next generation sequencing (NGS). The vaccinated animals, most of which displayed classical SIVΔnef viral replication kinetics of peak viremia at week 2 and viral set-points at or near undetectable levels (<30 copy Eq/ml) at week 8, were followed for 20 to 40 weeks before challenge (Fig. 1A). Plasma viral samples from the 12 vaccinated animals were sequenced at time points ranging from 1 week to 34 weeks post-vaccination.
RT-PCR was used to create four overlapping amplicons corresponding to the full-length SIV viral genome. Libraries were prepared from the pool of four amplicons and then sequenced using the Illumina MiSeq platform. Although there was already widespread sequence variation during the first 2 weeks after vaccination with SIVΔnef, none of the variants in any of the animals sequenced exceeded a threshold of 10% (Figs. 1C and S1), suggesting that these variants were not subject to significant selection in vivo during this initial observation period. However, by week 4 post-vaccination, as the viral load declined from peak viremia, there was a significant increase in the prevalence of nonsynonymous mutations that exceeded the 10% threshold. At this time point, plasma viral RNA samples from 12 animals displayed an average of 8 nonsynonymous mutations, each with a greater frequency than 10% (Figs. 1C and S2). To determine if sequence variation continued to increase between weeks 4 and 20, plasma viral RNA was sequenced from the 3 macaques with detectable viremia at week 20. Plasma viral RNA samples from these 3 animals at week 20 had an average of 28 nonsynonymous mutations with a frequency greater than 10%, a marked increase compared to an average of 8 nonsynonymous mutations at week 4 (Figs. 1C and S3).
To determine the extent of sequence evolution in SIVΔnef-vaccinated animals that had lower levels of plasma viremia, we analyzed lymph node CD4 T cell-associated viral RNA at week 20 from 3 animals. Follicular helper CD4 T cells, which are localized in secondary lymphoid tissues, are preferentially infected by SIVΔnef [33], and therefore have higher cell-associated viral RNA levels than CD4 T cells in peripheral blood [6]. Sequencing of lymph node CD4 T cell-associated viral RNA at week 20 revealed increased viral sequence variation even in SIVΔnef-infected animals with undetectable plasma viremia, albeit at lower levels of variation than observed for viremic animals. Lymph node associated viral RNA sequenced from animals with undetectable plasma viremia had an average of 15 nonsynonymous mutations with a greater frequency than 10% (Fig. 1D), a substantial increase over the 8 nonsynonymous mutations at week 4 for these animals (S2 Fig.).
SIVΔnef was longitudinally sequenced in 5 animals, from plasma virus at week 4 and either plasma virus or LN CD4 T cell-associated virus at week 20. Pooling all the variants (greater than 10% frequency) for all 5 animals, the average frequency of nonsynonymous mutations increased from 29% at week 4 to 47% at week 20 (p = 0.001) (Fig. 1B). Moreover, there was a significant accumulation in the number of nonsynonymous mutations, increasing from an average of 10 nonsynonymous mutations per animal at week 4 to an average of 24 at week 20 (p = 0.031) (Fig. 1E). In all animals whose virus was sequenced, there was extensive sequence evolution as measured by the number of nonsynonymous mutations and their frequencies, including in animals with undetectable plasma viral loads.
Evolution of the specificity of the CD8 T cell response in SIVΔnefvaccinated animals To determine if viral sequence evolution was accompanied by a shift in the epitopes targeted by the SIV-specific CD8 T cell response, we analyzed CD8 T cell responses at weeks 5 and 20 in the 12 vaccinated animals using IFN-γ ELISPOT assays using overlapping peptide pools spanning the entire SIV proteome. In contrast to increased protection between weeks 5 and 20, there was a 40% decrease in the magnitude of the SIV-specific CD8 T cell response in peripheral blood mononuclear cells (PBMC) between weeks 5 and 20 ( Fig. 2A). Similarly, the breadth of the CD8 T cell response did not correlate with increased protection between weeks 5 and 20, as the number of positive pools remained stable, with an average of 6 stimulating pools detected at week 5 and an average of 5.2 pools detected at week 20 in PBMCs (Fig. 2B).
A recent report demonstrated that the magnitude of the SIV-specific CD8 T cell response in lymph nodes at week 50 post-vaccination with SIVΔnef correlated with protection against intravenous challenge [6]. To determine if the kinetics of the SIV-specific CD8 T cell response in lymph nodes (LN) correlated with the increased protection observed between weeks 5 and 20, we looked at the magnitude and the breadth of the SIV-specific CD8+ T cell response in lymph nodes. Similar to the magnitude of the CD8 T cell response observed in peripheral blood, the magnitude of the SIV-specific CD8 T cell response in lymph nodes decreased modestly, by 25%, between weeks 5 and 20 (Fig. 2C). The breadth of the response, as defined by recognition of peptide pools of 30 peptides each, in the lymph nodes was likewise comparable to the breadth of the response in the periphery (Fig. 2C), with a slight decrease from 5 positive pools at week 5 to 4 positive pools at week 20.
We next examined CD8 T cell responses to Gag and Pol, the two most conserved proteins in lentiviruses, as a fraction of total responses to the SIV proteome, an indicator of the level of conservation of the SIV-specific CD8 T cell response. Intriguingly, there was a significant enrichment of the ratio of the magnitude of Gag and Pol responses to the total SIV-specific CD8 T cell response between weeks 5 and 20 in SIVΔnef-vaccinated animals (p = <0.0001) (Fig. 2D). The ratio of Gag and Pol responses to total SIV responses increased in every animal between week 5 and week 20, increasing from an average of 28% at week 5 to 47% at week 20.

Epitope mapping reveals a change in the CD8 T cell repertoire and immunodominance
To characterize the shift in CD8 T cell specificity at the level of individual epitopes, we mapped the CD8 T cell response to single peptides for 6 animals at weeks 5 and 20 post-vaccination with SIVΔnef. Epitope mapping of SIV-specific CD8 T cell responses was conducted on CD8 T cells from lymph nodes. For the 2 animals whose SIV-specific CD8 T cell responses were mapped in both lymph nodes and peripheral blood, there was no significant difference in the epitopes targeted by the CD8 T cell responses in these two sites (S1 Table).
Longitudinal epitope mapping in SIVΔnef-vaccinated animals revealed a change in the relative immunodominance of recognized epitopes between weeks 5 and 20 after SIVΔnef vaccination. Furthermore, the specific SIV epitopes recognized by CD8 T cell responses differed significantly between week 5 and week 20. Of the 36 CD8 T cell epitopes mapped in the 6 animals at week 20, only 14 were also recognized at week 5 (Figs. 3A and S4). Almost two thirds of the responses at week 20 are de novo responses or were not readily detectible at week 5.
As a first step in assessing whether these shifts in epitope specificity reflected increased anentropic specificity of the CD8 T cell response, we examined the conservation of targeted epitopes at different times after SIVΔnef vaccination. We assigned every mapped epitope a conservation score by calculating the mean entropy of nine amino acid stretches for all SIV proteins, as determined from the alignment of 10 SIVsmm predicted open reading frames (Fig. 4). SIVsmm strains were chosen because the virus clone SIVmac239 and the quasispecies SIV-mac251 are derived from SIVsmm. Interestingly, in 5 out of 6 animals, the immunodominant response targeted a more conserved epitope (lower entropy score) at week 20 than at week 5 (p = 0.039) (Fig. 3B). Moreover, 7 of the de novo CD8 T cell responses detected at week 20 targeted highly conserved epitopes with an entropy lower than 0.1.

Epitope escape drives CD8 T cell anentropic specificity
Having established the accumulation of sequence variants between weeks 5 and 20 and a simultaneous shift in the SIV epitopes targeted by the CD8 T cell response, we sought to determine whether the two events were causally related. Overlaying epitopes mapped onto the sequenced virus from the 6 animals vaccinated with SIVΔnef demonstrated that there was a high concordance between the occurrence of nonsynonymous mutations and identification of CD8 T cell epitopes. In macaque 225-97, for which we have the most comprehensive set of viral sequence data and longitudinal epitope mapping, almost all of the nonsynonymous mutations (7 out of 9 mutations) with a frequency greater than 50% fell within a mapped CD8 T cell epitope (Fig. 5A), implying a causal link between CD8 T cell pressure and viral escape. Moreover, the nonsynonymous mutations and mapped epitopes not only overlapped spatially but also overlapped temporally. In animals 225-97 and 256-00, epitope-specific CD8 T cell responses declined in frequency following escape of the epitope (Fig. 5B). Remarkably, escape epitopes in animal 225-97 were themselves overtaken by other epitope variants over time, suggesting CD8 T cell pressure on the escape epitopes as well (Fig. 5C).
However, we did not observe sequence variation in all identified CD8 T cell epitopes. As expected, we observed a strong bias favoring the escape of responses targeting highly variable epitopes, as determined by calculated entropy scores. An analysis of the top four responses mapped for each animal at weeks 5, 20 and 40 after SIVΔnef vaccination demonstrated that the escaped epitopes had significantly higher entropy scores than targeted epitopes that did not escape (p = 0.0013) (Fig. 6A). Similarly, SIV-specific CD8 T cells can be grouped into maintained responses, which target mostly conserved regions, and waning responses that target escaped epitopes. The magnitude of the variable responses, defined as targeting epitopes with an entropy above 0.25, decreased significantly between weeks 5 and 20 (p = .0101) and between weeks 5 and 40 (p = .0449), whereas the magnitude of the conserved responses was maintained (Fig. 6B).
To determine if the selective escape of variable epitopes, coupled with the decline of their CD8 T cell response frequencies, led to an immune response that was increasingly focused on more conserved epitopes, we developed a formula that estimates the conservation level of the immune response. The Response Conservation Index is the sum of all SIV-specific CD8 T cell response frequencies, weighted for entropy. The Response Conservation Index was calculated for every animal at both week 5 and 20 using the formula in Fig. 7. Importantly, in the six animals whose CD8 T cell responses were mapped, the Response Conservation Index increased significantly between weeks 5 and 20 (p = 0.0156) (Fig. 7).
To ascertain that the increased targeting of conserved epitopes in the vaccinated animals was not due to a genetic background enriched for protective MHC I molecules, we conducted comprehensive pyrosequencing of MHC I genotypes in the six animals whose CD8 T cell responses were longitudinally mapped. The cohort displayed a diverse set of MHC I molecules with three animals expressing the protective alleles Mamu-AÃ01 and Mamu-BÃ17, while the rest of the animals did not express any protective MHC I alleles (Table 1). Notably, increased anentropic specificity was documented in every animal between weeks 5 and 20 regardless of the MHC I background (Fig. 7).
Comprehensive MHC I genotyping identified major alleles expressed in the six animals whose CD8 T cell responses were mapped longitudinally. Major MHC I alleles were defined as

Epitope escape drives increased breadth and depth of the CD8 T cell response
The number of epitopes targeted by the CD8 T cell response did not differ between weeks 5 and 20, as measured by positive peptide pools (Fig. 2B, 2C) or the specific mapped epitopes (Fig. 8A). However, we hypothesized that escaped CD8 T cell responses were not deleted but rather maintained as a central memory T cell pool and, therefore, that the cumulative breadth of the CD8 T cell repertoire, encompassing all previously detected responses in a given animal, may be a more accurate measure of SIV-specific CD8 T cell breadth (Fig. 8B). Phenotypic analysis of SL8-specific CD8 T cells, targeting the highly mutable SL8 epitope, did reveal a shift to a central memory profile between weeks 5 and 20, characterized by increased cell surface expression of CD28 (p = 0.0087), CCR7 (p = 0.0006) and CD127 (p<0.0001), as well as a decrease in intracellular expression of perforin (p<0.0001) and Ki67 (p = 0.0194) (Fig. 8C). In Mamu-A1Ã001+ (previously known as Mamu-AÃ01) animals, SL8-specific CD8 T cell responses are frequently immunodominant during acute infection but quickly decline following the evolution of escape variants, where 82% and 94% of sequenced SL8 epitopes had mutated as early as week 4 in two SIVΔnef-vaccinated macaques (Fig. 8D). Although SL8-specific CD8 T cells reached nearly undetectable levels by week 20 as measured by IFN-γ ELISPOT and SL8-MHC tetramer staining, the responses rebounded in both animals following challenge with SIVmac251 (Fig. 8D). Given the maintenance of responses to escaped epitopes in SIVΔnef-vaccinated animals, all escaped responses were included in the calculation of cumulative breadth. In contrast to CD8 T cell breadth, the cumulative breadth of the CD8 T cell response increased significantly from an average of 6 responses per animal at week 5 to an average of 10 responses at week 20 (p = 0.0036) (Fig. 8B). Finally, SIVΔnef escape, in addition to allowing for anentropic specificity and increased cumulative breadth, contributes to the deepening of the CD8 T cell response. In macaque 256-00, the mapped Mamu-B17-restricted FW9 (FHEAVQAVW) epitope was intact at week 5 but the virus fully escaped by week 20 with 99.9% of sequenced reads containing the FW9-H831Y variant. To determine if de novo CD8 T cell responses to the escape variant arose, we tested both variants at week 5 and at the latest time point before challenge, week 38. Lymph node CD8 T cells from macaque 256-00 recognized FW9 at week 5 but not at week 38; conversely, the FW9-H831Y escape variant was recognized only at week 38 but not at week 5 (Fig. 9A). Similarly, the GY9-K76R epitope variant was recognized at week 38 but not at week 5 demonstrates increased depth of the CD8 T cell response after the escape of the mapped Mamu-AÃ02-restricted GY9 (GSENLKSLY) epitope (Fig. 9B). These data demonstrate that, at least in a subset of SIVΔnef-vaccinated animals, the depth of the CD8 T cell response increased between weeks 5 and 20 as new CD8 T cell responses were mounted against escape variants.

Discussion
Taken together, this data set supports a comprehensive model that provides new insights into the mechanism of SIV-specific immunity maturation after SIVΔnef vaccination. For the first time, we demonstrate that SIVΔnef sequence evolution is widespread after vaccination, not just confined to the previously documented highly mutable epitope of Tat SL8 [26], and that ongoing sequence evolution occurs even in animals with undetectable plasma viremia. Moreover, we show that the observed sequence evolution can be characterized as viral epitope escape, a consequence of SIV-specific CD8 T cell pressure on the virus as indicated by the many nonsynonymous mutations that appear within mapped CD8 T cell epitopes. Importantly, another direct consequence of immune CTL pressure is the disproportionate escape of highly variable epitopes, which leads to the waning of CD8 T cell responses targeting entropic epitopes and the accumulation of CD8 T cell responses specific for these epitopes. In addition, escaped SIVspecific CD8 T cell responses are maintained as central memory CD8 T cells and, as such, contribute to the increased cumulative breadth of the CD8 T cell response at week 20 post-vaccination with SIVΔnef. Finally, epitope escape and the ensuing generation of de novo CD8 T cell responses to these escape variants constitute a deepening of the CD8 T cell repertoire. Our work sheds new light on SIVΔnef and other LASIV vaccines and how their functional and mechanistic characteristics are fundamentally different from other vaccines. While many of the recent successful replication competent virus vector-based vaccines such as CMV and RRV are also persistent [34][35][36][37], we demonstrate that SIVΔnef is not only a persistent vaccine but that, as a lentivirus, it evolves with the immune response it induces and, ultimately, shapes the CD8 T cell response during the vaccination period. This unique characteristic of SIVΔnef to mutate allows for maturation of CTL specificity through increased targeting of conserved regions, deepening and broadening of the response. Furthermore, this new model accommodates earlier known properties of LASIV, such as the observation that the level of protection of LASIV vaccines inversely correlates with the level of attenuation of the vaccine [6,38]. Based on our current work, we hypothesize that LASIV vaccine efficacy depends on virus evolution, resulting in the concurrent CD8 T cell specificity shift from highly variable to conserved regions of the virus. Highly attenuated viruses, such as SIVΔ4, which exhibits lower efficacy compared to SIVΔnef [38], would have lower replication rates and less escape events, therefore less opportunity to shape the CD8 T cell response.
Despite three decades of HIV research, basic questions of vaccine design remain hotly debated, namely the optimal balance of CD8 T cell response magnitude, breadth and specificity that should be induced by a vaccine. In suggesting that anentropic specificity is a more important correlate of protection in SIVΔnef-vaccinated animals than response magnitude, our data corroborates previous work on long-term nonprogressors (LTNP), another important model of spontaneous lentiviral suppression. Research on LTNP has shown that MHC I molecules associated with protection restrict CD8 T cells that target regions within HIV Gag [39,40] as well as highly conserved epitopes [41]. However, unlike viral control in LTNP, which is associated with certain protective HLA molecules, we document increased targeting of highly conserved epitopes in all animals regardless of their MHC I genotypic background. Similar to our LASIVinduced protection, viral control in LTNP is correlated with the immunodominance of CD8 T cell responses that target highly conserved epitopes [42], given that immunodominant CD8 T cells exert more immune pressure [41]. Finally, studies of LTNP revealed the importance of targeting conserved regions during acute infection [41], before massive depletion of CD4 T cells has occurred [43]. Our research demonstrates that at week 20, whereas responses to variable regions have waned, the immunodominant SIV-specific CD8 T cells target more conserved epitopes than at week 5. The presence of immunodominant CD8 T cells exerting pressure on highly conserved regions during acute infection in animals challenged at week 20 post-vaccination, but not at week 5, could explain the observed difference in protection. Conversely, the absence of immunodominant responses targeting highly conserved regions in animals challenged at week 5 allows the virus to easily escape and the accumulation of anentropicspecific CD8 T cell responses would occur long after acute infection, if at all.
In addition to immunodominant SIV-specific CD8 T cells with anentropic specificity, we also show increased cumulative breadth and depth of the CD8 T cell response at week 20 compared to week 5. The increased cumulative breadth at week 20 is a function of waning escaped CD8 T cell responses that are maintained as central memory CD8 T cells and are recalled upon challenge. It is unclear how effective these responses are at long-term viral suppression, given that they have been escaped and mostly target variable regions of the virus. Still, a recent report suggests that acute-phase SIV-specific CD8 T cells targeting variable regions are important for viral control [44], indicating that the increased cumulative breadth of the CD8 T cell response at week 20 may also contribute to protection. Moreover, we detected de novo CD8 T cell responses specific to escaped epitope variants, illustrating that viral escape also contributes to increased immune response depth. Although LTNP generate HIV-specific CD8 T cells that cross-react with escape mutants [45], the generation of de novo CD8 T cell responses to escaped epitopes is rarely observed in adults infected with HIV [46]. In animal 256-00, we identified two such CD8 T cell responses that target variant epitopes at week 38 but were not detectable at week 5. We only assayed CD8 T cell responses against escape variants of known CD8 T cell epitopes with known MHC I restriction, and only if the MHC I is expressed by the infected animal. The relative frequency of de novo CD8 T cell responses in SIVΔnef-vaccinated animals, if verified in future studies, could be due to an intact CD4 T cell compartment that persists well into the chronic phase of LASIV infection.
Previous reports have documented rapid escape of variable epitopes during acute HIV-1 infection and maintenance of CD8 T cell responses specific for conserved epitopes into chronic infection [47][48][49]. Here, in the context of SIVΔnef, we demonstrate, in addition to rapid escape of variable epitopes and maintenance of existing CD8 T cell responses with anentropic specificity, de novo CD8 T cell responses against novel epitopes and escape variants during chronic infection. Furthermore, we correlate the accumulation of CD8 T cell responses with anentropic specificity and the increased repertoire depth and breadth to the kinetics of SIVΔnef vaccineinduced protection. Previous research on the coevolution of lentiviral epitopes and the CD8 T cell response has been conducted in the setting of wild-type infection, and it is unclear that our findings of increased anentropic specificity, cumulative breadth and repertoire depth during SIVΔnef vaccination would apply to wild-type lentiviral infection, in which the CD8 T cell response is likely to be impaired by the persistently high antigenic load, CD4 T cell depletion in the gut and microbial translocation-induced immune activation.
Notably, we did not see an increased magnitude of the SIV-specific CD8 T cell response in the lymph nodes between weeks 5 and 20. Fukazawa and colleagues [6] have recently shown that less attenuated LASIVs such as SIVΔnef have higher viral replication in the lymph nodes and, in turn, have higher lymph node CD8 T cell responses. Although these results corroborate existing data demonstrating that protection inversely correlates with LASIV attenuation [38], our data demonstrates that the magnitude of the CD8 T cell response does not correlate with the kinetics of protection, neither in the PBMC nor lymph node compartment. Although our current work supports the importance of conserved responses as a correlate for protection, others have shown that using conserved regions as immunogens may result in a reduction of the total magnitude of responses [50]. It is unknown whether vaccines inducing targeted responses toward conserved regions will make up for the loss of magnitude. Still it is encouraging that we observe increased anentropic specificity in all animals, suggesting that even nonprotective MHC I molecules can target highly conserved epitopes if properly primed.
It has also been proposed that the increased protection of LASIV vaccines is due to phenotypic maturation of the CD8 T cell response. It is important to delineate whether observed phenotypic changes in SIV-specific CD8 T cells correlate with LASIV protection and whether CD8 T cell phenotypic maturation is specific to LASIV or can be seen in other vaccines, such as adenovirus-, CMV-and RRV-based SIV vaccines. Finally, although this data set was focused on the CD8 T cell response and its evolution over time post-vaccination, we observed relatively few nonsynonymous mutations in env outside of mapped CD8 T cell epitopes, suggesting that antibody responses induced only modest immune selective pressure on SIVΔnef compared to CD8 T cell responses.
Although immunodominance, specificity, breadth and depth have been implicated for some time as determinants of a successful CD8 T cell response against HIV/SIV infection, these results demonstrate how such a CD8 T cell response may be induced by the highly protective SIVΔnef vaccine. This new model of SIVΔnef-mediated protection raises obvious implications for vaccine design, namely that a successful vaccine should focus CD8 T cell responses on conserved regions of the virus and their existing variants. The recent success of viral vectors such as CMV [34][35][36], RRV [37], and adenovirus [51] to elicit potent CD8 T cell responses that partially protect rhesus macaques from infection has been an important step forward in vaccine design. However, these vaccines are all derived from DNA viruses, which have significantly lower rates of mutation than lentiviruses, and thus they are unlikely to shape the immune response in a similar fashion to SIVΔnef. Our work suggests that immunogens should be carefully chosen for these vectors such that the induced responses preferentially target conserved regions of the virus, as well as existing variants of these regions to increase response breadth and depth simultaneously.
Supporting Information S1 Table. Epitope mapping of the CD8 T cell response in peripheral blood and lymph nodes. The mapped CD8 T cell response in peripheral blood (PBMC, left) and lymph nodes