Identification of immunodominant T cell epitopes induced by natural Zika virus infection

Abstract

Zika virus (ZIKV) is a flavivirus primarily transmitted by Aedes species mosquitoes, first discovered in Africa in 1947, that disseminated through Southeast Asia and the Pacific Islands in the 2000s. The first ZIKV infections in the Americas were identified in 2014, and infections exploded through populations in Brazil and other countries in 2015/16. ZIKV infection during pregnancy can cause severe brain and eye defects in offspring, and infection in adults has been associated with higher risks of Guillain-Barré syndrome. We initiated a study to describe the natural history of Zika (the disease) and the immune response to infection, for which some results have been reported. In this paper, we identify ZIKV-specific CD4+ and CD8+ T cell epitopes that induce responses during infection. Two screening approaches were utilized: an untargeted approach with overlapping peptide arrays spanning the entire viral genome, and a targeted approach utilizing peptides predicted to bind human MHC molecules. Immunoinformatic tools were used to identify conserved MHC class I supertype binders and promiscuous class II binding peptide clusters predicted to bind 9 common class II alleles. T cell responses were evaluated in overnight IFN-γ ELISPOT assays. We found that MHC supertype binding predictions outperformed the bulk overlapping peptide approach. Diverse CD4+ T cell responses were observed in most ZIKV-infected participants, while responses to CD8+ T cell epitopes were more limited. Most individuals developed a robust T cell response against epitopes restricted to a single MHC class I supertype and only a single or few CD8+ T cell epitopes overall, suggesting a strong immunodominance phenomenon. Noteworthy is that many epitopes were commonly immunodominant across persons expressing the same class I supertype. Nearly all of the identified epitopes are unique to ZIKV and are not present in Dengue viruses. Collectively, we identified 31 immunogenic peptides restricted by the 6 major class I supertypes and 27 promiscuous class II epitopes. These sequences are highly relevant for design of T cell-targeted ZIKV vaccines and monitoring T cell responses to Zika virus infection and vaccination.

1 Introduction

Flaviviruses including Dengue, West Nile, Yellow Fever and Zika are positive sense, single-stranded RNA viruses transmitted by mosquitoes in regions primarily located along the equator. Zika virus (ZIKV) was first identified in a monkey from the Ugandan Zika Forest in 1947, and the first confirmed human case was reported in Uganda in 1963 (1). Few confirmed human cases were found for the next several decades, until an outbreak in 2007 in the Yap Islands resulted in 49 confirmed cases; impressively, it is estimated that 75% of the population became infected during this outbreak (1). Individuals infected with ZIKV were asymptomatic or displayed only mild illness characterized by fever, rash, joint pain and/or conjunctivitis. A larger outbreak occurred in French Polynesia in 2013-14 in which 340 confirmed cases were found, but as many as 49% of the population of ~270,000 was estimated to have been infected based upon serological surveys (2). During this outbreak, a higher incidence of Guillain-Barré Syndrome (GBS) was observed (3). In 2014, autochthonous transmission of ZIKV was detected in Chile (Easter Island) coinciding with transmission in the Pacific Islands, and by 2015, circulation in Brazil was widespread and ZIKV infection was associated with cases of GBS as well as birth defects. Most concerning were associations with severe brain (microcephaly) and eye defects and neurodevelopmental abnormalities in infants born to women with ZIKV infections during pregnancy (1, 4). In one study, 1450 children born to women with laboratory evidence of ZIKV infection during pregnancy were evaluated: 203 (14%) had a ZIKV-associated neurologic abnormality or birth defect, and 84 (6%) had microcephaly (5). Though much has been learned in the past several years, the exact causes of ZIKV infection-induced neurologic complications during fetal development and adulthood have yet to be elucidated.

There are currently no FDA-approved treatments for ZIKV infection or preventive vaccines. ZIKV infection induces rapid and durable neutralizing antibody and T cell responses in humans. In an observational longitudinal study of ZIKV-infected persons, we reported that ZIKV-specific neutralizing antibody and T cell responses were elicited early during infection in all ZIKV-infected participants and these titers remained elevated throughout the duration of the study year (6). ZIKV neutralizing antibody titers were significantly higher in dengue virus (DENV)-experienced participants compared with DENV-naïve subjects, consistent with the known cross-reactive nature of antibodies induced by other flavivirus infections (7). ZIKV-specific antibodies, including those present in human sera after receiving formalin-inactivated Zika vaccine (ZPIV), can provide protection against ZIKV challenge in animal models (8). Antibody-dependent enhancement (ADE) has been reported during heterologous infections with different flavivirus infections (e.g., DENV) including in animal models of ZIKV infection (9, 10). However, dengue-naïve and dengue-exposed individuals exhibited similar viral loads and cytokine profiles during ZIKV infection (11). Thus, ADE remains an unlikely cause of severe disease during human ZIKV infection based on epidemiologic data. It is worth noting that most Zika vaccine development strategies focus on generation of neutralizing antibody responses.

We now report a detailed study of T cell responses induced by ZIKV infection in the cohorts described above. Flavivirus-specific T cells are induced rapidly in infected individuals and have been shown to mediate protection in animal models (12–14). In fact, CD8+ cells induced by DENV infection protect against ZIKV infection (15, 16) similar to findings from other heterologous flavivirus studies (17, 18). HLA-B7 transgenic mice vaccinated with DENV/ZIKV cross-reactive HLA-B7-restricted epitopes developed CD8+ T cell responses that were protective against ZIKV infection (19). Several ZIKV T cell epitopes have been identified using known DENV epitopes as the starting point for sequence-based comparisons. In the study presented here, we utilized advanced matrix-based MHC-binding prediction algorithms for 6 human class I and 9 human class II supertypes to identify T cell epitopes spanning the entire ZIKV open reading frame.

2 Materials and methods

2.1 Study design and participant description

Volunteers with suspected ZIKV infection after recent travel to Zika-endemic regions were recruited at US sites between July 2016 and September, 2017 as part of a natural history study. The study design and primary results were published previously (6). ZIKV infection was established by ZIKV PCR positive testing or ZIKV-specific IgM and neutralizing antibody test results indicative of infection (6). It was reported that 68.9% of the recruited Zika-infected donors were female, 28.9% Hispanic (71.1% non-Hispanic) and the median age was 44 (6). In this study, a subset of samples was evaluated from 14 ZIKV-naïve and 41 ZIKV-infected participants collected 1-6 months post-predicted infection. Of the 41 ZIKV-infected participants, 13 were considered flavivirus-experienced (6) based upon high baseline DENV-specific neutralizing antibody titers (≥250). Peripheral blood mononuclear cells (PBMC) were isolated from whole blood using standard techniques, cryopreserved, and stored in vapor phase nitrogen until use. Saliva samples were utilized for HLA typing (ProImmune), and 4-digit resolution typing results were mapped to class I and II supertypes based on work described by Southwood and Sidney et al (20, 21).

2.2 Immunoinformatics

Algorithms developed at EpiVax, Inc were used to identify CD4+ and CD8+ T cell epitopes from the reference ZIKV strain PRVABC59 (3,423 amino acids; Genbank protein accession number AMC13911.1). The immunoinformatic algorithms used in this work are available in the Interactive Screening and Protein Reengineering Interface (ISPRI) web-based toolkit, which provides seamless transition of data from silo to silo (22–30). Predicted binding for 9 common class II alleles (DRB1*0101, DRB1*0301, DRB1*0401, DRB1*0701, DRB1*0801, DRB1*0901, DRB1*1101, DRB1*1301, and DRB1*1501) expressed by >95% of the human population was assessed using the matrix-based MHC binding algorithm EpiMatrix (22–24). This algorithm scores individual 9-mer frames for predicted binding to each allele, and assessments with Z-scores above 1.64 were considered potential binders. Though Z-scores are a measure of binding probability rather than binding affinity, in practice higher Z-scores are correlated with both higher rates of binding and higher binding affinity. The binding groove of HLA class I and HLA class II molecules can accommodate, in most cases, 9 amino acids. The binding groove of the class I molecule is “closed-ended”. As a result, these molecules can only bind short peptides, usually ranging from 8 to 11 amino acids in length. The binding groove of the class II molecule is “open-ended” and can accommodate longer peptides. However, the core of the binding groove, the part that controls the interaction between the peptide ligand and the presenting MHC is still limited to just 9 amino acids. Class II ligands are longer, typically 15 to 25 amino acids in length. This is because class II ligands must be stabilized by the interaction of leading and trailing “flanking” residues with the outside of the class II HLA molecule. The modeling tool used at EpiVax, the EpiMatrix system, evaluates all possible 9-mer cores for binding. The Clustimer algorithm finds cores that are overlapped and constructs longer peptides that contain multiple binding cores and adds n- and c-terminal flanking residues in order to construct peptides that are appropriate for synthesis and testing. EpiMatrix was also employed to identify peptides predicted to bind alleles representative of the 6 major class I supertypes (A*0101, A*0201, A*0301, A*2402, B*0702, and B*4403) cumulatively expressed by >95% of humans. The ClustiMer algorithm (25, 26) was used to identify longer sequences (16-27 amino acids) predicted to bind multiple class II alleles, and cluster scores above 10 were further evaluated.

Complete sequences from 104 additional ZIKV strains available in 2016 were utilized to evaluate the conservation of epitopes identified from strain PRVABC59 in diverse clinical isolates (from Africa, Asia and the Americas) using the Conservatrix algorithm (27). JanusMatrix (28–30) is a unique EpiVax tool which was used to investigate potential cross-reactivity between predicted ZIKV epitopes and human proteins as well as between different flaviviruses. This tool is particularly useful in identifying (and triaging) “cross-conserved” epitopes which have human homology at the T cell receptor (TCR) face and are more likely to cause off-target effects (e.g., autoimmunity or regulatory T cell responses). Only peptides with low cross-conservation with human proteins (JanusMatrix homology score below 2) were selected. Selected peptides were synthesized to high purity (>80%) by 21^st Century Biochemicals, Inc, suspended to 1 mg/ml (each individual peptide) in dimethyl sulfoxide (DMSO) and frozen at -80°C in single-use aliquots until use. Cross-conservation analyses of confirmed epitopes were performed using a panel of DENV strains by JanusMatrix. Additionally, validated peptide sequences shown to be immunogenic in samples from ZIKV-infected participants expressing the appropriate HLA alleles were evaluated using the Immune Epitope Database (IEDB) (31) to determine the novelty of the identified epitopes.

2.3 Zika virus and peptide library

ZIKV strain PRVABC59 (BEI 50240) was propagated in Vero cells (ATCC CCL-81) and quantified using focus forming assays (FFA) as previously described (13, 32). Overlapping peptides (15-mers, overlapping by 10 amino acids) spanning the entire polyprotein of ZIKV strain PRVABC59 (GenBank: AMC13911.1) were synthesized and prepared into pools of 49 peptides each (except the final pool, N, which consisted of 46 peptides) by 21^st Century Biochemicals, Inc, as shown in Table 1. Peptides were suspended to 1 mg/ml (each individual peptide) in DMSO and stored at -80°C until use.

2.4 ELISPOT assays

PBMC were thawed, washed, and suspended in ELISPOT assay medium (RPMI1640 medium supplemented with 10% fetal bovine serum, L-glutamine, penicillin and streptomycin). IFN-γ ELISPOT assays were set up according to manufacturer recommendations (BD Biosciences #551849). Briefly, antibody-coated plates were washed and blocked with ELISPOT medium for at least 2 hours at room temperature. PBMC (3x10⁵/well) were stimulated with DMSO as a negative control (0.1% final concentration matched to peptide condition concentrations), peptide pools (1 µg/ml of each peptide), individual peptides (1 µg/ml), ZIKV (multiplicity of infection = 1), and Phorbol 12-myristate 13-acetate as a positive control (PMA; 5 ng/ml). After 22-24 hours of incubation, ELISPOT assays were developed according to manufacturer recommendations. Spots were visualized using aminoethyl carbazole (AEC) substrate solution, reactions stopped using multiple water rinses, and plates decontaminated by submersion in 2% formaldehyde. ELISPOT plates were scanned using a CTL S6 Universal-V Analyzer (ImmunoSpot 6.0 plate reader; ImmunoCapture Version 6.6) and spots enumerated with ImmunoSpot CTL cell counting software (version 1.3).

2.5 Statistics and analysis

Proprietary immunoinformatic tools developed at EpiVax were utilized for selection of putative T cell epitopes and epitope clusters. The numbers of spot-forming cells (SFC) per million PBMC are shown in graphs generated using GraphPad Prism, version 9. Statistical analyses of composite data from different cohorts of participants were evaluated using Wilcoxon rank-sum tests. The threshold for positive responses to a specific peptide were determined by being a) ≥ 10 SFC/million PBMC, b) >2 standard deviations (SD) above the mean of the sample’s DMSO response, and c) > 2 SD above the mean of all samples’ DMSO responses. Hit rates for specific peptides were calculated as the number of participants with a positive response divided by the number of participants tested.

3 Results

3.1 Immunoinformatic selection of human T cell epitopes expressed by ZIKV

Two strategies employing IFN-γ ELISPOT assays were developed to evaluate T cell responses and their specificities (molecular targets) induced by ZIKV infection: targeted immunoinformatics vs. unbiased whole genome approaches. The overall strategy to identify putative T cell epitopes is shown in Figure 1A. The full-length polyprotein sequence from ZIKV strain PRVABC59 was parsed into 9-mer frames and scored for predicted binding to 9 common class II DRB1 alleles cumulatively expressed by 95% of the population using the previously described and validated EpiMatrix algorithm (33). For each allele, frames scoring >1.64 on the EpiMatrix Z-scale are determined as likely MHC binders. In parallel, the Conservatrix algorithm was used to evaluate the conservation of each scored 9-mer frame within a panel of 104 ZIKV isolates. The ClustiMer tool was used to identify overlapping 9-mer frames that contain predicted class II epitopes, and generated T cell clusters which are longer, epitope dense sequences with broad promiscuity (25, 26). These T cell clusters were examined for conservation with human proteins utilizing the JanusMatrix algorithm which can identify sequences predicted to bind HLA and have similar residues interacting with the TCR face. Clusters with high homology (JanusMatrix >2.0) to human proteins overall and at the TCR face were identified and excluded from further study, to avoid identification of epitopes with potential to induce regulatory/suppressive T cells or autoimmunity (28–30). The top clusters selected for synthesis and in vitro testing ranged from 16-27 amino acids in length, contained multiple predicted MHC ligands (combined Z-scores of 10.9 – 48.5; average of 25.7), and were highly conserved in the panel of 104 ZIKV isolates (95% on average at ≥80% sequence identity). An example of an EpiMatrix Cluster Report for Cluster DR1 is shown (Figure 1B) with the locations of the selected promiscuous DR clusters (DR) in the ZIKV genome (Figure 1C). A similar approach was utilized to identify putative CD8+ T cell targets. The ZIKV polyprotein sequence was parsed into 9- and 10-mer frames, and predicted binding to each of 6 HLA alleles representative of major class I supertypes (A1, A2, A3, A24, B7 and B44) was performed using EpiMatrix as described above. Selected epitopes on average were conserved at ≥80% amino acid homology in >95% of the 104 strains and were present throughout the length of the entire ZIKV polyprotein (Figure 1C).

Figure 1

3.2 Evaluation of T cell responses in volunteers previously infected with ZIKV

PBMCs from ZIKV^pos and ZIKV^neg participants were stimulated in overnight IFN-γ ELISPOT assays with live ZIKV, pools of overlapping ZIKV peptides (A-N), and pools of peptides predicted to bind promiscuously to HLA-DR alleles, or to each of the 6 MHC class I supertypes. PBMCs were only stimulated with peptide pools matching their HLA class I typing results (e.g.,HLA A1+, A2+, B7+, B44+ PBMCs stimulated with A1, A2, B7, and B44 peptide pools). Shown in Figure 2 are results from four different ZIKV^pos volunteers. All four of the volunteers responded to the promiscuous panDR clusters (DR) and at least one pool predicted to bind to one of the volunteer’s MHC I supertypes. Responses to live ZIKV and the overlapping 15-mer peptide pools were varied, and T cell responses to the overlapping peptides were lower, on average, than the maximal responses induced by the pools predicted using immunoinformatic tools.

Figure 2

These observations are further demonstrated across the cohort of 41 ZIKV^pos participants (Figures 3A, B). A total of 14 participants determined to be ZIKV^neg (uninfected) were also studied in these assays, and very little background was observed in these individuals. Statistically significant differences were observed between ZIKV positive and negative participants for all the tested overlapping and predicted peptide pools, as well as for live ZIKV (Figures 3A, B). We next determined whether live ZIKV and ZIKV peptide-specific T cell responses were different in Zika-infected, but previously flavivirus-naïve participants, versus responses from previously flavivirus-experienced, Zika-infected participants (determined by dengue serology). Shown in Figures 3C, D are results from ZIKV^pos participants segregated by previous flavivirus experience. Responses to overlapping and predicted ZIKV peptide pools were numerically higher in ZIKV^pos but previously flavivirus naïve subjects than in ZIKV^pos and flavivirus-experienced participants, with some response differences achieving statistical significance (overlapping pools D and J, and HLA-A2- and panDR-predicted pools).

Figure 3

3.3 Identification of individual epitopes eliciting responses in PBMCs from ZIKV-infected participants

After determining the peptide pools which induced positive responses in each volunteer, we next sought to identify the individual peptides responsible for the positive responses. For example, shown in Figure 4 are peptide pool screening assay results from participants that showed responses to the B44 pool (A) or the B7 pool (C). In the follow-up ELISPOT assays, PBMCs from these same volunteers were stimulated with the same peptide pool (e.g., B44 and B7 pools) as well as with each of the individual peptides. As shown in panel B, the majority of T cells from the volunteer responding to the B44 peptide pool were directed to a single peptide: B44-08 (NS5₂₈₁₆ DENHPYRTW). T cells from the other volunteer responsive to the B7 peptide pool responded to only two individual peptides contained within the B7 pool (panel D; B7-06: NS5₂₉₁₁ RPRVCTKEEF and B07-07: NS2A₁₂₉₀VPRTDNITL).

Figure 4

Similar assays were conducted with additional volunteers, and positive responses were determined based on calculations described in Methods. Shown in Figure 5 are the positive response rates of volunteers tested for each class I supertype (panels A-F) and promiscuous class II cluster sets (panel G). Only 2 epitopes were identified using PBMCs from 10 participants expressing HLA-A1 (panel A). In contrast, 8 individual HLA-A24-restricted T cell epitopes were identified (panel F), using PBMCs from just 4 participants expressing HLA-A24 supertype alleles. Diverse responses to the class II panDR clusters were observed, with 27 unique peptides inducing responses in at least 1 of the 34 participant PBMCs tested. Many class I supertype responses in different volunteers focused on a single or select few of the predicted epitopes. In addition, we found many examples where at least 50% of volunteers expressing alleles within a given supertype responded to the same immunodominant epitopes. In contrast, response rates to the different DR clusters were in general more diverse, with lower frequencies of responses to individual epitopes than those observed for class I. In total, 31 MHC I-restricted T cell epitopes and 27 MHC II-restricted epitope clusters were identified, as summarized in Tables 2–4. The identified epitopes are located throughout the ZIKV proteome in both structural and non-structural gene products and conserved throughout various ZIKV isolates. It is worth noting that many of the identified epitopes have been reported previously, and nearly all of the T cell epitopes identified using our strategy are ZIKV-specific as they are not conserved in other flaviviruses, including DENV 1-4 isolates (Supplemental Table 1).

Figure 5

4 Discussion

Human B cell and antibody responses to ZIKV infection have been well characterized by multiple groups, including ours (6, 34–36), and over the past few years much has been learned about ZIKV-specific T cells (6, 37–41). In our current work, we corroborate findings from other studies demonstrating strong ZIKV-specific effector T cell responses after ZIKV infection (Figure 3). In addition to identifying potent T cell responses to live ZIKV stimulation, we also found heterogenous T cell responses to overlapping peptide pools and pools of peptides predicted to bind class I supertypes and common class II HLA molecules (Figures 2, 3). These responses are consistent with those reported elsewhere in which ZIKV-infected participant T cells responded to ZIKV peptide megapools containing multiple overlapping (or predicted) ZIKV peptides (42, 43). Using these immunoinformatically predicted peptides and a two-step ELISPOT method, we identified and validated 58 T cell epitopes (31 MHC I-restricted CD8+ T cell epitopes, and 27 promiscuous MHC II-restricted CD4+ T cell epitope clusters, summarized in Tables 2–4) which warrant further study as potential components of novel T cell-targeted ZIKV vaccines.

The T cell studies presented here extend the data previously published using samples from the same cohort of Zika-infected participants (6). Both our initially published data and the work reported here showed that ZIKV-specific CD4+ and CD8+ T cell responses were targeted to sequences across the ZIKV polyprotein. Also, in both reports, T cell responses in acutely ZIKV-infected persons were higher in previously flavivirus-naïve participants than in previously flavivirus-experienced subjects. It is possible that acute ZIKV infection results in higher T cell responses in flavivirus-naïve individuals because these persons develop higher viral loads in the absence of pre-existing cross-reactive neutralizing antibody responses. This hypothesis is consistent with previous reports indicating that despite significant cross-reactivity between ZIKV and Dengue comparing B cell epitopes, very little cross-reactivity between ZIKV and other flavivirus has been found studying T cell epitopes (42). Our results presented here do provide distinct and novel information compared with our earlier publication in that: 1) the assays used were different (IFN-γ ELISPOT versus direct ex vivo intracellular cytokine staining), 2) the cytokines used to detect positive responses were different (IFN-γ alone associated with protective effector functions versus at least one of three cytokines measured by ICS including IFN-γ, IL-2 or TNF, which may detect less differentiated T cells which may not contribute to optimal protective effects), and 3) we focus here on more broadly relevant T cell epitopes across HLA supertypes that might provide better population level coverage if included in future T cell-targeted vaccines.

In general, it appears that peptide pools generated using advanced algorithms out-performed overlapping peptide approaches (Figure 3). Median magnitudes of T cell responses to pools of predicted epitopes were up to four times higher than responses to overlapping peptide pools. This suggests that the immunoinformatic prediction-based strategy may be more useful for identification of relevant epitopes. While all 9-10-mer class I epitopes would be represented in at least one peptide in the 15-mer overlapping peptide pools, it is possible that some class II peptides would not be ideal for binding to class II molecules (due to truncated or off-center binding motifs and/or missing flanking positions). The immunoinformatic predictions included epitopes shown to induce responses shared by up to 60% of all volunteers tested expressing the same class I supertype. This may be due to the inherent focus of our prediction strategy which specifically targets identification of broadly relevant T cell epitopes across all individuals sharing the same supertype. The described tools predict HLA binding potential for a panel of 9 common class II HLA alleles. These alleles were selected because they are common within the human population worldwide and relatively distinct from each other. EpiVax refers to these alleles as supertypes. Each supertype is functionally equivalent or nearly equivalent to many additional family member alleles (i.e., alleles within a supertype family share a set of common peptide binding preferences). Southwood, Sette and Sidney first described the supertype phenomenon for both class I and class II HLA (21, 44, 45). Taken collectively, the 9 class II supertype alleles, along with their respective family members, cover well over 95% of HLA types present in most human population groups. Similar population coverage was determined for class I. Epitopes were identified based on predicted binding to a panel of six common class I alleles: A*0101, A*0201, A*0301, A*2402, B*0702, and B*4403. These alleles are supertypes. Each one is functionally equivalent to or nearly equivalent to many additional family member alleles. Taken collectively, these 6 supertype alleles, along with their respective family members, cover well over 98% of the human population, and some (such as HLA-A02) are expressed by nearly half of the human population. T cell epitopes identified using an overlapping approach may be restricted by rare human class I individual alleles and not broadly across all alleles of a supertype and would thus not be as relevant for development of vaccines for protection of diverse populations. Similarly, the panDR T cell epitope clusters studied here are likely to be broadly relevant, as most of the T cell clusters studied induced responses in multiple volunteers. We note that both strategies (overlapping vs. predictions) have been utilized to successfully identify immunogenic epitopes (46–49). In fact, we were able to identify the same immunodominant T cell epitope using a single participant’s PBMC with both the overlapping peptide/matrix-based deconvolution approach (VAHYMYLIPGLQAAA) and the immunoinformatic predictive strategy described here (A24-7: HYMYLIPGL) (data not shown). However, our direct comparisons in this study demonstrate that the predictive algorithms may outperform overlapping peptide pools specifically for identifying more epitopes broadly relevant across the population.

Several assay methods are used by investigators to identify human MHC-restricted T cell epitopes. Direct ex vivo assays using human PBMCs, such as the ones utilized in this study, identify medium to high frequency effector T cells, and may favor the detection of immunodominant T cell responses over subdominant responses. In vitro assays including measurements of T cell expansion after activation (e.g., cultured ELISPOT assays) allow for detection of T cell responses capable of both proliferation and effector functions including central memory T cells, which after activation and proliferation differentiate into effector T cells. Thus, these expanded ELISPOT assays generally identify higher numbers of T cell epitopes (50, 51). However, it is unknown whether direct ex vivo or expanded T cell assays are better at predicting T cell responses capable of providing long-term protective immunity. In fact, one of the limitations of the current study is that most of the epitopes identified in our work have not been confirmed to be involved in protective immunity. Many investigators utilize HLA transgenic mice for identification of T cell epitopes, which also offer the opportunity to evaluate in vivo protective immunity for certain pathogens (19, 49, 52). Some of the epitopes identified here have been previously reported to be targets of T cells in humans and/or transgenic mice expressing human MHC. For example, ZIKV peptides B7-4, B7-6, and B7-7 identified here (Figure 5C) have been shown to be immunogenic after ZIKV infection in HLA-B7 transgenic mice (19). Furthermore, peptide-based vaccines which included two of these epitopes were shown to provide some protective immunity (significant reduction of viremia) against live ZIKV challenge in HLA-B7 transgenic animals (19). The additional epitopes identified in our study could allow development of protective ZIKV vaccines for a higher proportion of humans expressing highly diverse MHC class I and II molecules.

Our results demonstrate different patterns of class I immunodominance at the supertype level, with different volunteers expressing the same supertypes having a different hierarchy of responses to individual supertype peptides (Figure 2). However, we also observed high frequency immunodominance at the population level across participants expressing certain supertypes when identifying individual epitopes from immunogenic MHC I supertype peptide pools (Figure 5). For example, in 9 participants tested for responses to individual B7 peptides, only 5 epitopes induced a response in at least a single participant (Figure 5C). However, PBMC from 6 of 9 HLA-B7-expressing participants responded to 2 of the B07 peptides (B7-06: NS5₂₉₁₁ RPRVCTKEEF and B7-07: NS2A₁₂₉₀ VPRTDNITL). Similar immunodominance patterns were identified for additional class I supertypes including HLA-A1, -A2, and -B44 (Figures 5A, B, E). It is not known whether this immunodominance is due to different levels of antigen expression, specific fragment proteolysis and processing, or more avid binding of the peptides to MHC or TCR. Additional studies are needed to determine the reasons for such immunodominance.

Responses to the promiscuous ZIKV panDR clusters were more diverse (Figure 5F) than those observed for class I. However, this was expected because each of the ZIKV clusters were predicted to contain epitopes restricted by multiple commonly expressed HLA-DR molecules (ZIKV clusters were identified for promiscuous binding to 9 common HLA-DR supertype alleles). Some of these clusters induced responses in only one or two individuals (of 34 tested), while several induced responses in >25% of participants (e.g., cluster DR3 with a >40% response rate). In total, T cell responses to 27 of the tested 33 ZIKV DR clusters were detected in at least a single volunteer. Eleven ZIKV DR clusters induced responses in >10% of the tested participants. These clusters are highly relevant for development of T cell-targeted vaccines for use in genetically diverse human populations.

Multiple ZIKV vaccine candidates are in various stages of pre-clinical and clinical development and most of these are designed to generate humoral immunity. Vaccine trials using purified formalin-inactivated ZIKV (ZPIV and TAK-426) have demonstrated induction of ZIKV neutralizing antibody responses (8, 53), though durability and magnitude of responses were found to be suboptimal in at least one report (54). Additional vaccine strategies including DNA (55, 56) and adenoviral vectors (57) expressing preM-ENV have also been evaluated in clinical trials and were observed to also induce anti-ZIKV antibody responses. The majority of T cell epitopes identified in this study are present in proteins not expressed in the most advanced recombinant vaccines studied in humans, since our methods allowed identification of epitopes from both structural and non-structural proteins. Our strategy of identifying immunodominant T cell epitopes should provide critical information towards the development of T cell-targeted ZIKV vaccines that invoke immunity in populations expressing diverse HLA. Future studies should be conducted to engineer multi-epitope vaccines expressing these ZIKV-conserved T cell epitopes devoid of homology to human proteins.

References

• 1

  BaudDGublerDJSchaubBLanteriMCMussoD. An update on Zika virus infection. Lancet (2017) 390(10107):2099–109. doi: 10.1016/S0140-6736(17)31450-2

  • CrossRef
  • Google Scholar

• 2

  AubryMTeissierAHuartMMerceronSVanhomwegenJRocheCet al. Zika virus seroprevalence, French Polynesia, 2014-2015. Emerg Infect Dis (2017) 23(4):669–72. doi: 10.3201/eid2304.161549

  • CrossRef
  • Google Scholar

• 3

  Cao-LormeauVMBlakeAMonsSLastereSRocheCVanhomwegenJet al. Guillain-Barre Syndrome outbreak associated with Zika virus infection in French Polynesia: a case-control study. Lancet (2016) 387(10027):1531–9. doi: 10.1016/S0140-6736(16)00562-6

  • CrossRef
  • Google Scholar

• 4

  AntoniouEOrovouESarellaAIliadouMRigasNPalaskaEet al. Zika virus and the risk of developing microcephaly in infants: A systematic review. Int J Environ Res Public Health (2020) 17(11). doi: 10.3390/ijerph17113806

  • CrossRef
  • Google Scholar

• 5

  RiceMEGalangRRRothNMEllingtonSRMooreCAValencia-PradoMet al. Vital signs: zika-associated birth defects and neurodevelopmental abnorMalities possibly associated with congenital zika virus infection - U.S. Territories and freely associated states, 2018. MMWR Morb Mortal Wkly Rep (2018) 67(31):858–67. doi: 10.15585/mmwr.mm6731e1

  • CrossRef
  • Google Scholar

• 6

  El SahlyHMGorchakovRLaiLNatrajanMSPatelSMAtmarRLet al. Clinical, virologic, and immunologic characteristics of zika virus infection in a cohort of US patients: prolonged RNA detection in whole blood. Open Forum Infect Dis (2019) 6(1):ofy352. doi: 10.1093/ofid/ofy352

  • CrossRef
  • Google Scholar

• 7

  PriyamvadaLQuickeKMHudsonWHOnlamoonNSewatanonJEdupugantiSet al. Human antibody responses after dengue virus infection are highly cross-reactive to Zika virus. Proc Natl Acad Sci U.S.A. (2016) 113(28):7852–7. doi: 10.1073/pnas.1607931113

  • CrossRef
  • Google Scholar

• 8

  ModjarradKLinLGeorgeSLStephensonKEEckelsKHde la BarreraRAet al. Preliminary aggregate safety and immunogenicity results from three trials of a purified inactivated Zika virus vaccine candidate: phase 1, randomised, double-blind, placebo-controlled clinical trials. Lancet (2018) 391(10120):563–71. doi: 10.1016/S0140-6736(17)33106-9

  • CrossRef
  • Google Scholar

• 9

  DejnirattisaiWSupasaPWongwiwatWRouvinskiABarba-SpaethGDuangchindaTet al. Dengue virus sero-cross-reactivity drives antibody-dependent enhancement of infection with zika virus. Nat Immunol (2016) 17(9):1102–8. doi: 10.1038/ni.3515

  • CrossRef
  • Google Scholar

• 10

  GargHYehRWattsDMMehmetoglu-GurbuzTResendesRParsonsBet al. Enhancement of Zika virus infection by antibodies from West Nile virus seropositive individuals with no history of clinical infection. BMC Immunol (2021) 22(1):5. doi: 10.1186/s12865-020-00389-2

  • CrossRef
  • Google Scholar

• 11

  TerzianACBSchanoskiASMotaMTOda SilvaRAEstofoleteCFColomboTEet al. Viral load and cytokine response profile does not support antibody-dependent enhancement in dengue-primed zika virus-infected patients. Clin Infect Dis (2017) 65(8):1260–5. doi: 10.1093/cid/cix558

  • CrossRef
  • Google Scholar

• 12

  WenJWangYTValentineKMDos Santos AlvesRPXuZRegla-NavaJAet al. CD4(+) T cells cross-reactive with dengue and zika viruses protect against zika virus infection. Cell Rep (2020) 31(4):107566. doi: 10.1016/j.celrep.2020.107566

  • CrossRef
  • Google Scholar

• 13

  HassertMWolfKJSchwetyeKEDiPaoloRJBrienJDPintoAK. CD4+T cells mediate protection against Zika associated severe disease in a mouse model of infection. PloS Pathog (2018) 14(9):e1007237. doi: 10.1371/journal.ppat.1007237

  • CrossRef
  • Google Scholar

• 14

  Elong NgonoASyedTNguyenAVRegla-NavaJASusantonoMSpasovaDet al. CD8(+) T cells mediate protection against Zika virus induced by an NS3-based vaccine. Sci Adv (2020) 6(45). doi: 10.1126/sciadv.abb2154

  • CrossRef
  • Google Scholar

• 15

  HassertMHarrisMGBrienJDPintoAK. Identification of protective CD8 T cell responses in a mouse model of zika virus infection. Front Immunol (2019) 10:1678. doi: 10.3389/fimmu.2019.01678

  • CrossRef
  • Google Scholar

• 16

  Regla-NavaJAElong NgonoAViramontesKMHuynhATWangYTNguyenATet al. Cross-reactive Dengue virus-specific CD8(+) T cells protect against Zika virus during pregnancy. Nat Commun (2018) 9(1):3042. doi: 10.1038/s41467-018-05458-0

  • CrossRef
  • Google Scholar

• 17

  LiJGaoNFanDChenHShengZFuSet al. Cross-protection induced by Japanese encephalitis vaccines against different genotypes of Dengue viruses in mice. Sci Rep (2016) 6:19953. doi: 10.1038/srep19953

  • CrossRef
  • Google Scholar

• 18

  TarbeMDongWHuGXuYSunJGrayoSet al. Japanese encephalitis virus vaccination elicits cross-reactive HLA-class I-restricted CD8 T cell response against zika virus infection. Front Immunol (2020) 11:577546. doi: 10.3389/fimmu.2020.577546

  • CrossRef
  • Google Scholar

• 19

  WenJTangWWSheetsNEllisonJSetteAKimKet al. Identification of Zika virus epitopes reveals immunodominant and protective roles for dengue virus cross-reactive CD8(+) T cells. Nat Microbiol (2017) 2:17036. doi: 10.1038/nmicrobiol.2017.36

  • CrossRef
  • Google Scholar

• 20

  SidneyJPetersBFrahmNBranderCSetteA. HLA class I supertypes: a revised and updated classification. BMC Immunol (2008) 9:1. doi: 10.1186/1471-2172-9-1

  • CrossRef
  • Google Scholar

• 21

  SouthwoodSSidneyJKondoAdel GuercioMFAppellaEHoffmanSet al. Several common HLA-DR types share largely overlapping peptide binding repertoires. J Immunol (1998) 160(7):3363–73. doi: 10.4049/jimmunol.160.7.3363

  • CrossRef
  • Google Scholar

• 22

  De GrootASArditoMMcClaineEMMoiseLMartinWD. Immunoinformatic comparison of T-cell epitopes contained in novel swine-origin influenza A (H1N1) virus with epitopes in 2008-2009 conventional influenza vaccine. Vaccine (2009) 27(42):5740–7. doi: 10.1016/j.vaccine.2009.07.040

  • CrossRef
  • Google Scholar

• 23

  De GrootASMartinW. Reducing risk, improving outcomes: bioengineering less immunogenic protein therapeutics. Clin Immunol (Orlando Fla) (2009) 131(2):189–201. doi: 10.1016/j.clim.2009.01.009

  • CrossRef
  • Google Scholar

• 24

  MoiseLMcMurryJABuusSFreySMartinWDDe GrootAS. In silico-accelerated identification of conserved and immunogenic variola/vaccinia T-cell epitopes. Vaccine (2009) 27(46):6471–9. doi: 10.1016/j.vaccine.2009.06.018

  • CrossRef
  • Google Scholar

• 25

  MeisterGERobertsCGBerzofskyJADe GrootAS. Two novel T cell epitope prediction algorithms based on MHC-binding motifs; comparison of predicted and published epitopes from Mycobacterium tuberculosis and HIV protein sequences. Vaccine (1995) 13(6):581–91. doi: 10.1016/0264-410X(94)00014-E

  • CrossRef
  • Google Scholar

• 26

  MoiseLGutierrezAKibriaFMartinRTassoneRLiuRet al. iVAX: An integrated toolkit for the selection and optimization of antigens and the design of epitope-driven vaccines. Hum Vaccines immunotherapeutics (2015) 11(9):2312–21. doi: 10.1080/21645515.2015.1061159

  • CrossRef
  • Google Scholar

• 27

  De GrootASBosmaAChinaiNFrostJJesdaleBMGonzalezMAet al. From genome to vaccine: in silico predictions, ex vivo verification. Vaccine (2001) 19(31):4385–95. doi: 10.1016/S0264-410X(01)00145-1

  • CrossRef
  • Google Scholar

• 28

  LiuRMoiseLTassoneRGutierrezAHTerryFESangareKet al. H7N9 T-cell epitopes that mimic human sequences are less immunogenic and may induce Treg-mediated tolerance. Hum Vaccines immunotherapeutics (2015) 11(9):2241–52. doi: 10.1080/21645515.2015.1052197

  • CrossRef
  • Google Scholar

• 29

  LosikoffPTMishraSTerryFGutierrezAArditoMTFastLet al. HCV epitope, homologous to multiple human protein sequences, induces a regulatory T cell response in infected patients. J Hepatol (2015) 62(1):48–55. doi: 10.1016/j.jhep.2014.08.026

  • CrossRef
  • Google Scholar

• 30

  WadaYNithichanonANobusawaEMoiseLMartinWDYamamotoNet al. A humanized mouse model identifies key amino acids for low immunogenicity of H7N9 vaccines. Sci Rep (2017) 7(1):1283. doi: 10.1038/s41598-017-01372-5

  • CrossRef
  • Google Scholar

• 31

  ZhangQWangPKimYHaste-AndersenPBeaverJBournePEet al. Immune epitope database analysis resource (IEDB-AR). Nucleic Acids Res (2008) 36(Web Server issue):W513–8. doi: 10.1093/nar/gkn254

  • CrossRef
  • Google Scholar

• 32

  BrienJDLazearHMDiamondMS. Propagation, quantification, detection, and storage of West Nile virus. Curr Protoc Microbiol (2013) 31:15D 3 1–D 3 8. doi: 10.1002/9780471729259.mc15d03s31

  • CrossRef
  • Google Scholar

• 33

  SchaferJRJesdaleBMGeorgeJAKouttabNMDe GrootAS. Prediction of well-conserved HIV-1 ligands using a matrix-based algorithm, EpiMatrix. Vaccine (1998) 16(19):1880–4. doi: 10.1016/S0264-410X(98)00173-X

  • CrossRef
  • Google Scholar

• 34

  AndradePGimblet-OchiengCModirianFCollinsMCardenasMKatzelnickLCet al. Impact of pre-existing dengue immunity on human antibody and memory B cell responses to Zika. Nat Commun (2019) 10(1):938. doi: 10.1038/s41467-019-08845-3

  • CrossRef
  • Google Scholar

• 35

  BhaumikSKPriyamvadaLKauffmanRCLaiLNatrajanMSChoAet al. Pre-existing dengue immunity drives a DENV-biased plasmablast response in ZIKV-infected patient. Viruses (2018) 11(1). doi: 10.3390/v11010019

  • CrossRef
  • Google Scholar

• 36

  CollinsMHTuHAGimblet-OchiengCLiouGAJadiRSMetzSWet al. Human antibody response to Zika targets type-specific quaternary structure epitopes. JCI Insight (2019) 4(8). doi: 10.1172/jci.insight.124588

  • CrossRef
  • Google Scholar

• 37

  SchouestBGrifoniAPhamJMateusJSydneyJBrienJDet al. Pre-existing T cell memory against zika virus. J Virol (2021) 95(12). doi: 10.1128/JVI.00132-21

  • CrossRef
  • Google Scholar

• 38

  CrookeSNOvsyannikovaIGKennedyRBPolandGA. Identification of naturally processed Zika virus peptides by mass spectrometry and validation of memory T cell recall responses in Zika convalescent subjects. PloS One (2021) 16(6):e0252198. doi: 10.1371/journal.pone.0252198

  • CrossRef
  • Google Scholar

• 39

  MapalagamageMWeiskopfDSetteADe SilvaAD. Current understanding of the role of T cells in chikungunya, dengue and zika infections. Viruses (2022) 14(2). doi: 10.3390/v14020242

  • CrossRef
  • Google Scholar

• 40

  PardyRDRicherMJ. Protective to a T: the role of T cells during zika virus infection. Cells (2019) 8(8). doi: 10.3390/cells8080820

  • CrossRef
  • Google Scholar

• 41

  LimaNSRollandMModjarradK. Trautmann L. T cell immunity and zika virus vaccine development. Trends Immunol (2017) 38(8):594–605. doi: 10.1016/j.it.2017.05.004

  • CrossRef
  • Google Scholar

• 42

  GrifoniAVoicHDhandaSKKiddCKBrienJDBuusSet al. T cell responses induced by attenuated flavivirus vaccination are specific and show limited cross-reactivity with other flavivirus species. J Virol (2020) 94(10). doi: 10.1128/JVI.00089-20

  • CrossRef
  • Google Scholar

• 43

  LaiLRouphaelNXuYNatrajanMSBeckAHartMet al. Innate, T-, and B-cell responses in acute human zika patients. Clin Infect Dis (2018) 66(1):1–10. doi: 10.1093/cid/cix732

  • CrossRef
  • Google Scholar

• 44

  BuiHHSidneyJDinhKSouthwoodSNewmanMJSetteA. Predicting population coverage of T-cell epitope-based diagnostics and vaccines. BMC Bioinf (2006) 7:153. doi: 10.1186/1471-2105-7-153

  • CrossRef
  • Google Scholar

• 45

  SetteASidneyJ. Nine major HLA class I supertypes account for the vast preponderance of HLA-A and -B polymorphism. Immunogenetics (1999) 50(3-4):201–12. doi: 10.1007/s002510050594

  • CrossRef
  • Google Scholar

• 46

  LehmannPVSuwansaardMZhangTRoenDRKirchenbaumGAKarulinAYet al. Comprehensive evaluation of the expressed CD8+ T cell epitope space using high-throughput epitope mapping. Front Immunol (2019) 10:655. doi: 10.3389/fimmu.2019.00655

  • CrossRef
  • Google Scholar

• 47

  LehmannAAKirchenbaumGAZhangTRechePALehmannPV. Deconvoluting the T cell response to SARS-coV-2: specificity versus chance and cognate cross-reactivity. Front Immunol (2021) 12:635942. doi: 10.3389/fimmu.2021.635942

  • CrossRef
  • Google Scholar

• 48

  GrifoniAWeiskopfDRamirezSIMateusJDanJMModerbacherCRet al. Targets of T cell responses to SARS-coV-2 coronavirus in humans with COVID-19 disease and unexposed individuals. Cell (2020) 181(7):1489–501 e15. doi: 10.1016/j.cell.2020.05.015

  • CrossRef
  • Google Scholar

• 49

  EickhoffCSTerryFEPengLMezaKASakalaIGVan AartsenDet al. Highly conserved influenza T cell epitopes induce broadly protective immunity. Vaccine (2019) 37(36), 5371–5381. doi: 10.1016/j.vaccine.2019.07.033

  • CrossRef
  • Google Scholar

• 50

  MeyersLMGutierrezAHBoyleCMTerryFMcGonnigalBGSalazarAet al. Highly conserved, non-human-like, and cross-reactive SARS-CoV-2 T cell epitopes for COVID-19 vaccine design and validation. NPJ Vaccines (2021) 6(1):71. doi: 10.1038/s41541-021-00331-6

  • CrossRef
  • Google Scholar

• 51

  MalavigeGNJonesLBlackAPOggGS. Varicella zoster virus glycoprotein E-specific CD4+ T cells show evidence of recent activation and effector differentiation, consistent with frequent exposure to replicative cycle antigens in healthy immune donors. Clin Exp Immunol (2008) 152(3):522–31. doi: 10.1111/j.1365-2249.2008.03633.x

  • CrossRef
  • Google Scholar

• 52

  EickhoffCSVan AartsenDTerryFEMeymandiSKTrainaMMHernandezSet al. An immunoinformatic approach for identification of Trypanosoma cruzi HLA-A2-restricted CD8⁺ T cell epitopes. Hum Vaccines immunotherapeutics (2015) 11(9):2322–8. doi: 10.1080/21645515.2015.1061160

  • CrossRef
  • Google Scholar

• 53

  HanHHDiazCAcostaCJLiuMBorkowskiA. Safety and immunogenicity of a purified inactivated Zika virus vaccine candidate in healthy adults: an observer-blind, randomised, phase 1 trial. Lancet Infect Dis (2021) 21(9):1282–92. doi: 10.1016/S1473-3099(20)30733-7

  • CrossRef
  • Google Scholar

• 54

  StephensonKETanCSWalshSRHaleAAnselJLKanjilalDGet al. Safety and immunogenicity of a Zika purified inactivated virus vaccine given via standard, accelerated, or shortened schedules: a single-centre, double-blind, sequential-group, randomised, placebo-controlled, phase 1 trial. Lancet Infect Dis (2020) 20(9):1061–70. doi: 10.1016/S1473-3099(20)30085-2

  • CrossRef
  • Google Scholar

• 55

  TebasPRobertsCCMuthumaniKReuschelELKudchodkarSBZaidiFIet al. Safety and immunogenicity of an anti-zika virus DNA vaccine. N Engl J Med (2021) 385(12):e35. doi: 10.1056/NEJMoa1708120

  • CrossRef
  • Google Scholar

• 56

  GaudinskiMRHouserKVMorabitoKMHuZYamshchikovGRothwellRSet al. Safety, tolerability, and immunogenicity of two Zika virus DNA vaccine candidates in healthy adults: randomised, open-label, phase 1 clinical trials. Lancet (2018) 391(10120):552–62. doi: 10.1016/S0140-6736(17)33105-7

  • CrossRef
  • Google Scholar

• 57

  SalischNCStephensonKEWilliamsKCoxFvan der FitsLHeerweghDet al. A double-blind, randomized, placebo-controlled phase 1 study of ad26.ZIKV.001, an ad26-vectored anti-zika virus vaccine. Ann Intern Med (2021) 174(5):585–94. doi: 10.7326/M20-5306

  • CrossRef
  • Google Scholar