key: cord-0986516-xngn52jh
authors: He, Chenfeng; Malone, Michael J.; Wendel, Ben S.; Ma, Ke-Yue; Del Alcazar, Daniel; Weiner, David B.; De Jager, Philip L.; Del Río-Estrada, Perla M.; Ablanedo-Terrazas, Yuria; Reyes-Terán, Gustavo; Su, Laura F.; Jiang, Ning
title: Transcriptome and TCR Repertoire Measurements of CXCR3(+) T Follicular Helper Cells Within HIV-Infected Human Lymph Nodes
date: 2022-05-06
journal: Front Immunol
DOI: 10.3389/fimmu.2022.859070
sha: 93f6a2be671146d0b13196a3b9f00690544b15a3
doc_id: 986516
cord_uid: xngn52jh

Follicular-helper T cells (T(FH)) are an essential arm of the adaptive immune system. Although T(FH) were first discovered through their ability to contribute to antibody affinity maturation through co-stimulatory interactions with B cells, new light has been shed on their ability to remain a complex and functionally plastic cell type. Due to a lack sample availability, however, many studies have been limited to characterizing T(FH) in mice or non-canonical tissue types, such as peripheral blood. Such constraints have resulted in a limited, and sometimes contradictory, understanding of this fundamental cell type. One subset of T(FH) receiving attention in chronic infection are CXCR3-expressing T(FH) cells (CXCR3(+)T(FH)) due to their abnormal accumulation in secondary lymphoid tissues. Their function and clonal relationship with other T(FH) subsets in lymphoid tissues during infection, however, remains largely unclear. We thus systematically investigated this and other subsets of T(FH) within untreated HIV-infected human lymph nodes using Mass CyTOF and a combination of RNA and TCR repertoire sequencing. We show an inflation of the CXCR3(+)T(FH) compartment during HIV infection that correlates with a lower HIV burden. Deeper analysis into this population revealed a functional shift of CXCR3(+)T(FH) away from germinal center T(FH) (GC-T(FH)), including the altered expression of several important transcription factors and cytokines. CXCR3(+)T(FH) also upregulated cell migration transcriptional programs and were clonally related to peripheral T(FH) populations. In combination, these data suggest that CXCR3(+)T(FH) have a greater tendency to enter circulation than their CXCR3(-) counterparts, potentially functioning through distinct modalities that may lead to enhanced defense.

Follicular-helper T cells (T FH ) provide B cells with the necessary costimulatory signals for their affinity maturation of broadly neutralizing antibodies (BNAbs) (1), a critical component of host defense. Despite high levels of circulating IgG in HIV-infected patients (2, 3) , BNAbs and natural viral control are uncommon (4, 5) . Given the importance of T FH contribution in a successful antibody response, it is possible their inadequacy during HIV infection leads to a lack of effective and timely B cell help.

Surprisingly, sheer numbers of T FH do not appear to be the problem. In fact, T FH exist at high frequencies during untreated HIV infection (6, 7) , suggesting a qualitative, rather than quantitative, impotence. In fact, we and others have shown that T FH in HIV + patients are functionally and phenotypically distinct compared to healthy donors (8, 9) , even during effective antiretroviral therapy (ART). In addition to their altered cell state, HIV-specific T FH appear to be slow-acting, requiring prolonged exposure to antigen likely during a period of uncontrolled viral replication (10) . Further complicating matters, T FH are readily infected by HIV. Several studies have implicated their role as an important reservoir for latent infection (11) (12) (13) (14) and subsequent viral persistence. This complex interplay between host and virus have made it a great challenge to understand the pitfalls of the immune system during HIV infection.

A high degree of cellular plasticity in the T FH compartment (15) (16) (17) has also led to our incomplete understanding of its role in HIV infection. T FH are generally identified by their expression of CXCR5, a chemokine receptor necessary for proper follicular localization in the lymph node (LN) (18, 19) . However, several less understood subsets have been characterized based on other cell surface markers. For example, in experimental autoimmune encephalomyelitis (EAE), a mouse disease model for multiple sclerosis, IL17a-expressing T FH have been found to infiltrate the central nervous system and exacerbate inflammation and disease (20) . In contrast, antigen-specific CCR6-expressing T FH were observed at high levels in SARS-COV-2 infection, correlating with less severe disease (21) . Interestingly, CCR6 + T FH in this context, normally described as TH17-like, seldom expressed IL17a. In an opposite functional context, regulatory T FH have also been described as FOXP3 + CXCR5 + T FH that appear to suppress the germinal center reaction (22) . Recently, we described an activated CXCR5 -T FH -like population in LNs that accumulates in HIV infection. These cells specifically coexpress high levels of ICOS and PD-1, are clonally related to T FH and can promote the production of HIV-specific antibodies from B cells in vitro (9) . We speculate that these cells, having a heightened migratory gene program, may infiltrate inflamed LNs and modify the B cell response; however, our understanding of their role in infection and relationship with other CXCR5 + T FH populations remains insufficient.

Another T FH subset that has received a fair amount of attention expresses the bona fide T H1 marker, CXCR3. A recent study depicting the clonal relationships between T FH subsets in the tonsils and blood of healthy humans found a strong clonal connection between CXCR3-expressing T FH (CXCR3 + T FH ) and circulating T FH (cT FH ) (23) . In agreement with previous publications (24, 25) , the authors posit that the origin of germinal center T FH (GC-T FH ) may be CXCR3 + T FH that initially encountered antigen in the periphery. Interestingly, other studies have shown CXCR3 + T FH to be enriched during chronic SIV infection in non-human primates (26, 27) , in the blood during HIV infection (28, 29) , and after influenza vaccination (30) . Given their inflation in non-steady state conditions, it is possible that studying CXCR3 + T FH in these contexts may reveal a clearer picture of their role leading up to and during their B cell interactions. Furthermore, with an altered cell state appearing to be a hallmark of T FH in HIV disease progression (8) , it is possible that understanding these alternative T FH subsets will bring light to novel mechanisms of disease and/or host defense. A systematic study of CXCR3 + T FH during perturbation, across multiple tissues in humans is thus of critical need to fully understand their function within the adaptive immune system.

Although blood studies in humans and work done with model organisms have laid the foundation for our understanding of T FH , their results have been contradictory at times. For example, Velu et al. observed similar helper function between LN CXCR3 + T FH and CXCR3 -T FH to induce SIVspecific IgG in vitro (27) , while other studies suggested that CXCR3 + T FH from human blood provide poor B cell help (31, 32) . It is likely that context, such as tissue type and disease state, has a great impact on the role and function of T FH . As such, our grasp remains limited without clear delineation of the relationship between cT FH and their native, lymphoidresident counterparts in both health and disease. In this study, we aim at narrowing this deficit in knowledge using a combination of Mass CyTOF and TCR repertoire and RNA sequencing on several different T cell populations from paired LNs and peripheral blood samples in untreated HIV + patients. Importantly, we find that CXCR3 + T FH are both inflated, phenotypically distinct, and correlate with lower HIV burden. We thus emphasize our analysis on this population and its relationship to several better characterized T cell populations within the human LN. Our multi-parameter approach revealed enhanced proliferative potential, an upregulation of cell migration pathways, and strong cross-tissue and crossphenotype clonal relationships within CXCR3 + T FH . These data suggest the potential of CXCR3 + T FH to expand within the LN, enter circulation and possibly contribute to host defense through alternative processes than their canonical GC-T FH counterparts. Further investigation into the impact of CXCR3 + T FH on viral load and spread over time with a larger, longitudinally tracked cohort will be important to unravel specific mechanisms that these may operate through.

Subjects included in this study are from a subset of patients recruited for our previous study (8) . Briefly, LN samples from HIV + donors were excised from palpable cervical LNs for clinical diagnostic workup in Mexico. HD samples were de-identified mesenteric or inguinal lymph nodes from the Cooperative Human Tissue Network (CHTN). Twenty-two HIV + donor LNs and nine HD LNs were used for Mass CyTOF experiments. Considering sample availability and size, LN and PBMC samples from an independent, but semioverlapping cohort (Tables S1, S5) were collected from LN and PBMC samples were collected from seven non-ART and one ART-treated HIV + patients. Sample sizes were not prespecified but dictated by the availability of the sample. All samples from HIV+ patients were de-identified and were obtained in accordance with the Declaration of Helsinki after obtaining written informed consent of participants and as part of protocol B03-16, which was reviewed and approved by the Research Committee and the Ethics in Research Committee of the National Institute of Respiratory Diseases "Ismael Cosıó Villegas", Mexico City.

Mass CyTOF data were collected as described in a previous publication (8) , except that the cells were not stimulated prior to staining. Data were collected and normalized as described there. In analysis specific to this paper, cells were subset by drilling down on the gates defined in Figure S1 with FlowJo ™ v10.8.0 Software (BD Life Sciences) and exported as individual FCS files. FCS files were then read into R using the "flowCore" package and analyzed using custom scripts. Dimensional reduction and clustering were done using the R package "Seurat". Cluster annotations and rational are described in Table S2 . Boxplots and UMAP plots were generating using the R package "ggplot2". The heatmap in Figure 1 was made with the R package "pheatmap" (33) . Statistical analyses were done using the R package "rstatix". Final plots were compiled in Biorender.

Cell staining and sorting were performed as previously described (9) . Briefly, cryopreserved LN cell suspensions were freshly thawed. A BD FACSARIA was used to sort live CD4 + T cells as CD3 + CD4 + CD8 -CD14 -CD19 -. Then, from the live CD4 + T cell population, naïve (CD45RO -CXCR5 -CD57 -CCR7 + ), CXCR3 + T FH (CD45RO + CXCR5 + CCR7 -PD1 + CXCR3 + ), CXCR3 -T F H ( C D 4 5 R O + C X C R 5 + C C R 7 -P D 1 + C X C R 3 -) , G C -T F H (CD45RO + CXCR5 + CD57 + ) and CXCR5 -T FH (CD45RO + CXCR5 -PD1 + ICOS + ) populations were sorted directly into lysis buffer and stored at -80°C until further use. 

Total RNA was purified using AllPrep DNA/RNA Mini and RNeasy MiniElute kits (Qiagen) according to the manufacturer's instructions. Purified RNA quality was evaluated using Agilent RNA 6000 Pico Kit. Samples were then either supplemented with RNase inhibitor (RNase OUT, Thermofisher Scientific) and stored at -80°C or taken directly to reverse transcription.

TCRb repertoire sequencing libraries were prepared and sequenced as described previously (8) . Briefly, 30% of purified RNA was used for reverse transcription. Molecular identifiers (MIDs) were added to cDNA templates during reverse transcription using a TRBCbinding primer with 12 random nucleotides and a partial Illumina adapter (Table S4) . PCR1 was carried out using multiplexed TRBVbinding primers. A second round of PCR was then used to append final Illumina sequencing adapters to the TCRb junctioncontaining inserts. Final libraries were sequenced on an Illumina paired end 150x150 configuration at a minimum depth of 10 reads/ cell. After sequencing, reads were clustered based on their MID and TCRb sequence similarity. Consensus sequences were then built to correct PCR and sequencing errors as described (34) . TCRb sequencing information is summarized in Table S5 . The CDR3 blast module MIGEC (35) was used for TRBV/J and CDR3 annotation. Bhattacharyya coefficient (36) was used to measure TCR repertoire similarity. Circos plots were generated using circlize R package (37) .

RNA-Seq libraries were prepared using a protocol modified from SMART-seq2 (SSII) (38) . Briefly, 2ul of purified total RNA was reverse transcribed using a poly-T primer fused to the ISPCR handle (Table S6 ; RT_dT30VN). Second-strand synthesis was then done using the SSII template-switching oligo (Table S6 ; SSII_TSO). The following program was used: 42°C RT for 90min, 10 cycles of 50°C for 2min, 42°C for 2min and then 70°C for 15min. cDNA amplification was then done using KAPA HIFI and the ISPCR handle (Table S6 ; RT_TSPCR) with the following program: 98°C initial denature for 3min, 16 cycles of 98°C denature for 20sec, 64°C annealing for 30sec, 72°C extension for 6min and 72°C final extension for 5min. The PCR product was then purified using AmpureXP beads (Agencourt), according to the manufacturer's protocol. Purified PCR product was then diluted to 0.1 -0.3 ng/ul and tagmented using the Nextera XT kit (Illumina) with a reduced reaction volume. Briefly, 1.25 ul of diluted sample was used in a 5ul total reaction volume and fragmented for 10min at 55°C, and then held at 10°C. The reaction was then neutralized by adding 1.25 ul of NT buffer. Final libraries were then generated from the tagmentation product using Nextera adapters (Table S6 ) and the following program: 72°C for 3 min, 95°C denature for 3sec, 12 cycles of 95°C denature for 10sec, 55°C annealing for 30sec, 72°C extension for 1min and 72°C final extension for 5min. Indexed libraries were then purified using Ampure XP beads according to manufacturer's instructions and quantified with an Agilent High Sensitivity DNA kit. The final libraries were pooled and sequenced on Illumina HiSeq (150x150) with a minimum depth of about 200 reads per cell. Each library was split into two technique replicates and sequenced independently to ensure the reliability of RNA sequencing results. Sequencing reads were aligned to human reference genome GRCh38.p5 using RSEM (39) . Differentially expressed genes were analyzed using DESeq2 (40) . GSVA, a non-parametric unsupervised method that quantifies the relative enrichment of selected pathways, was performed using the R package, GSVA (41) . Gene Set Enrichment Analysis (GSEA) on selected pathways or gene sets were performed with the R package fgsea (42) . Tonsil Tfh signatures (GSE50391, CXCR5 high CD45RO + versus CXCR5tonsil samples) were created using the GEO2R online tool (https://www.ncbi.nlm. nih.gov/geo/geo2r/). Metascape (43) and NetworkAnalyst (44) was used to determine GO pathway enrichment.

HIV transcripts were quantified from bulk RNA-Seq following the stepwise procedures ( Figure S19 ) similar to the method outlined in Chang et al. (45) . Briefly, reads that could not be aligned to the human genome were mapped to the HIV genome (NC_001802.1) with the gapped aligner software HISAT2 (46), accounting for splice junction events. Quantification of HIV transcript expression was calculated by normalizing the HIV genome mappable reads to human genome mappable reads and the input cell number. For a more reliable quantification of HIV copies in each library, the correlation between two replicates was first evaluated ( Figure S20 ), then the HIV copy numbers of both technical replicates were averaged.

A paired sample two-tailed Student's t-test was used for pairwise comparisons, with a P value less than 0.05 as a cutoff to determine statistical significance. These comparisons include CXCR3 + /CXCR3 -T FH percentages in Figure S10C , TCR repertoire similarities in Figures 4B, C and HIV copies in Figure 5A . For gene set enrichment analysis, CAMERA (47) (Correlation Adjusted MEan RAnk gene set test) was used to calculate enrichment significance. Pearson correlation was used to determine the degree of association.

We first set out to determine if, in concordance with SIV infection, CXCR3 + T FH were inflated in human LN during HIV infection. We used a 29 marker Mass CyTOF panel to survey T cells from 22 HIV + and 9 healthy donor LNs (Characteristics of the donors are listed in Table S1 ). We compared the frequencies of four T FH populations using the gating strategies described in Figure S1 . Two of these populations, denoted as MEM T FH and CXCR3 -T FH , are defined as CD45RO + CXCR5 + and CD45RO + CXCR5 + PD1 + CXCR3 -CD4 + T cells, respectively. These two gating schemes are frequently used when characterizing T FH (8, 32, 48) . The remaining two populations, C X C R 3 + T F H a n d C X C R 5 -T F H , are defi n e d a s CD45RO + CXCR5 + PD1 + CXCR3 + and CD45RO + CXCR5 -PD1 + ICOS + CD4 + T cells, respectively. CXCR3 + T FH have been characterized in rhesus macaques (RMs) and human blood (27, 31, 32) , while CXCR5 -T FH were recently described as functional T FH in human LN and blood during HIV infection (9) . Amongst the four populations analyzed, CXCR3 + T FH and CXCR5 -T FH were significantly enriched in HIV + patient LNs versus healthy donor LNs ( Figure 1A) . It should be noted, however, that all LNs from HIV + donors were cervical, whereas HD LNs were derived from a mix of distinct anatomical locations (Table S1) . Regardless, given their high frequency over steady state, we suspected these two T FH subsets may play a role in HIV infection. To directly address this question, we compared the frequencies of T FH populations with each patient's plasma viral load (pVL) ( Figure S2) . Interestingly, the only metric correlated to pVL is the ratio of CXCR3 + T FH to CXCR3 -T FH ( Figure 1B) , suggesting that donors with greater CXCR3 + T FH populations may be better at controlling virus.

To gain insight into the role and interplay of these T FH populations in HIV infection, we sought to delve deeper into the Mass CyTOF data. To do so, we subset the data by drilling down on four T cell populations (CXCR3 + T FH , CXCR3 -T FH , CXCR5 -T FH and Naïve) described in Figure S1 , effectively sorting these populations in silico. We then projected the in silico sorted populations onto two dimensions using Uniform Manifold and Projection (UMAP). At a first glance, cells derived from HIV + donors occupied several distinct locations on the UMAP ( Figure S3) . Additionally, although modest heterogeneity was observed between the three in silico sorted T FH populations, each tended to occupy different regions at different densities ( Figure S4 ). We thus decided to use unsupervised clustering to gain unbiased insight into each of the populations and their relationships with one another. Twelve clusters emerged, several representing canonical CD4 + T cell and T FH phenotypes ( Figure 1C ). Select signature markers ( Figure S5 ), as well as the heatmap in Figure 1C were used to annotate clusters 0-8. Clusters 9-11 were not considered in downstream analyses due to their insignificant sizes. To evaluate the impact of ART on phenotype distribution, we next calculated the frequencies of ART + and ARTpatients in each cluster ( Figure S6 ). Although some clusters show some trending, we found no significant phenotypic skewing and thus proceeded to group both ART + and ARTpatients as HIV + in subsequent analyses of this dataset.

To determine unique features of HIV infection, we compared the distribution of T cells in HIV + and healthy LN-derived cells in each cluster ( Figure 1D ). In agreement the elevated expansion of the T FH compartment seen in the literature (6, 7), we found that cluster 5, denoted as proliferating T FH based on high expression of Ki67 and BLyS, was significantly enriched in HIV + donor LN cells (p < 0.001). Surprisingly, even though cluster 2 had the highest and most homogeneous expression of CXCR3 ( Figure S5 ), we found no difference between HIV + and healthy cells. This might reflect that many of the CXCR3 + T FH in HIV infection take on diverse phenotypic programs. Specifically, given the proximity of cluster 2 with the NAÏVE cluster, it is possible that the CXCR3 + T FH within cluster 2 represent early T FH entering the LN as previously described (20) , whereas CXCR3 + T FH taking on alternative phenotypic programs may represent other stages of the T FH lifecycle. In line with this argument, we reasoned that cluster 5, in addition to being enriched in HIV + donors, might also have signatures unique to HIV infection. To test this hypothesis, we subsampled cluster 5 to have equal numbers of HIV + and healthy donor cells and ran a likelihood-ratio test on each marker (49) . We found that CXCR3 expression in cluster 5 was in fact significantly higher in HIV + than healthy donors ( Figure 1E) . Given the HIV-intrinsic upregulation of CXCR3 in cluster 5, we reasoned that the in silico sorted populations might have differential distributions within cluster 5. As expected, CXCR3 + T FH were significantly enriched in cluster 5 over all other populations in HIV + patients, but not healthy donors ( Figure S7) . Concordantly, CXCR3 + T FH from HIV + donors occupied cluster 5 at a significantly higher frequency than from healthy donors, whose predominant population was contrastingly CXCR3 -T FH ( Figure S8 ). In addition to CXCR3, Granzyme A, a marker indicative of cytotoxic function in T-helper cells (50) , was also enriched in cells from HIV + donors within cluster 5 ( Figure 1E ). It is possible that, although CXCR3 + T FH generally reside in a more quiescent, immature state in the steady state, HIV infection triggers them to activate, proliferate in the LN, and take on phenotypic signatures unique to HIV infection.

To better understand the relationship between the phenotypic signatures that the in silico sorted T FH populations ( Figure S1 ) take on specifically in HIV infection, we did pairwise comparisons of the distributions of CXCR3 + , CXCR3 -, and CXCR5 -T FH populations within HIV+ donors across each cluster ( Figure 1F ). Reflecting the cellular plasticity of T FH , gating down on only a handful of markers unsurprisingly revealed modest heterogeneity within each population. As expected, however, cluster 2, defined by the highest expression of CXCR3, bore a significantly large proportion of CXCR3 + T FH compared to all other populations. In line with previous studies (8, 9) , CXCR5 -T FH were significantly skewed toward activated states in clusters 5 and 6 (Proliferating T FH and EOMES + T FH , respectively), while CXCR3 + and CXCR3 -T FH both tended to exist in cluster 7 (Activated GC-T FH ). Although statistical significance was only reached between CXCR3 + T FH and CXCR3 -T FH in cluster 2, we noticed that CXCR3 + T FH often existed at a frequency between CXCR3 -T FH and CXCR5 -T FH , suggesting this population might exist as a transitional state between the canonical T FH and the recently characterized CXCR5 -T FH . Given this possibility, we next evaluated the median position of each HIV+ donor's populations on the same UMAP coordinates as in Figure 1C ( Figure 1G ). At first glance, the Naïve population is clearly distinct from the T FH populations, which group together on the right side of the UMAP. We also noted that, as suggested in Figure 1F , CXCR3 + T FH appear at an intermediate position between CXCR3 -T FH and CXCR5 -T FH . Statistical analysis of each population also revealed that all populations, except CXCR3 + T FH and CXCR5 -T FH , were distinct over UMAP_1 (Table S3) . Additionally, hierarchical clustering of each population from each donor revealed a grouping of CXCR3 + T FH with CXCR5 -T FH more than with CXCR3 -T FH ( Figure S9 ). In concordance with our hypothesis, these data also posit CXCR3 + T FH as an interim population between CXCR3 -T FH and CXCR5 -T FH .

Our findings in the Mass CyTOF dataset prompted us to gain deeper insight into the similarities and differences of these T FH populations in HIV infection. Although Mass CyTOF can be useful for delineating cell states due to its ability to analyze large numbers of single cells rapidly and robustly for a selected panel of surface markers, it failed to provide an unbiased analysis of the transcriptome differences between samples. Additionally, due to the destructive nature of Mass CyTOF, TCR sequences, and thus clonal relationships, cannot be obtained. To circumvent these limitations, we sorted five populations of T cells for bulk RNA sequencing (RNA-Seq) from 7 HIV + , non-ART-treated donors using the gating strategy illustrated in Figure S10 . We excluded the one ART-treated donor in this section as we were unsure how treatment would affect a highly sensitive assay such as RNA-seq. We focused specifically on five populations. The first two, CXCR3 + T FH (CD45RO + CXCR5 + CCR7 -PD1 + CXCR3 + ), and CXCR3 -T FH (CD45RO + CXCR5 + CCR7 -PD1 + CXCR3 -), were sorted to delineate the specific effects of CXCR3 expression on the state of LN T FH . GC-T FH (CD45RO + CXCR5 + PD1 + CD57 + ) and CXCR5 -T FH (CD45RO + CXCR5 -ICOS + PD1 + ) were sorted as two distinct reference points between two subsets that have been demonstrated to provide B cell help (8, 9) . CD57-expressing T FH were sorted as a proxy for GC-T FH as they have been shown to be a major subset of active GC cells (51) (52) (53) , which serves as a positive control for T FH function. An initial analysis on CyTOF and Flow Cytometry data shows comparable expression of CD57 between CXCR3 + T FH and CXCR3 -T FH populations ( Figures S11, S12) . Naïve CD4+ T cells CD45RO -CXCR5 -CCR7 + ) were sorted as a non-T FH control.

Given a clear enrichment of CXCR3 + T FH within HIV + patient LNs in the Mass CyTOF dataset, we wanted to understand their transcriptional characteristics in HIV infection. Of particular concern, existing studies with animal LN (27) and human blood (31) show contradicting results on their potential for B cell help. We thus set out to explore the transcriptional landscape of CXCR3 + T FH in relation to the other, better characterized, T FH populations described above using RNA sequencing (RNA-Seq). We projected the RNA-Seq data from each population onto 2dimensional space using Principal Component Analysis (PCA). As expected, the naïve population was distinct from all the other populations based on PC1. The T FH populations, although similar on PC1, were separated based on PC2 (Figure 2A ). In concordance with our analysis on the Mass CyTOF data, CXCR5 -T FH were distinguishable from GC-T FH and CXCR3 -T FH on PC2, while CXCR3 + T FH situated in the middle. These data again suggest CXCR3 expression on T FH may be indicative of a transitional state between the canonical GC-T FH and the recently described motile CXCR5 -T FH (9) . Further advocating for CXCR3 + T FH as a transitional state, we also noticed CXCR3 expression in both flow cytometry and Mass CyTOF was higher in CXCR5 -T FH than CXCR5 + T FH ( Figure S13 ). To better understand the transcriptome differences among various T FH populations, we next used Metascape (43) to evaluate the pathways enriched in genes contributing to PC2 ( Figure S14) . Interestingly, genes that mostly negatively correlated with PC2 (upregulated in CXCR5 -T FH while downregulated in GC-T FH ) were enriched in cell division and migration. However, genes that most positively correlated with PC2 (upregulated in GC-T FH while downregulated in CXCR5 -T FH ) were enriched on pathways inhibiting cell proliferation and migration (Table S7) . For example, CTLA4 transmits inhibitory signal to T cells proliferation (54) , SMAD3 mediates inhibition of CD4 T-cell proliferation (55) , MCC blocks cell cycle progression from G0/ G1 to S phase (56) . In support of our CyTOF dataset, CXCR3 + T FH and CXCR5 -T FH appear to be in a more activated, motile state than their canonical GC-T FH counterpart.

To evaluate the differences between CXCR3 + and CXCR3 -T FH more directly, we next analyzed differentially expressed genes (DEGs, Figure 2B ) between the two populations. In total, 86 out of 12,515 genes were identified as significant DEGs (Table S8) . Hierarchical clustering based on those DEGs further demonstrated the similarity of CXCR3 -T FH with GC-T FH and CXCR3 + T FH with CXCR5 -T FH ( Figure 2C) . A significant upregulation of canonical T H1 genes (CXCR3, IFNG, CXCR6 and GBP1) in CXCR3 + T FH also indicated their unique T H1 -like program. Additionally, upregulation of NELL2 and MKI67 in CXCR3 + T FH suggested an increased proliferative capacity compared to their CXCR3counterpart. Interestingly, CXCR3 + T FH also upregulated ITGB7, encoding a subunit of the integrin a4b7, that plays a role in leukocyte adhesion and can serve as a homing receptor bound by HIV (27) . Thus, it is possible that CXCR3 + T FH could be similarly, or even more, susceptible to HIV infection than other T FH subsets. Surprisingly, IL10, encoding an unconventional cytokine (IL10) in T H1 cells, was also highly expressed by CXCR3 + T FH . Since previous studies support IL10 as a key player in the establishment and perpetuation of HIV persistence (57) , the upregulation of IL10 in CXCR3 + T FH may result in enhanced persistence of HIV after infection.

To infer the functional potential of human LN-derived CXCR3 + T FH , we next compared them to CXCR3 -T FH specifically on T FH -related genes. We created a T FH signature gene set using previously published RNA-Seq data from human tonsil samples (31) . Both CXCR3 + T FH , CXCR3 -T FH and GC-T FH appeared similar when hierarchically clustered on tonsil T FH signature genes ( Figure S15A ). GSEA and GSVA based on the same gene set also revealed no significant differences between these two cell populations (Figures S15B, C) . Detailed analysis on several manually curated T FH -related genes also suggested CXCR3 + and CXCR3 -T FH bear similar T FH marker gene expression patterns ( Figure S15D ). For example: FOXO1, a negative regulator of BCL6 (58) , was downregulated in both populations. Additionally, BCL6, MAF and CD84, key transcriptional regulators of T FH differentiation, were similar in both CXCR3 + and CXCR3 -T FH . B3GAT1 (an enzyme necessary for the production of CD57), which is specifically expressed by active GC-T FH (8, (51) (52) (53) , was also comparably expressed in both populations. As expected, CXCR5 and its upstream regulator ASCL were both similarly expressed in CXCR3 + and CXCR3 -T FH . Taken together, these observations reveal that CXCR3 + T FH have a similar T FH transcriptional program to CXCR3 -T FH , suggesting a functional overlap between the two populations.

In summary, global transcriptome analysis depicts CXCR3 + T FH unique from GC-T FH , biasing toward a T H1 -like program. It also revealed CXCR3 + T FH are similar to the recently described CXCR5 -T FH . Focusing specifically on T FH marker genes, significant differences between CXCR3 + and CXCR3 -T FH were not observed, supporting the paradigm that human LN-derived CXCR3-expressing T FH may still bear T FH function.

To investigate deeper into the functional program of CXCR3 + T FH , we used NetworkAnalyst (44) to identify GO biological pathways enriched within this population. Only two GO pathways were significantly enriched within CXCR3 + T FH when compared with CXCR3 -T FH . Interestingly, both pathways were related to T cell migration ( Figure 3A) . We visualized the differentially expression genes of these two pathways in detail, as expected, most of the genes were upregulated in CXCR3 + T FH ( Figure 3B; Figure S16 ), which suggested a possibility of this cell subset to be more motile. Our previous study (9) suggested that T FH can downregulate CXCR5 expression and accumulate as CXCR5 -T FH in LNs during HIV infection. The same study also found that these CXCR5 -T FH provide B cell help and have a propensity to migrate into the periphery. The combination of CXCR3 + T FH also bearing a cell migration signature, as well as existing in an intermediate cell state between CXCR5 -T FH and GC-T FH based on both Mass CyTOF and RNA-Seq data, led us to hypothesize that CXCR3 + T FH may be the immediate relative of CXCR5 -T FH . To test this hypothesis more directly, we compared the transcriptomes of CXCR3 + and CXCR3 -T FH using GC-T FH and CXCR5 -T FH as reference gene sets. Specifically, we first identified 938 DEGs from 12,559 total A B D C FIGURE 3 | CXCR3 + T FH uniquely express cell migratory genes and transcriptomic-ally similar with CXCR5 -T FH . DEGs between CXCR3 + T FH and CXCR3 -T FH were enriched on T cell migration pathways (A) with migratory genes upregulated in CXCR3 + T FH (Leukocyte chemotaxis, GO:0030595) (B) GSEA analysis indicates genes elevated in CXCR5 -T FH (Gene set 1) were upregulated in CXCR3 + T FH (C) while genes upregulated in GC-T FH (Gene set 2) were also elevated in CXCR3 -T FH (D) GSEA were performed with R package 'fgsea', specifically, genes were pre-ranked by their p-values, which indicates whether one gene is highly expressed in CXCR3 + T FH or CXCR3 -T FH . The bar on the X-axis indicates one gene from a selected gene set (i.e., either Gene set 1 or Gene set 2). The curve is a running sum of the bars on the X-axis. shared genes between GC-T FH and CXCR5 -T FH . Among these DEGs, 514 were upregulated in CXCR5 -T FH (Gene set 1) and 424 were upregulated in GC-T FH (Gene set 2). We then performed GSEA analysis to compare CXCR3 + and CXCR3 -T FH on Gene Set 1 and 2 ( Figures 3C, D) . As hypothesized, CXCR3 + T FH followed the CXCR5 -T FH program, while CXCR3 -T FH more closely followed the GC-T FH program. Together with our observations in Mass CyTOF, global transcriptome, differential gene expression, gene set enrichment, and pathway analyses, we conclude that CXCR3 + T FH likely exist in T Cell Receptor Repertoire Demonstrates That CXCR3 + T FH Are Clonally Related to CXCR5 -T FH and Peripheral T FH After observing that CXCR3 + T FH upregulated cell migratory pathways and skew phenotypically toward CXCR5 -T FH , we suspected that changes in CXCR3 expression on T FH might facilitate transitions to and from GC-T FH and CXCR5 -T FH cell states. To test this hypothesis, we used bulk T cell receptor (TCR) repertoire sequencing to measure their in vivo clonal relationship with various T cell populations. Because of the immense diversity generated by V(D)J recombination, TCR sequences can be thought as unique 'ID-Cards' (59) , where cells sharing a given TCR sequence are likely to be the progenies of the same ancestor. By comparing the overlap of TCR sequences among samples, a clonal lineage across tissues and cell states can be inferred. We thus compared the TCR repertoire similarity across the five LN T cell populations using Bhattacharyya Coefficient (36) as the similarity index. Bhattacharyya coefficient (BHA) measures the overlap of clonotypes between two T cell populations, while taking size of clones into consideration. In contrast to previous findings in healthy human tonsils (23), CXCR3 + T FH , CXCR3 -T FH and GC-T FH were strikingly similar in their repertoires across all donors ( Figures 4A, C; Figure S17 ). This suggests that these three populations are likely capable of transitioning in and out of each of these states with ease, which may be a phenomenon intrinsic to HIV infection. In concordance with transcriptomic data, however, CXCR3 + T FH were significantly more related to CXCR5 -T FH than CXCR3 -T FH across all donors ( Figures 4A, B) . For two out of the seven donors, we were able to compare the CXCR3 + and CXCR3 -T FH repertoires with the cT FH repertoires in the blood (60) (Figures 4D-F; Figure S18 ). Circos plots representing the overlap of TCRs between cT FH and LN populations indicate that CXCR3 + T FH are also more clonally related to cT FH than their CXCR3 -T FH counterpart ( Figure 4E ), similar to a recent finding in healthy tonsils (23) . A similar trend was observed considering clone size weighted BHA index ( Figure 4F ). In combination, high clonal overlap amongst all CXCR5-expressing T FH in the LN, and the elevated relationship of CXCR3 + T FH with CXCR5 -T FH and cT FH , suggests that this population sits as a bridge between secondary lymphoid tissue and the periphery.

CXCR3 + T FH in the Human LN Are a Reservoir for HIV T FH are major reservoirs for HIV (61) , complicating their role in controlling virus through the coordination of GC reactions. Furthermore, CXCR3 + T FH have been speculated to maintain a dynamic HIV reservoir due to high expression of HIV co-receptors (i.e. CCR5 and a4b7) (62) . High viral load in this highly motile cell type could either exacerbate systemic spread of virus or expose viral antigens to the immune system. Given our result that CXCR3 + T FH frequency negatively correlates with pVL ( Figure 1C ), we speculated a high level of HIV copies in CXCR3 + T FH might be a mechanism for improving clearance of virus.

To investigate this idea, we mapped the bulk RNA-Seq reads from each T cell population to the HIV genome to quantify their relative HIV infection intensities (see methods for details). We detected comparable levels of HIV transcripts in both CXCR3 + and CXCR3 -T FH , although both were significantly higher than Naïve cells ( Figure 5A ; Figure S21 ). Similar to experiments in RMs (27), HIV abundance also negatively correlated with repertoire diversity ( Figure 5B) , reflecting the propensity of proliferating and antigenexperienced cells to bear a higher HIV burden. Given that the two T FH subsets discussed here share a similarly high HIV burden, the consequence of the migratory potential of CXCR3 + T FH is likely a complex balance between limiting systemic spread of virus and exposing it to novel immune compartments.

T FH are paramount to the elicitation of BNAbs. Whether induced by vaccination or natural infection, BNAbs targeting invading pathogens have been shown in a variety of contexts to mitigate, eliminate, and even prevent disease (63, 64) . Some viruses, however, can escape from or disrupt this response, often leading to uncontrolled viral replication and chronic infection. In HIV infection, natural BNAbs are seldom produced and almost never lead to viral control (4, 5) . Given that T FH are a major site of HIV infection (11) (12) (13) , it is possible that alternative phenotypic and functional programs within the T FH compartment are taken on to overcome this deficit in humoral immunity. Unfortunately, sampling T FH from humans in their native lymphoid tissue in untreated HIV infection poses an immense challenge in the field. Here, we have overcome this obstacle, combining several high-throughput assays and analyses on T FH from HIV + patient LNs to gain deeper insight into their role in disease.

Mounting evidence points to a high degree of cellular plasticity within the T FH compartment (15) (16) (17) . We show here the inflation of CXCR3 + T FH in HIV + patient derived LNs compared to healthy donors that correlates with a lower viral load. Further investigation into this population revealed its abnormal cell state specifically within HIV+ donors that positions it in an intermediate phenotypic state between GC-T FH (8, 31) and the recently described CXCR5 -T FH (9) . Deeper analysis into the clonal relationship of CXCR3 + T FH corroborated these results, revealing its heightened repertoire similarity with both CXCR5 -T FH and cT FH . Additionally, transcriptome analysis revealed a propensity for CXCR3 + T FH to upregulate cell migration pathways. Taken together, we posit CXCR3 + T FH as a bridge between lymphoid tissue and the periphery. Given that these CXCR3 + T FH bear a high HIV burden, may be primed for cellular migration, and are affiliated with a lower viral load, it is possible that their movement in and out of lymphoid tissues leads to better viral control by exposing antigens to novel immune compartments. If so, it is possible that the upregulation of CXCR3 in T FH may be an effort of the immune system to drive viral reservoirs out of hiding. Given the propensity of low levels of HIV to remain through ART (12) and the increased likelihood that BNAbs are produced during high antigen exposure (10), a mass exodus of T FH may be beneficial to the host during early infection or ART. Furthermore, targeting these migration pathways has the potential to be a viable therapeutic. Of course, it is also possible CXCR3 upregulation in T FH is not HIV-intrinsic, but rather a product of chronic infection. In fact, previous research has implicated CXCR3 + T FH in HCV infection (60) , Zika virus infection (65) , and acute febrile malaria (66) and vaccination (67) in children. Interestingly, however, while CXCR3 + T FH are affiliated with neutralizing antibodies in HCV and Zika virus infection, they appear to lead to poor prognosis in malaria. Understanding the differences and similarities between CXCR3 + T FH in each of these disease contexts, specifically in lymphoid tissues as addressed in this study, will shed light on different immune-intrinsic properties of chronic infection.

An important limitation of this study was our inability to directly measure T FH function on the populations of interest. Although we were able to point to specific genes, pathways, and signatures that suggest their ability to provide B cell help, definitive knowledge of the functional capacity of CXCR3 + T FH in human LNs will be important to understanding their role in HIV infection. Future studies should aim at accomplishing this. Additionally, the spectral limitations of FACS also prevented us from being able to sort more specific T FH subsets that may be important to resolve the unique phenotypic programs of CXCR3 + T FH in HIV infection. For example, our analyses on the Mass CyTOF data suggested that CXCR3 + T FH take on a unique proliferative program specific to HIV + patients that was not robustly seen in healthy donors, a phenomenon that may explain some of the repertoire discrepancies between this study and others (20) . However, since we were unable to accommodate more markers when sorting for TCR and RNA sequencing, especially with the limited number of cells within each LN sample, these unique HIV-intrinsic factors may have been diluted. Recent innovations in antibody barcoding and single cell RNA-seq coupled with creative approaches to isolate these unique T FH populations may lead to a clearer understanding of their phenotypic and clonal relationships.

In summary, we have evaluated several LN-blood-matched T FH from a small, albeit rare cohort of untreated HIV + patients using Mass CyTOF and RNA and bulk TCR repertoire sequencing. Our analyses revealed a phenotypic shift of CXCR3 + T FH from a GC-T FH cell state toward an unconventional CXCR5 -T FH state. CXCR3 + T FH also upregulated migratory transcriptional programs and were clonally related to peripheral T FH cell populations. Altogether, these data suggest that CXCR3 + T FH may be transitional state between their CXCR3lymphoid and peripheral counterparts. Future work aimed at delineating the temporal relationships of these T FH populations, drilling deeper into the function of more specific phenotypic niches, will be pertinent to fully understanding their implication in HIV infection.

Publicly available datasets along with new sequencing data were analyzed and can be found here: https://www.ncbi.nlm.nih.gov/ projects/gap/cgi-bin/study.cgi?study_id=phs001548.v2.p1.

Subjects included in this study are from a subset of patients recruited for our previous study (8) . All samples were deidentified and obtained with IRB regulatory approval from the University of Pennsylvania. The patients/participants provided their written informed consent to participate in this study.

CH and MJM analyzed data, performed research, and wrote the manuscript; MJM, BSW and K-YM performed sequencing experiments; DDA performed cell sorting experiments; DBW and PLDJ helped with data interpretation; PMDR-E, YA-T, and GR-T established the infrastructure to recruit HIV + patients, provided HIVinfected samples, and associated clinical information. NJ and LFS designed the study; NJ directed the study. All authors contributed to the article and approved the submitted version. 

Tfh Cells and HIV Bnabs, an Immunodominance Model of the HIV Neutralizing Antibody Generation Problem

Relationship of Serum Immunoglobulin and IgG Subclass Levels to Race, Ethnicity and Behavioral Characteristics in HIV Infection

Serum Immunoglobulin G Subclass Dysbalances in the Lymphadenopathy Syndrome and Acquired Immune Deficiency Syndrome

The Characterization of Non-Progressors: Long-Term HIV-1 Infection With Stable CD4+ T-Cell Levels

Dissecting the Neutralizing Antibody Specificities of Broadly Neutralizing Sera From Human Immunodeficiency Virus Type 1-Infected Donors

Spatial Alterations Between CD4 + T Follicular Helper, B, and CD8 + T Cells During Simian Immunodeficiency Virus Infection: T/B Cell Homeostasis, Activation, and Potential Mechanism for Viral Escape

Expansion of HIV-Specific T Follicular Helper Cells in Chronic HIV Infection

The Receptor Repertoire and Functional Profile of Follicular T Cells in HIV-Infected Lymph Nodes

Mapping the Lineage Relationship Between CXCR5+ and CXCR5-CD4+ T Cells in HIV-Infected Human Lymph Nodes

Quality and Quantity of TFH Cells Are Critical for Broad Antibody Development in SHIVAD8 Infection

Multiple Origins of Virus Persistence During Natural Control of HIV Infection

Peripheral T Follicular Helper Cells Are the Major HIV Reservoir Within Central Memory CD4 T Cells in Peripheral Blood From Chronically HIV-Infected Individuals on Combination Antiretroviral Therapy

Germinal Centre CD4+ T Cells Are an Important Site of HIV Replication In Vivo

Follicular Helper T Cells Serve as the Major CD4 T Cell Compartment for HIV-1 Infection, Replication, and Production

Heterogeneity of Phenotype and Function Reflects the Multistage Development of T Follicular Helper Cells

Helper Cell Heterogeneity: Time, Space, and Function

Mapping the Diversity of Follicular Helper T Cells in Human Blood and Tonsils Using High-Dimensional Mass Cytometry Analysis

Helper Cell Differentiation, Function, and Roles in Disease

Helper Cell Biology: A Decade of Discovery and Diseases

Role of TFH Cells in Promoting T Helper 17-Induced Neuroinflammation

Antigen-Specific Adaptive Immunity to SARS-CoV-2 in Acute COVID-19 and Associations With Age and Disease Severity

Follicular Regulatory T Cells Can Access the Germinal Center Independently of CXCR5

CD4+ T Follicular Helper Cells in Human Tonsils and Blood Are Clonally Convergent But Divergent From Non-Tfh CD4+ Cells

Circulating Precursor CCR7loPD-1hi CXCR5+ CD4+ T Cells Indicate Tfh Cell Activity and Promote Antibody Responses Upon Antigen Reexposure

Do Memory CD4 T Cells Keep Their Cell-Type Programming: Plasticity Versus Fate Commitment?: Complexities of Interpretation Due to the Heterogeneity of Memory CD4 T Cells, Including T Follicular Helper Cells

Severe Malaria Infections Impair Germinal Center Responses by Inhibiting T Follicular Helper Cell Differentiation

Induction of Th1-Biased T Follicular Helper (Tfh) Cells in Lymphoid Tissues During Chronic Simian Immunodeficiency Virus Infection Defines Functionally Distinct Germinal Center Tfh Cells

Frequencies of Circulating Th1-Biased T Follicular Helper Cells in Acute HIV-1 Infection Correlate With the Development of HIV-Specific Antibody Responses and Lower Set Point Viral Load

Circulating CXCR5+CXCR3+PD-1lo Tfh-Like Cells in HIV-1 Controllers With Neutralizing Antibody Breadth

Induction of ICOS+CXCR3+CXCR5+ T H Cells Correlates With Antibody Responses to Influenza Vaccination

Human Circulating PD-1+CXCR3-CXCR5+ Memory Tfh Cells Are Highly Functional and Correlate With Broadly Neutralizing HIV Antibody Responses

Human Blood CXCR5+CD4+ T Cells Are Counterparts of T Follicular Cells and Contain Specific Subsets That Differentially Support Antibody Secretion

Immune Repertoire Sequencing Using Molecular Identifiers Enables Accurate Clonality Discovery and Clone Size Quantification

Towards Error-Free Profiling of Immune Repertoires

On A Measure of Divergence Between Two Statistical Populations Defined by Their Probability Distributions

Circlize Implements and Enhances Circular Visualization in R

Full-Length RNA-Seq From Single Cells Using Smart-Seq2

RSEM: Accurate Transcript Quantification From RNA-Seq Data With or Without a Reference Genome

Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data With Deseq2

Gene Set Variation Analysis for Microarray and RNA-Seq Data

An Algorithm for Fast Preranked Gene Set Enrichment Analysis Using Cumulative Statistic Calculation

Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets

NetworkAnalyst 3.0: A Visual Analytics Platform for Comprehensive Gene Expression Profiling and Meta-Analysis

Next-Generation Sequencing Reveals HIV-1-Mediated Suppression of T Cell Activation and RNA Processing and Regulation of Noncoding RNA Expression in a CD4 + T Cell Line

HISAT: A Fast Spliced Aligner With Low Memory Requirements

Camera: A Competitive Gene Set Test Accounting for Inter-Gene Correlation

Rapid Induction of Antigen-Specific CD4+ T Cells Is Associated With Coordinated Humoral and Cellular Immunity to SARS-CoV-2 mRNA Vaccination

Data Exploration, Quality Control and Testing in Single-Cell qPCR-Based Gene Expression Experiments

Intratumoral CD4+ T Cells Mediate Anti-Tumor Cytotoxicity in Human Bladder Cancer

Subspecialization of CXCR5+ T Cells: B Helper Activity Is Focused in a Germinal Center-Localized Subset of CXCR5+ T Cells

Human CD57+ Germinal Center-T Cells Are the Major Helpers for GC-B Cells and Induce Class Switch Recombination

Inhibition of Human T Cell Proliferation by CTLA-4 Utilizes CD80 and Requires CD25+ Regulatory T Cells

SMAD3 Prevents Graft-Versus-Host Disease by Restraining Th1 Differentiation and Granulocyte-Mediated Tissue Damage

Cytoplasmic Protein That Blocks Cell Cycle Progression From the G0/G1 to S Phase

The Role of IL-10 in Regulating Immunity to Persistent Viral Infections

ICOS Coreceptor Signaling Inactivates the Transcription Factor FOXO1 to Promote Tfh Cell Differentiation

Circulating CXCR3+ Tfh Cells Positively Correlate With Neutralizing Antibody Responses in HCV-Infected Patients

Follicular CD4 T Helper Cells as a Major HIV Reservoir Compartment: A Molecular Perspective

Tfh1 Cells in Germinal Centers During Chronic HIV/SIV Infection

Distinct Antibody and Memory B Cell Responses in SARSCoV-2 Naïve and Recovered Individuals Following mRNA Vaccination

CDC Guidance for Evaluating Health-Care Personnel for Hepatitis B Virus Protection and for Administering Postexposure Management

ZIKV Infection Induces Robust Th1-Like Tfh Cell and Long-Term Protective Antibody Responses in Immunocompetent Mice

Circulating Th1-Cell-Type Tfh Cells That Exhibit Impaired B Cell Help Are Preferentially Activated During Acute Malaria in Children

CXCR3+ T Follicular Helper Cells Induced by Co-Administration of RTS, S/AS01B and Viral-Vectored Vaccines are Associated With Reduced Immunogenicity and Efficacy Against Malaria

The authors would like to thank Dr. Jessica Podnar at the Genomic Sequencing and Analysis Facility at UT Austin for helping with transcriptome and TCR repertoire sequencing runs; Dr. Pengyu Ren at UT Austin for providing computational resources. We also thank the Cooperative Human Tissue Network (CHTN), a National Cancer Institute supported resource, for providing tissue samples.

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.859070/ full#supplementary-material Conflict of Interest: NJ is a Scientific Advisor and holds equity interest in ImmuDX, LLC, and Immune Arch, Inc. In the interest of full disclosure, DBW reports that he serves on Advisories for AstraZeneca, Geneos, Advaccinepharma, he participates in BOD service for Inovio. Remuneration received by DBW for these services includes SRA funding, direct payments, stock or stock options disclosed here.The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.