
==== Front
JCI Insight
JCI Insight
JCI Insight
JCI Insight
2379-3708
American Society for Clinical Investigation

38861490
180667
10.1172/jci.insight.180667
Research Article
Preservation of naive-phenotype CD4+ T cells after vaccination contributes to durable immunity
Pan Yi-Gen turtle1009@gmail.com
1
Bartolo Laurent laurent.cqfd@gmail.com
1
Xu Ruozhang tammyxrz@gmail.com
12
Patel Bijal V. bijp101@gmail.com
12
Zarnitsyna Veronika I. veronika.i.zarnitsyna@emory.edu
3
Su Laura F. laurasu@upenn.edu
12
1 Department of Medicine, Division of Rheumatology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, USA.
2 Corporal Michael J. Crescenz VA Medical Center, Philadelphia, Pennsylvania, USA.
3 Department of Microbiology and Immunology, Emory University, Atlanta, Georgia, USA.
Address correspondence to: Laura F. Su, University of Pennsylvania, 421 Curie Blvd. BRBII/III 311, Philadelphia, Pennsylvania 19104, USA. Phone: 215.898.4181; Email: Laurasu@upenn.edu.
11 6 2024
22 7 2024
11 6 2024
9 14 e18066728 2 2024
5 6 2024
© 2024 Pan et al.
2024
Pan et al.
https://creativecommons.org/licenses/by/4.0/ This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
This article is available online at https://insight.jci.org/articles/view/180667
Memory T cells are conventionally associated with durable recall responses. In our longitudinal analyses of CD4+ T cell responses to the yellow fever virus (YFV) vaccine by peptide-MHC tetramers, we unexpectedly found CD45RO–CCR7+ virus-specific CD4+ T cells that expanded shortly after vaccination and persisted months to years after immunization. Further phenotypic analyses revealed the presence of stem cell–like memory T cells within this subset. In addition, after vaccination T cells lacking known memory markers and functionally resembling genuine naive T cells were identified, referred to herein as marker-negative T (TMN) cells. Single-cell TCR sequencing detected expanded clonotypes within the TMN subset and identified TMN TCRs shared with memory and effector T cells. Longitudinal tracking of YFV-specific responses over subsequent years revealed superior stability of TMN cells, which correlated with the longevity of the overall tetramer+ population. These findings uncover additional complexity within the post-immune T cell compartment and implicate TMN cells in durable immune responses.

Naive CD4+ T cells persist after vaccination and may contribute to durable protective immunity.

Immunology
Vaccines
Memory
T cells
National Institute of HealthR01AI66358 National Institute of HealthR01AI134879 U.S. Department of Veterans Affairs https://doi.org/10.13039/100000738 I01CX001460 U.S. Department of Veterans Affairs https://doi.org/10.13039/100000738 I01BX005422 Additional funding sources
==== Body
pmcIntroduction

Functional immunological memory underlies the protective efficacy of vaccines against subsequent infections (1, 2). However, the reason protection from some vaccines lasts decades whereas protection from others wanes after a few months remains unknown. A crucial aspect of immune memory involves CD4+ T cells (3). CD4+ T cells provide key signals for B cell maturation and high-affinity antibody production (4). They are also needed to support the expansion and maintenance of functional CD8+ T cells and can directly contribute to antiviral effects (4–6). Past studies in mice and humans have identified naive-like antigen-experienced T cells with superior longevity and plasticity as a source of durable memory (7–9). Broadly categorized as stem cell–like memory T (TSCM) cells, these cells phenotypically resemble naive T cells by positive CCR7 and CD45RA or negative CD45RO expression, yet they display differentiation markers such as CD95, CXCR3, and CD49d (9–11). In people immunized with the highly efficacious and durable yellow fever virus (YFV) vaccine, class I tetramer analyses identified TSCM as the predominant phenotype of virus-specific CD8+ T cells greater than 8 years after vaccination (10, 12).

The durability of CD4+ T cell memory is less understood. Although capable of differentiating into TSCM cells (13–16), CD4+ T cells are generally less responsive to homeostatic cytokines IL-7 and IL-15 (17–19), which augment TSCM differentiation in cultured CD8+ T cells (20). Here, we examined virus-specific CD4+ T cells after YFV vaccination to delineate key features of durable CD4+ T cell responses. YFV-specific CD4+ T cells were identified and tracked longitudinally by direct ex vivo class II peptide-MHC (pMHC) tetramer staining. We showed the presence of various memory subsets within YFV-specific CD4+ T cells several months after vaccination, including TSCM cells. Unexpectedly, about a quarter of CD45RO–CCR7+ tetramer-labeled T cells lacked CD95, CXCR3, CD11a, and CD49d expression, distinguishing them from TSCM cells. These marker-negative T (TMN) cells were a part of expanded clonotypes and shared TCR sequences with memory and effector T cells, suggesting in vivo antigen responses. TMN cells showed minimal decay over years and correlated with the stability of tetramer+ populations. Our findings expand the current definition of antigen-experienced T cells to include those that retain an undifferentiated phenotype. These TMN cells may contribute to durable immunity.

Results

Detection of naive-like CD4+ T cells after YFV vaccination.

We had previously performed a longitudinal study of YFV-specific CD4+ T cells to evaluate the impact of the preexisting repertoire on T cell responses to primary immunization with the YFV vaccine (21). Starting with this data set, we examined features of memory T cells that developed at least 7 months after vaccination. This showed that approximately half of the YFV-specific memory pool consisted of central memory T (TCM) cells, with about 21% of tetramer+ cells retaining a naive-like CD45RO–CCR7+ phenotype (Figure 1, A and B). Proportionally, the abundance of CD45RO–CCR7+ subset was highest before vaccination, decreased initially after vaccination, and then reaccumulated several months later (Figure 1C). Quantified as cells per million CD4+ T cells, CD45RO–CCR7+ tetramer+ cells quickly increased and reached a peak approximately 1 month after vaccination (Figure 1D). The frequency of the CD45RO–CCR7+ T cell subset did not differ significantly by donor age and was associated with the robustness of the response (Supplemental Figure 1, A–C; supplemental material available online with this article; https://doi.org/10.1172/jci.insight.180667DS1). At the memory time point, CD45RO–CCR7+ YFV-specific T cells were more abundant in tetramer+ populations that reached a higher frequency and correlated with efficient recruitment of memory cells (Figure 1, E and F). These data suggest that the presence of CD45RO–CCR7+ CD4+ T cells is a feature of an effective T cell response.

Post-immune T cells are heterogeneous and include a differentiation marker–negative subset.

We hypothesized that the post-vaccine CD45RO–CCR7+ subset largely consisted of TSCM cells as in CD8+ T cells (10, 12). To test this, we performed tetramer staining on 28 YFV-specific CD4+ populations from 7 individuals, recognizing 16 unique epitopes with antibodies against TSCM-associated markers, CXCR3, CD95, CD11a, and CD49d (Supplemental Tables 1 and 2). Staining with this broader antibody panel on blood collected 7 to 48 months after vaccination indeed identified CD45RO–CCR7+tetramer+ T cells that expressed 1 or more TSCM markers. However, we noted that a portion of CD45RO–CCR7+ CD4+ T cells remained negative for CXCR3, CD95, CD11a, and CD49d expression (Figure 2A and Supplemental Figure 2, A–D). To gain further insights into the heterogeneity within the CD45RO–CCR7+ subset, we combined 1,465 YFV-specific CD4+ T cells from 1 donor and visualized combinatorial antibody staining on uniform manifold approximation and projection (UMAP) using the Spectre pipeline (22). This identified regions with low CD45RO and high CCR7 signals, which encompassed a CXCR3+ (cluster 0) and a TSCM marker–negative population (cluster 4) (Figure 2, B–D). We defined CD45RO–CCR7+ cells lacking any measured differentiation markers as TMN cells and classified those expressing at least 1 of CXCR3, CD95, CD11a, or CD49d as TSCM cells (Supplemental Table 3). We then performed manual gating and used Boolean combinations of these gated cells to subdivide the CD45RO–CCR7+ subset into TMN and various TSCM combinations (Supplemental Figure 2, B–D). On average, 27% of the CD45RO–CCR7+ subset consisted of TMN cells (Figure 2E). Among TSCM cells, approximately a quarter expressed only 1 differentiation marker, with CXCR3 being the most common (Figure 2, E and F). Finding antigen-specific T cells that do not express known memory or TSCM markers after a clear prior exposure was unexpected. To test if TMN cells functionally behave like antigen-experienced T cells despite lacking surface markers of differentiation, we treated post-vaccine PBMCs with phorbol myristate acetate (PMA) and ionomycin for 4 to 5 hours. Antigen-specific T cells were captured by tetramers, divided into distinct phenotypic subsets, and analyzed for TNF-α and IFN-γ production. This showed that the post-immune TMN subset produced significantly fewer cytokines compared with memory T cells within the same tetramer+ population (Figure 2, G and H). Thus, YFV vaccination induced a diverse post-immune repertoire that included CD4+ TSCM cells and a naive-like TMN population that lacked phenotypic and functional features of antigen experience.

Virus-specific TMN cells respond to antigens.

We were intrigued by the existence of virus-specific T cells that retained a naive functional phenotype after vaccination. Past studies have identified nonstimulatory TCR interactions that decoupled T cell activation from ligand binding (23). The impaired ability to respond productively to antigens may be one reason why some tetramer-labeled T cells retained a naive phenotype. To investigate this possibility, we quantified TMN, TSCM, and TCM cells for differences in their functional avidity by peptide stimulation. YFV-specific T cell clones were generated using samples from 2 donors obtained 7 to 8 months after YFV vaccination. Among the 48 clones that grew, 40 clones (90%) had the correct specificity by tetramer restaining and/or response to peptides (Figure 3, A and B, and Supplemental Figure 3A). We did not identify peptide-nonresponsive T cells as all clones that were stained with tetramers responded to peptide stimulation. To determine if TMN cells might be harder to activate because of a lower functional avidity, we divided the clones according to their direct ex vivo phenotype and selected 5 clones each from TCM, TSCM, and TMN groups for further analyses. YFV-specific clones were stimulated with decreasing concentrations of the cognate peptide and analyzed for response by cytokine production (Supplemental Figure 3B). TMN cell–derived clones responded similarly to peptides by TNF-α production, with no significant differences in maximal effective peptide concentration (EC50) values between groups (Figure 3, C and D). T cell clones, regardless of their ex vivo phenotypes, also produced similar levels of IFN-γ and IL-2 and had comparable TNF-α+IFN-γ+IL-2+ coexpression (Figure 3E). In addition, we evaluated the proliferative capacity of TMN, TSCM, and TCM cell–derived clones by CellTrace Violet (CTV) dilution and observed no significant differences in the proliferative response to peptide stimulation (Figure 3, F and G, and Supplemental Figure 3C). Thus, TCR-ligand engagement is likely intact for vaccine-specific T cells that retained a naive phenotype after vaccination.

Antigen-experienced TMN cells.

Although TMN cells respond well to antigens in vitro, it remains possible for them to be less competitive in resource-limiting environments. To investigate this, we reason that we can use TCR sequences to infer stimulation and proliferative response in vivo. Because T cell progenies originating from a T cell express identical TCR sequences, we can further leverage these sequences as molecular barcodes to investigate the clonal relationship between distinct phenotypic subsets. However, capturing sufficient numbers of TMN cells was challenging because of their limited number within the available blood samples. To overcome this problem, we generated new tetramers using affinity-matured HLA-DR monomers containing mutations that enhanced CD4 binding to improve the overall capture efficiency (24). When compared with the wild-type (WT) DR, these tetramers stained a larger population of T cells without significantly skewing the phenotypic proportions (Supplemental Figure 4, A–C). In total, we sorted single cells from 5 tetramer-labeled populations and obtained TCR sequences from 607 YFV-specific CD4+ T cells after amplification and sequencing (Figure 4A and Supplemental Table 4). Consistent with clonal expansion after vaccination, over 70% of the sequences were identified in multiple tetramer-sorted T cells. Among expanded sequences, 25% to 52% were abundant and found in at least 10 individual T cells (Figure 4B). Most T cells displayed a TCM or TEM phenotype based on antibody staining at the time of sorting. TMN phenotype was infrequent, expressed by 3% to 4% of sequenced T cells and confined to the 2 most extensively sequenced populations recognizing YF45. Consistent with in vivo expansion, TMN cells did not preferentially express unique TCRs, but rather, they were distributed across various clone sizes (Figure 4C). We focused the subsequent analyses on YF45-specific T cells that included the TMN subset. Early post-vaccine measurements of YF45-specific T cells from HD2 and HD3 showed that both populations had generated robust responses to the YFV vaccine (Figure 4D) (21). In agreement with an antigen-driven response, TMN cells contained expanded clonotypes and shared overlapping sequences with various memory subsets (Figure 4, E and F, and Supplemental Figure 4E). In separately generated T cell clones from the same individuals, TMN cell–derived clones expressed TCRs that matched the sequences from ex vivo–sorted T cells of diverse clone sizes and phenotypes (Supplemental Figure 4F).

The presence of shared TCR sequences with memory T cells, together with clonal expansion, suggest that TMN cells had encountered and responded to antigens. Alternatively, clonotype sharing between memory and naive cells may be explained by the presence of multiple naive T cells with the same TCR, where a portion did not encounter YFV antigen and remained naive. To investigate this possibility, we examined the prevaccination repertoire of YF45-specific T cells in these individuals to determine if preexposure T cells expressing TMN-associated TCRs were abundant before vaccination (21). The changes in clonal dynamics were assessed by tetramer staining, sorting, and sequencing of the TCRs of YF45-specific T cells from blood collected 14 days after vaccination. Of note, the sampling depth at day 14 was much shallower because of limited sample availability. Ten million to 20 million PBMCs were used for tetramer staining whereas over 100 million CD4+ T cells were used to capture rare precursor T cells in the prevaccination sample. In total, we examined TCR sequences from 129 precursor T cells and 238 effector T cells (Figure 5A and Supplemental Table 5). Before vaccination, no pre-vaccine TCRs matched TMN-derived TCRs from HD2, and only 1 sequence was identified in HD3. This shared TCR mapped to a unique sequence and not to the expanded preexisting clonotypes in this individual (Figure 5, B and C). Fourteen days following vaccination, matched TCR frequencies increased to 7% in HD2 and 23% in HD3, capturing a total of 7 clonotypes (Figure 5, B and C). The majority of matched clonotypes were expanded and expressed by a variety of differentiated cells (Figure 5D). Among T cells without a TMN match, 2 cells in HD3 were TMN like and lacked differentiation marker expression. One of these cells was found as a part of an expanded clone at day 14 and shared the same TCR with multiple memory T cells in the day 210 sample (Figure 5E and Supplemental Table 5). Additionally, the post-vaccine day 14 blood from this donor contained 3 CD45RO–CCR7+ T cells without additional differentiation marker staining (Supplemental Table 5).

Human T cells stimulated with cytokines can maintain a naive phenotype (25). Our data suggest that this can also occur with antigen-specific responses. To broaden the analyses, we compared the frequency of CD45RO–CCR7+ T cells in various YFV-specific populations before and after YFV vaccination (Figure 5F). This showed an increase in CD45RO–CCR7+tetramer+ cells after exposure, which included diverse ratios of TSCM and TMN cells (Figure 5G). We focused on the TMN subset as they resembled unstimulated naive T cells. If there is no mechanism to replenish the naive repertoire or hold on to a naive phenotype, we expected a decrease in naive cell frequency following an immune response as responding T cells acquired differentiated states. However, contrary to this, the averaged TMN cell frequency after vaccination was comparable to the frequency of naive cells before vaccination (Figure 5H). To examine population-level differences, we divided the post-vaccine TMN frequency of each population by its initial naive T cell frequency before vaccination. This revealed lower TMN frequencies in two-thirds of the populations, while the remaining third showed an increase (Figure 5I). The absence of a reduction in naive T cells did not indicate a lack of responsiveness to the vaccine. Instead, among the 10 populations that gained naive T cells, there was a higher average fold-change in total tetramer+ frequency before and after vaccination (Figure 5J). Collectively, these data suggest the TMN subset comprises T cells that had responded to vaccination.

TMN cells contribute to durable memory.

While memory T cells are essential for generating rapid recall responses, naive T cells are known for their long lifespan and regenerative potential (26, 27). Therefore, we hypothesized that TMN cells promote the durability of antigen-specific responses after vaccination. To test this idea, we analyzed additional time points from 5 donors with longitudinal PBMCs collected up to 6.7 years after YFV vaccination (Figure 6A and Supplemental Figure 5A). Past modeling of cellular turnover suggests that different phenotypic subpopulations undergo separate and distinct in vivo dynamics (28, 29). To evaluate the stability of individual phenotypic subsets, we subdivided 19 YFV-specific populations according to TMN, TSCM, TCM, TEM, and TEMRA phenotypes based on CD45RO, CCR7, CD95, CXCR3, CD11a, and CD49d expression. Their time-dependent change was quantified as a fitted slope using a mixed effects exponential decay model. This revealed different rates of decay between cells in distinct differentiation states. CD4+ TEM cells had the largest negative slope, indicating the greatest decrease over time. In contrast, TMN cells exhibited remarkable stability, with no discernible decline observed during the follow-up period. The stability of the TMN subset significantly surpassed that of other phenotypic subsets, including TSCM and TCM cells, which are typically considered long-lived (Figure 6B).

Next, we examined the decay kinetics of the overall tetramer+ populations. Because some data were generated before switching to modified DR, paired analyses by WT and modified tetramers on the same blood sample were used to generate an equation for normalizing the frequencies across experiments (Supplemental Figure 4D). Among the 5 donors followed longitudinally, 2 received 1 YFV dose as typical for the YFV vaccine, whereas 3 received a second YFV vaccine 7 months to a year after the initial dose (Supplemental Table 6). Revaccination was originally a part of an IRB-approved protocol to study recall response. Although this strategy did not effectively induce acute T cell responses, likely due to YFV neutralization by antibodies generated from the first vaccination, we investigated whether it might still affect long-term T cell dynamics. We divided tetramer+ populations based on donors’ vaccine doses. Consistent with the longevity of YFV vaccine–mediated protection, YFV-specific CD4+ T cells displayed an average half-life (t1/2) of close to 4 years after 1 YFV immunization (Figure 6C). This durable response was further stabilized after reexposure in the 2-dose group (Supplemental Figure 5B).

TMN cell frequencies did not significantly differ between the 1-dose and 2-dose groups (Supplemental Figure 5C). Given the wide variation in TMN cells, we divided tetramer+ populations based on TMN frequency into top and bottom halves (Figure 6D). This showed that populations with more TMN cells were more stable compared with populations in the bottom TMN group (Figure 6E). The TMN frequency within a given virus-specific population positively correlated with the stability of the overall population (Figure 6F). By contrast, we did not find significant differences between high and low groups based on TSCM, TCM, TEM, and TEMRA frequencies (Supplemental Figure 5D). On the phenotypic level, all tetramer+ populations contained various memory subsets, but the top TMN group was more phenotypically diverse. We divided populations based on the first TMN frequency obtained within the 1–2 years after YFV vaccination and showed that populations with more TMN cells exhibited greater diversity of differentiation states over time as measured by the Shannon diversity index (Figure 6, G and H). Collectively, these data highlight the stability of the TMN subset and uncover their association with durable and diverse T cell memory after YFV vaccination.

Discussion

A T cell is typically referred to as naive if it has not yet encountered its specific cognate antigen(s). However, a clear antigenic history is often not available in human studies, and thus specific surface markers are commonly used to infer antigen experience (11, 30). Here, we examined CD4+ T cell memory to YFV vaccination to define key features of durable memory by direct ex vivo class II tetramer staining and enrichment. This showed a diverse memory pool composed of various differentiation states after YFV vaccination. Surprisingly, we also uncovered an increase in naive-phenotype YFV-specific T cells after vaccination.

Deeper phenotypic analyses showed that a portion of CD45RO–CCR7+ cells were previously described TSCM cells. The remaining cells did not express known surface markers of T cell differentiation that were tested. We evaluated several potential explanations for their origin. Our initial idea that TMN cells were unresponsive was incorrect. Peptide-stimulated TMN cell–derived T cell clones produced robust levels of cytokines at comparable EC50 as T cell clones generated from memory cells. The presence of expanded TCR clones within the TMN subset further suggests that some TMN cells had encountered and responded to antigens in vivo. The observation that TMN cells, memory T cells, and effector T cells have overlapping TCR repertoires suggests shared precursors among these subsets and provides additional support to the idea that some TMN cells are antigen experienced, as their common precursors divided at least once in response to antigen. This extent of overlap would not be expected if clonal expansion within the TMN subset were driven by different processes, such as cytokine-mediated and/or bystander response. Our interpretation that some expanded antigen-experienced cells exist within the TMN population is further supported by increased naive T cell frequency following immunization. We show that, despite recruitment into effector responses, the frequency of naive T cell frequency recognizing YFV is largely unchanged. Among populations analyzed before and after immunization, approximately one-third gained more naive-phenotype T cells in the form of TMN cells above their pre-vaccine baseline. It is also unlikely we would capture TMN cells 14 days after vaccination if they had not expanded, since the same clonotype was not detected before vaccination despite analyzing many more cells.

In addition to an antigen-experienced subset, some TMN cells could be genuine naive T cells. As an alternative possibility, there might be multiple naive cells of the same clonotype, but only a subset responded to the vaccine. Incomplete response may have occurred because some YFV-specific T cells were not at the right place or time to respond to the immunogen. Other cells may not have received enough of the appropriate signal to undergo expansion and/or differentiation (31). This model could be consistent with our TCR data but does not explain why certain specificities contain higher frequencies of naive T cells after vaccination. An intriguing alternative for this increase might involve ongoing thymic activities. Although thymic output for new T cells declines with age, recent thymic emigrants remain detectable in adults (32–34). Newly produced T cells might add to the peripheral TMN pool to increase the abundance of YFV-specific naive T cells after vaccination. TMN cells from various sources could contribute to the naive repertoire after antigen exposures, varying in their impact depending on antigen specificity, the type of stimulation, age, and other individual characteristics.

Only a few select vaccines are capable of mediating lifelong protection. How durable immunological memory is maintained remains a key unresolved question. While memory T cells are the cornerstone of protective immunity by virtue of their ability to rapidly initiate a functional response to pathogen rechallenge, naive T cells possess superior self-renewal capacity and differentiation plasticity (3, 26, 27, 35). Irrespective of how TMN cells might have originated, we asked whether TMN cells are related to the longevity of T cell response. Our findings revealed remarkable stability of the TMN population, showing minimal decay for nearly 7 years. These stable T cells could potentially support the longevity of the overall immune response, extending it beyond the lifespan of individual memory T cells. Consistent with this model, TMN cells are more abundant in durable CD4+ populations that are stable over time. Based on the diverse memory phenotypes in TMN-enriched populations, we further speculate that TMN cells have the potential to differentiate into multiple states, thereby contributing to the phenotypic diversity of T cell memory.

A limitation of our study is that we did not further subset TMN cells to evaluate potential heterogeneity in gene expression, epigenetic landscape, or functional attributes because of the limitations on cell numbers. Future studies using models allowing manipulation of TMN cells will be needed to build on their association with durable immunity and establish causation. As our analyses are focused on CD4+ T cell responses to YFV in a few healthy individuals, larger studies involving additional vaccines will be important to strengthen our findings and to understand how the TMN subset influences vaccine durability across different ages and disease settings.

In summary, our findings highlight the complexity within the post-immune T cell compartment and add to our understanding of the diverse spectrum of T cells exhibiting naive features (27, 36). Using an antigen-specific approach, our data suggest that some cells considered naive by phenotypic criteria are actually antigen experienced. As a whole, post-immune TMN cells are remarkably stable. Understanding how TMN cells are generated, are maintained, and work alongside memory cells in providing long-lasting immune protection could aid the future development of improved vaccine strategies.

Methods

Sex as a biological variable

Both male and female participants were included in this study. Study volunteers were recruited in the order of participation without restrictions based on the sex.

Human samples

This study used cryopreserved cells stored in FBS with 10% DMSO from an ongoing vaccine study at the University of Pennsylvania (21). This study includes 7 healthy adult participants with no prior YFV exposure who received 1 or 2 doses of the 17D live-attenuated YFV vaccine (YF-VAX, Sanofi Pasteur). Individuals older than 65 were excluded. As YFV vaccine provides lifelong protection against YFV, revaccination with the attenuated 17D-204 strain is not expected to increase risk. Five participants were followed longitudinally for 2 to 6.7 years after vaccination. Participant characteristics are shown in Supplemental Table 1.

Cell lines

Hi5 cells (Thermo Fisher Scientific) were maintained by insect cell culture medium (ESF921, Expression Systems) supplemented with 0.02% gentamicin at 28°C.

Protein expression and tetramer production

His-tagged HLA-DRA/B1*0301, 0401, 0407, and 1501 protein monomers of WT sequence or with L112W, S118H, V143M, or T157I mutations (24) were produced by Hi5 insect cells and extracted from culture supernatant using Ni-NTA (QIAGEN). HLA-DR monomers were biotinylated overnight at 4°C using BirA biotin ligase (Avidity) and purified by size-exclusion chromatography using Superdex 200 size-exclusion column (AKTA, GE Healthcare, now Cytiva). Biotinylation was confirmed by gel-shift assay. Peptide exchange and tetramerization for WT and modified affinity-matured DR were performed using standard protocols as previously described (37, 38). In brief, HLA-DR proteins were incubated with thrombin (MilliporeSigma) at room temperature for 3–4 hours and exchanged with peptides of interest in 50-fold excess at 37°C for 16 hours. Peptide-loaded HLA-DR monomers were incubated with fluorochrome-conjugated streptavidin at 4–5:1 ratio for 2 minutes at room temperature, followed by a 15-minute incubation with an equal volume of biotin-agarose slurry (MilliporeSigma). Tetramers were buffer exchanged into PBS, concentrated using Amicon ULTRA 0.5 mL 100 kDa (MilliporeSigma), and kept at 4°C for no more than 2 weeks prior to use.

Ex vivo T cell analyses and cell sorting

Phenotypic analyses and frequency quantification.

Tetramer staining was performed on at least 10,000,000 PBMCs with 5 μg of tetramers in a 100 μL reaction for 1 hour at room temperature as previously described (21, 38, 39). Tetramer-tagged cells were enriched by adding anti-fluorochrome and anti-HIS MicroBeads (Miltenyi Biotec). The mixture was passed through LS columns (Miltenyi Biotec). Column-bound cells were washed and eluted according to manufacturer protocol. For antibody staining, the enriched samples were stained with viability dyes and exclusion markers (anti-CD19 and anti-CD11b, Supplemental Table 7), along with combinations of surface markers as specified in the experiments (anti-CD3, anti-CD4, anti-CD45RO, anti-CCR7, anti-CD11a, anti-CD95, anti-CD49d, and anti-CXCR3, Supplemental Table 7), in 50 to 100 μL of FACS buffer (PBS plus 2% FCS, 2.5 mM EDTA, 0.025% sodium azide) for 30 minutes at 4°C. Samples were fixed with 2% paraformaldehyde and acquired by flow cytometry using LSRII (BD). Data analyses were performed by FlowJo (BD). The Boolean tool was used to define TSCM and TMN cells within the CD45RO–CCR7+ subset. Frequency calculation was obtained by mixing one-tenth of samples with 200,000 fluorescence beads (Spherotech) for normalization.

For longitudinal experiments involving both WT and modified DR, paired data from WT and modified DR, with a minimum of 2 data points per time point for each specificity, were used to derive the equation for normalization: log2(Freqmodified) = 3.72 + 0.35 × log2(FreqWT) (Supplemental Figure 4D). Frequencies generated by WT tetramers that were below the normalized values were adjusted. Mixed effects exponential decay models were used to analyze longitudinal changes in antigen-specific T cell populations and estimate the corresponding slopes. These models were implemented in MonolixSuite 2021R1 (Lixoft) and fitted to data after vaccination. Initial T cell specificity values were lognormally distributed, exponential decay rates were normally distributed, and lognormal multiplicative error was used. The estimation of the population parameters was performed using the Stochastic Approximation Expectation-Maximization algorithm. Half-lives were calculated as ln(2)/k, where the corresponding k values represented the estimated exponential decay rate constants. Estimated decay rates were converted into slopes as –k.

For multidimensional analyses, a total of 1,465 manually gated tetramer+ cells were exported from FlowJo, read into R by flowCore, and combined into 1 data set for subsequent data processing and analyses using the Spectre package in R (22). Staining intensities were converted using Arcsinh transformation with a cofactor of 200. Batch alignment was performed using CytoNorm (40). Clustering was performed using Phenograph with nearest neighbors set to 55 (k = 55) (41). UMAP was used for dimensional reduction and visualization (42).

Function response.

T cells were rested overnight, followed by 4–5 hours of stimulation by PMA (5 ng/mL, MilliporeSigma) and ionomycin (500 ng/mL, MilliporeSigma) in the presence of monensin (2 μM, MilliporeSigma) and Brefeldin A (5 μg/mL, MilliporeSigma). Tetramer and surface antibody staining were performed as above. Intracellular staining with antibodies against TNF-α, IFN-γ, CD3, and CD4 (Supplemental Table 7) was performed following BD Cytofix/Cytoperm Fixation/Permeabilization Kit according to manufacturer protocol.

Cell sorting.

Around 60,000,000 CD3+ or CD4+ T cells were used and stained with up to 10 μg of each tetramer in a 100 μL reaction. Antibody staining was performed as above without fixation. Individual tetramer-labeled cells were isolated for TCR sequencing or T cell cloning by index sorting using the purity mode on FACSAria (BD).

Generation and stimulation of T cell clones

Clone generation.

Cells were stained with tetramers and enriched with magnetic beads as described above. Single-tetramer-stained CD4+ T cells were sorted into individual wells in a round-bottom, 96-well plate containing 105 irradiated PBMCs, 104 JY cell line (Thermo Fisher Scientific), PHA (1:100, Thermo Fisher Scientific), IL-7 (25 ng/mL, PeproTech), and IL-15 (25 ng/mL, PeproTech). IL-2 (50 IU/mL, PeproTech) was added on day 5 and replenished every 3–5 days. Cells were resupplied with fresh medium with IL-2 (50 IU/mL), PHA (1:100), and 105 irradiated PBMCs every 2 weeks.

DC generation.

Monocytes from HLA-DR allele-matched donors were isolated using negative enrichment kits (RosetteSep Human Monocyte Enrichment Cocktail, StemCell Technology). A total of 10,000,000 cryopreserved monocytes were cultured in 15 mL DC media (RPMI 1640 plus glutamine, 10% FCS, 1× penicillin/streptomycin, 10 mM HEPES) in the presence of 100 ng/mL GM-CSF and 500 U/mL IL-4. Three days later, half the culture media was replaced with fresh DC media with 100 ng/mL GM-CSF, 500 U/mL IL-4, and 0.05 mM 2-mercaptoethanol. Cells in suspension were harvested at 5 to 6 days and added to a flat-bottom, 96-well plate at 25,000 DCs per well. DCs were treated with 100 ng LPS and peptides (0.00001 μg/mL to 10 μg/mL) for 16 hours and replenished with fresh media before coculturing with T cells.

Stimulation of T cell clones.

T cell clones were rested overnight in fresh media without IL-2 and added to wells containing matured DCs at 1:1 ratio in the presence of monensin (2 μM, MilliporeSigma) and Brefeldin A (5 μg/mL, MilliporeSigma). After 5 hours, cells were transferred into a new 96-well, round-bottom plate, washed once with FACS buffer, and stained with viability dyes and exclusion markers (anti-CD19 and anti-CD11b) for 30 minutes at 4°C. Intracellular staining with antibodies against TNF-α, IFN-γ, IL-2, CD3, and CD4 (Supplemental Table 7) was performed following BD Cytofix/Cytoperm Fixation/Permeabilization Kit according to manufacturer protocol. EC50 was determined using the percentage of T cell clones that produced TNF-α in response to decreasing peptide concentrations (10, 1, 0.1, 0.01, 0.001, 0.0001, and 0.00001 μg/mL). A nonlinear fit without constraint was applied to log-transformed concentration using the equation Y = Bottom + (Top – Bottom)/(1 + 10[(LogEC50 – X) × HillSlope]) in Prism (GraphPad). For the proliferation assay, T cell clones were labeled with 1:1,000 diluted CTV (Thermo Fisher Scientific) following manufacturer protocol. The CTV-stained cells were rested in fresh media without IL-2 for 16 hours. A total of 25,000 rested T cells were cocultured with DCs pulsed with 10 μg/mL cognate peptides or treated with PHA as a positive control (1:100, Thermo Fisher Scientific). After 5 days, cells were harvested and stained with viability dyes and surface antibodies (anti-CD19, anti-CD11b, anti-CD3, and anti-CD4 for 30 minutes at 4°C followed by fixation with 2% paraformaldehyde. Samples were acquired by flow cytometry using LSRII and analyzed by FlowJo.

Single-cell TCR sequencing and analyses

Single-cell TCR sequencing by nested PCRs was performed using the primer sets and the protocol as previously described (21, 43). In brief, reverse transcription was performed with CellsDirect One-Step qRT-PCR kit according to the manufacturer’s instructions (CellsDirect, Invitrogen) using a pool of 5′ TRVB region–specific primers and 3′ C region primers. The cDNA library was amplified using a second set of multiple internally nested V-region and C-region primers with HotStarTaq DNA polymerase kit (QIAGEN). The final PCR was performed on an aliquot of the second reaction using a primer containing common base sequence and a third internally nested Cβ primer. PCR products were gel purified (QIAGEN) and sequenced on a NovaSeq 6000 platform (Illumina). TCR sequences were pre-processed as previously described (21). In brief, forward and reverse reads were converted into 1 paired-end read using pandaseq (44). Data were demultiplexed by the unique combination of plate, row, and column barcodes. Consensus TCRβ sequences were identified using the V(D)J alignment software MiXCR (45). A threshold of a read count of 200 reads per sequence was applied to the consensus sequences. If more than 1 TCRα or TCRβ chain passes this criterion, we retain the dominant TCRβ and the 2 TCRα chains with the highest read count. For data obtained from cells several months after vaccination, we additionally require phenotypic annotation based on antibody staining from index sort data. Data were excluded if phenotypic information was not retained or was ambiguous. For downstream analyses, data wrangling was performed using the tidyverse package. TCRs were matched by TCRβ if only the β chain was available or by TCRβ plus at least 1 TCRα if α chain(s) were called. Circos plots were made using the circlize package of R software (46).

Statistics

Normality was assessed using D’Agostino-Pearson test. Spearman’s was used if either of the 2 variables being correlated was non-normal. Otherwise, Pearson’s was used to measure the degree of association. Least squares linear regression was used to calculate the best fitting line. Statistical comparisons were performed using 2-tailed Student’s t test, paired t test, Welch’s 1-way ANOVA, repeated measures 1-way ANOVA, 2-way ANOVA, or mixed effects model. A P value of less than 0.05 was used as the significance level and adjusted using Tukey’s multiple comparisons test or Dunn’s multiple comparison test, as indicated in the figure legends, if multiple comparisons were performed. Statistical analyses were performed using GraphPad Prism. Lines and bars represent the mean, and variability is represented by the SEM.

Study approval

All participants have given written informed consent. All samples were deidentified and obtained with IRB regulatory approval from the University of Pennsylvania. The study was approved by the IRB at the University of Pennsylvania (approval 820884).

Data availability

All data needed to evaluate the conclusions in the paper are present in the paper or the supplement. Analyses are performed using standard analysis packages. All data points are reported in the Supporting Data Values file.

Author contributions

LFS was responsible for conceptualization, YP for experimentation, LB for sequence analyses, RX for high-dimensional phenotypic analyses, BVP for study recruitment, VIZ for modeling and statistical support, and LFS for supervision; LFS, VIZ, LB, and YP were responsible for manuscript preparation.

Supplementary Material

Supplemental data

Supporting data values

We thank our study participants for their participation. Funding sources were NIH R01AI134879 (to LFS), NIH R01AI66358 (to LFS), VA Merit Award I01CX001460 (to LFS), VA COVID Award I01BX005422 (to LFS), and CHU fund (to LFS).

Version 1 06/11/2024

In-Press Preview

Version 2 07/22/2024

Electronic publication

Figure 1 Identification of CD45RO–CCR7+ virus-specific CD4+ T cells after YFV vaccination.

(A) Direct ex vivo tetramer and antibody staining of a representative YFV tetramer+ (tet+) population using blood collected about 7 months after YFV vaccination. (B) The percentage of YFV tet+ T cells with the indicated combination of CD45RO and CCR7 expression. Plot summarizes data from 36 specificities 7 to 34 months after YFV vaccination from 7 donors. (C and D) The abundance of CD45RO–CCR7+YFV tet+ CD4+ T cells in 7 healthy participants was quantified as a percentage of tetramer+ cells (C) or by frequency (D). Each symbol represents data from a distinct YFV-specific population. Experiments were repeated an average of 3.3 times. (E) The correlation between CD45RO–CCR7+YFV tet+ T cell frequency and the overall tet+ frequency of the same population at least 7 months after vaccination. (F) The correlation between the CD45RO–CCR7+YFV tet+ T cell frequency and the fold-change from the pre-vaccine baseline to the memory time point. CD45RO–CCR7+tet+ frequencies were determined using samples taken 7 to 34 months after vaccination. n = 33, and populations with no post-vaccine naive cells were excluded. Repeated measures (RM) 1-way ANOVA (B) or mixed effects analysis (C and D) was performed and corrected with Tukey’s multiple comparisons test. (E and F) Pearson’s correlation was computed. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001.

Figure 2 Post-vaccine CD4+ T cells are heterogeneous and include naive-like subsets.

(A) FACS plots show the expression of the indicated marker on a representative YFV-specific population. The tet+ population is overlaid onto tet– bulk CD4+ T cells. (B) UMAP displays phenograph-defined clusters. Data combine 1,465 CD4+ cells labeled by 7 YFV tetramers from HD3. (C and D) The staining intensity of individual markers is shown on a heatmap for each cluster (C) or displayed on the UMAP (D). (E and F) The relative abundance of CD45RO–CCR7+ YFV-specific T cells by the indicated numbers of markers (E) or the type of markers (F). Frequency in F combines all cells positive for a particular marker within the CD45RO–CCR7+ subset. Marker combinations were determined using Boolean operators on manually defined gates. Each symbol represents a tetramer+ population (n = 28). Experiments were repeated an average of 2.5 times. (G) PBMCs were stimulated for 4–5 hours by PMA and ionomycin and assayed for cytokine production by intracellular cytokine staining. The plot shows representative TNF-α and IFN-γ expression by TMN cells and memory T cells (non-CD45RO–CD28+) from the same tetramer-labeled population. (H) T cell responses by TNF-α and IFN-γ production for the indicated phenotypic subset. Each population was identified with a pool of 5–7 tetramers of the same DR allele, using cells from 3 donors. (E and F) RM 1-way ANOVA was performed and corrected with Tukey’s multiple comparison test. (H) The Friedman test was performed and corrected using Dunn’s multiple comparison test. *P < 0.05, ***P < 0.001, ****P < 0.0001.

Figure 3 TMN-derived T cell clones respond to antigen stimulation.

(A) Schematics of single-cell T cell cloning. Post-vaccine T cells from HD2 and HD3 were stained with 1501-YF45 tetramers; sorted based on TMN, TSCM, or TCM phenotypes; and expanded for 2 to 3 weeks in culture. (B) In vitro–expanded T cell clones were restained with tetramers and cultured with vehicle- (DMSO) or peptide-treated monocyte-derived DCs. Representative plots show tetramer staining and cytokine production by intracellular cytokine staining. (C and D) T cell clones were stimulated with decreasing concentrations of YFV peptides. The response was measured by TNF-α production (C) and quantified by EC50 values after subtracting the background signal from vehicle-treated control (D). (E) Peptide dose response of T cell clones by IFN-γ, IL-2, and IL-2+TNF-α+IFN-γ+ production. (F) Representative histograms show CTV dilution in response to 10 μg/mL of peptide stimulation. (G) Plot summarizes the frequency of the CTVlo population after a 5-day culture for clones in each phenotypic group. All experiments were repeated at least twice with n = 5 in each group. (C and E) RM 2-way ANOVA was performed and corrected with Tukey’s multiple comparison test. (D and G) Kruskal-Wallis test and Dunn’s multiple comparison test were used. PHA, phytohemagglutinin.

Figure 4 TMN cells are clonally related to memory T cells.

(A) The plot summarizes the number of cells sequenced from the indicated specificities and donors collected 242 (HD2) or 210 (HD3) days after primary YFV vaccination. (B) Clone size and phenotypic distribution of YFV tet+ populations in A. Phenotypic data were obtained by index sorting: TMN (CD45RO–CCR7+CXCR3–CD95–CD11a–CD49d–), TSCM (CD45RO–CCR7+ and positive for at least 1 of CXCR3, CD95, CD11a, or CD49d), TCM (CD45RO+CCR7+), TEM (CD45RO+CCR7–), TEMRA (CD45RO–CCR7–). Cells with ambiguous phenotypes were excluded. (C) Distribution of phenotypes in B by clonotype frequency, ranked from largest to unique clonotypes. (D) The frequencies of YF45-specific T cells before vaccination and at the indicated days following YFV immunization. (E) Each circos plot represents TCRs from YF45 tet+ cells obtained 210–242 days after vaccination, separated by the associated indexed phenotypes. Cells are ordered by frequency within each arc. Gray marks cells expressing unique TCRs; other colors represent expanded or shared sequences. Shared TCRβ or TCRα and β, when a TCRα is available, is connected by a line across distinct phenotypic subsets. Blue lines highlight TCRs from TMN cells that are shared with cells expressing other phenotypes. (F) The percentage of TCRs in each memory subset that matched TMN-derived sequences. TEM, effector memory; TEMRA, effector memory cells reexpressing CD45RA.

Figure 5 TMN cells are clonally related to effector T cells and contribute to the maintenance of the naive repertoire after vaccination.

(A) The number of TCRs from YF45 tet+ T cells from the indicated donors, before and 14 days after YFV vaccination. (B) TCRs from cells in A were compared to TMN-derived TCRs in the day 210–242 post-vaccine samples. Bar graphs show the percentages of cells with or without a match. (C) Each circos plot represents TCRs from TMN YF45 tet+ cells and TCRs from the same specificities, before and 14 days after YFV vaccination, from the same donors. Cells are ordered by frequency within each arc. Gray marks cells expressing unique TCRs; other colors represent expanded or shared sequences. Connecting lines highlight shared TCRs between a TMN cell and cells in a previous time point. (D) The number and phenotypes of T cells in the day 14 sample that matched a TMN-derived clonotype from the memory time point. (E) Clonal dynamics of a TCR expressed by a TMN cell in the 14-day post-vaccine sample. Plot shows the number and phenotypes of T cells that expressed the same TCR sequence at the indicated time points. ND, not detected. (F and G) The frequencies of CD45RO–CCR7+ YFV-specific T cells before and after vaccination for the corresponding populations (F). Post-vaccine CD45RO–CCR7+tet+ cells were separated into TSCM or TMN subsets based on CXCR3, CD95, CD11a, or CD49d staining (G). (H) Plot compares the frequencies of CD45RO–CCR7+tet+ cells before vaccination with the frequencies of TMN cells in the corresponding population from the post-vaccine sample. (I) The post-vaccine TMN frequency of each tet+ population was divided by its initial CD45RO–CCR7+ frequency before vaccination. The ratio is defined as a gain if above 1 and a loss if below 1. Pie chart shows the numbers of populations that had gained or lost naive T cells after vaccination. (J) Plot shows the fold-change in tet+ frequency, from the pre-vaccine baseline to the memory time point, for populations that had gained or lost naive cells. For F, H, and J, n = 27 tetramer+ populations. Welch’s t test was performed. *P < 0.05, **P < 0.01.

Figure 6 TMN cells are stable and associated with durable T cell memory.

(A) Representative plots show YFV-specific CD4+ T cells over the indicated time points from HD2. (B) Each tet+ population of a given specificity was subdivided according to phenotypes. The change over time for each phenotypic subset was quantified by the estimated slope using a mixed effects exponential decay model (n = 19 tetramer+ populations from 5 donors). (C) A mixed effects exponential decay model fitted to the dynamics of YFV-specific CD4+ T cells after a single YFV vaccination (n = 8 populations, combined from donors 4 and 5). The estimated decay (blue line) was used for calculating the t1/2. (D) Ranking of tet+ populations by the averaged frequency of TMN cells within each population across all time points. (E) Plot summarizes the estimated slopes of individual tet+ populations, divided into top and bottom halves by TMN frequency in D. (F) The correlation between slopes characterizing the change over time for the overall tetramer+ populations and their corresponding averaged TMN frequencies. (G) Pie charts show the distribution of memory subsets. Populations were divided into top and bottom groups by the first measured TMN frequency obtained within 1–2 years after YFV vaccination. (H) Phenotypic diversity of each tet+ population was quantified using Shannon diversity index, categorized into top or bottom groups based on TMN frequency as in G. Each symbol represents 1 tetramer+ population. Experiments were repeated an average of 2.3 times. Data are represented as mean ± SEM. (B) RM 1-way ANOVA was performed and corrected with Tukey’s multiple comparisons test. (E and H) Welch’s t test was performed. (F) Spearman’s correlation was performed. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001.

Conflict of interest: The authors have declared that no conflict of interest exists.

Copyright: © 2024, Pan et al. This is an open access article published under the terms of the Creative Commons Attribution 4.0 International License.

Reference information: JCI Insight. 2024;9(14):e180667.https://doi.org/10.1172/jci.insight.180667.
==== Refs
1 Kathryn AF et al Preexisting high frequencies of memory CD8+ T cells favor rapid memory differentiation and preservation of proliferative potential upon boosting Immunity 2013 39 1 171 183 10.1016/j.immuni.2013.07.003 23890070
2 Martin JR et al Pathogen-specific inflammatory milieux tune the antigen sensitivity of CD8(+) T cells by enhancing T cell receptor signaling Immunity 2013 38 1 140 152 10.1016/j.immuni.2012.09.017 23260194
3 Wherry EJ Barouch DH T cell immunity to COVID-19 vaccines Science 2022 377 6608 821 822 10.1126/science.add2897 35981045
4 Crotty S T follicular helper cell differentiation, function, and roles in disease Immunity 2014 41 4 529 542 10.1016/j.immuni.2014.10.004 25367570
5 Laidlaw BJ et al The multifaceted role of CD4(+) T cells in CD8(+) T cell memory Nat Rev Immunol 2016 16 2 102 111 10.1038/nri.2015.10 26781939
6 Swain SL et al Expanding roles for CD4+ T cells in immunity to viruses Nat Rev Immunol 2012 12 2 136 148 10.1038/nri3152 22266691
7 Gattinoni L et al A human memory T cell subset with stem cell-like properties Nat Med 2011 17 10 1290 1297 10.1038/nm.2446 21926977
8 Zhang Y et al Host-reactive CD8+ memory stem cells in graft-versus-host disease Nat Med 2005 11 12 1299 1305 10.1038/nm1326 16288282
9 Pulko V et al Human memory T cells with a naive phenotype accumulate with aging and respond to persistent viruses Nat Immunol 2016 17 8 966 975 10.1038/ni.3483 27270402
10 Fuertes Marraco SA et al Long-lasting stem cell-like memory CD8+ T cells with a naïve-like profile upon yellow fever vaccination Sci Transl Med 2015 7 282 282ra48 10.1126/scitranslmed.aaa3700 25855494
11 Caccamo N et al Atypical human effector/memory CD4+ T cells with a naive-like phenotype Front Immunol 2018 9 2832 10.3389/fimmu.2018.02832 30559746
12 Akondy RS et al Origin and differentiation of human memory CD8 T cells after vaccination Nature 2017 552 7685 362 367 10.1038/nature24633 29236685
13 Long HM et al MHC II tetramers visualize human CD4+ T cell responses to Epstein-Barr virus infection and demonstrate atypical kinetics of the nuclear antigen EBNA1 response J Exp Med 2013 210 5 933 949 10.1084/jem.20121437 23569328
14 Buzon MJ et al HIV-1 persistence in CD4+ T cells with stem cell-like properties Nat Med 2014 20 2 139 142 10.1038/nm.3445 24412925
15 Mpande CAM et al Functional, antigen-specific stem cell memory (TSCM) CD4+ T cells are induced by human mycobacterium tuberculosis infection Front Immunol 2018 9 324 10.3389/fimmu.2018.00324 29545791
16 Pawel M et al Th17 cells are long lived and retain a stem cell-like molecular signature Immunity 2011 35 6 972 985 10.1016/j.immuni.2011.09.019 22177921
17 Tan JT et al Interleukin (IL)-15 and IL-7 jointly regulate homeostatic proliferation of memory phenotype CD8+ cells but are not required for memory phenotype CD4+ cells J Exp Med 2002 195 12 1523 1532 10.1084/jem.20020066 12070280
18 Cho JH et al T cell receptor-dependent regulation of lipid rafts controls naive CD8+ T cell homeostasis Immunity 2010 32 2 214 226 10.1016/j.immuni.2009.11.014 20137986
19 Guimond M et al Interleukin 7 signaling in dendritic cells regulates the homeostatic proliferation and niche size of CD4+ T cells Nat Immunol 2009 10 2 149 157 10.1038/ni.1695 19136960
20 Cieri N et al IL-7 and IL-15 instruct the generation of human memory stem T cells from naive precursors Blood 2013 121 4 573 584 10.1182/blood-2012-05-431718 23160470
21 Pan YG et al Vaccination reshapes the virus-specific T cell repertoire in unexposed adults Immunity 2021 54 6 1245 1256 10.1016/j.immuni.2021.04.023 34004140
22 Ashhurst TM et al Integration, exploration, and analysis of high-dimensional single-cell cytometry data using Spectre Cytometry A 2022 101 3 237 253 10.1002/cyto.a.24350 33840138
23 Sibener LV et al Isolation of a structural mechanism for uncoupling T cell receptor signaling from peptide-MHC binding Cell 2018 174 3 672 687 10.1016/j.cell.2018.06.017 30053426
24 Sugata K et al Affinity-matured HLA class II dimers for robust staining of antigen-specific CD4+ T cells Nat Biotechnol 2021 39 8 958 967 10.1038/s41587-021-00836-4 33649568
25 Soares MV et al IL-7-dependent extrathymic expansion of CD45RA+ T cells enables preservation of a naive repertoire J Immunol 1998 161 11 5909 5917 10.4049/jimmunol.161.11.5909 9834071
26 Gattinoni L et al T memory stem cells in health and disease Nat Med 2017 23 1 18 27 10.1038/nm.4241 28060797
27 van den Broek T et al The full spectrum of human naive T cells Nat Rev Immunol 2018 18 6 363 373 10.1038/s41577-018-0001-y 29520044
28 Asquith B et al Lymphocyte kinetics: the interpretation of labelling data Trends Immunol 2002 23 12 596 601 10.1016/S1471-4906(02)02337-2 12464572
29 Gossel G et al Memory CD4 T cell subsets are kinetically heterogeneous and replenished from naive T cells at high levels Elife 2017 6 e23013 10.7554/eLife.23013 28282024
30 Mahnke YD et al The who’s who of T-cell differentiation: human memory T-cell subsets Eur J Immunol 2013 43 11 2797 2809 10.1002/eji.201343751 24258910
31 Evavold BD Allen PM Separation of IL-4 production from Th cell proliferation by an altered T cell receptor ligand Science 1991 252 5010 1308 1310 10.1126/science.1833816 1833816
32 den Braber I et al Maintenance of peripheral naive T cells is sustained by thymus output in mice but not humans Immunity 2012 36 2 288 297 10.1016/j.immuni.2012.02.006 22365666
33 Junge S et al Correlation between recent thymic emigrants and CD31+ (PECAM-1) CD4+ T cells in normal individuals during aging and in lymphopenic children Eur J Immunol 2007 37 11 3270 3280 10.1002/eji.200636976 17935071
34 Ryan DK et al Homeostasis of the naive CD4+ T cell compartment during aging J Immunol 2008 180 3 1499 1507 10.4049/jimmunol.180.3.1499 18209045
35 Kunzli M Masopust D CD4+ T cell memory Nat Immunol 2023 24 6 903 914 10.1038/s41590-023-01510-4 37156885
36 Thomson Z et al Trimodal single-cell profiling reveals a novel pediatric CD8αα+ T cell subset and broad age-related molecular reprogramming across the T cell compartment Nat Immunol 2023 24 11 1947 1959 10.1038/s41590-023-01641-8 37845489
37 Day CL et al Ex vivo analysis of human memory CD4 T cells specific for hepatitis C virus using MHC class II tetramers J Clin Invest 2003 112 6 831 842 10.1172/JCI200318509 12975468
38 Su LF et al Virus-specific CD4(+) memory-phenotype T cells are abundant in unexposed adults Immunity 2013 38 2 373 383 10.1016/j.immuni.2012.10.021 23395677
39 Bartolo L et al SARS-CoV-2-specific T cells in unexposed adults display broad trafficking potential and cross-react with commensal antigens Sci Immunol 2022 7 76 eabn3127 10.1126/sciimmunol.abn3127 35857619
40 Van Gassen S et al CytoNorm: a normalization algorithm for cytometry data Cytometry A 2020 97 3 268 278 10.1002/cyto.a.23904 31633883
41 Levine JH et al Data-driven phenotypic dissection of AML reveals progenitor-like cells that correlate with prognosis Cell 2015 162 1 184 197 10.1016/j.cell.2015.05.047 26095251
42 Becht E et al Dimensionality reduction for visualizing single-cell data using UMAP Nat Biotechnol 2019 37 38 44 10.1038/nbt.4314 30531897
43 Han A et al Linking T-cell receptor sequence to functional phenotype at the single-cell level Nat Biotechnol 2014 32 7 684 692 10.1038/nbt.2938 24952902
44 Masella AP et al PANDAseq: paired-end assembler for illumina sequences BMC Bioinformatics 2012 13 31 10.1186/1471-2105-13-31 22333067
45 Bolotin DA et al MiXCR: software for comprehensive adaptive immunity profiling Nat Methods 2015 12 5 380 381 10.1038/nmeth.3364 25924071
46 Gu Z et al circlize implements and enhances circular visualization in R Bioinformatics 2014 30 19 2811 2812 10.1093/bioinformatics/btu393 24930139
